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The task has been completed of compiling colored plane mosaics corresponding to all 
15 symmetry groups, using elementary "tiles" with straight-line outlines only. 


In our basic article on this question [1],weused 
as an illustration of the group P4,; two four-color 
mosaics (I, II) which were formally correct, but not 
very successful. With the crystallographically 
equivalent two sorts of 4; axis present (as required 
by group theory), we found that at the point of emer- 
gence at one-half of each of these axes, eight ele- 
mentary triangular "tiles" converged, while for the 
other half there were only four. This was a natural 
result of choosing an elementary figure containing 
only one right angle. A shape with two rightangles 
is desirable. To achieve this we could use either a 
quadrilateral of the most general type (Fig.1la on 
the attached color table), or a rectangle (Fig. 1b on 
the same table). 

In the same article, we took an erroneous path 
from the two tetragonal mosaics with the centered 
group 14,md = F4,dm to the corresponding ortho- 
rhombic mosaics with the group Fdd; elongation and 
contraction of the cell must take place along d, the 
diamond planes, i.e., along the diagonal of the 
square shown in Fig.V, so that we then obtain Fig. 
VII and Fig. VIII in the form in which they are 
shown in the color table in the present article. 

Although mosaic IV is made up of strictly 
square tiles, the screw and colored 4; and 43 axes 
on it nonetheless only lie at the meeting places of 
four tiles, because (see Fig. 1) through the center 
of each tile there passes only an ordinary, two-fold 
rotation axis. Here, using the same square outline 
and with the same conditions imposed (4; and 43 axes 
at intersections and two-fold rotation axes at the 
centers of the tiles), we can obtain another kind of 
elementary figure by losing part of the area of each 


component square tile in the direction of one mirror 
plane and gaining the same area from two adjacent 
square tiles along the perpendicular mirror plane. 
The resulting mosaic (Fig. IVa in the coler table) 
has an elementary figure showing the re-entrant 
angles which are very common in Moresque decora- 
tive work [2]. This corresponds to the group I4;md= 
F4,dm, and an orthorhombic version is also pos- 
sible, with symmetry Fdd (Fig. VIIa on the color 
table). 

If each of the one-color squares in mosaic IV is 
left, in accordance with symmetry mm2, with only 
two triangles colored out of the four formed in each 
square by its two diagonals, then the empty tri- 
angles (see Fig.2) can be arranged in a different 
fashion, but with retention of the same mosaic sym- 
metry of I4,md = F4,dm. 

These can be colored in the same color as the 
triangles of adjacent squares, which gives mosaic 
V, again with square tiles, but with the areas of 
these only half those of the original tiles, and with 
their edges lying at 45° to the edges of the original 
square tiles. The centers of these smaller tiles do 
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not contain the two-fold axes, which now emerge at 
the intersections of four tiles, where they alternate 
with 4, and 43 axes (Fig. V). 

In the second version the empty triangles, which 
are equal in number to the retained colored tri- 
angles, are colored in four new colors so that the 
latter are related by exactly the same symmetry 
I4,;md = F4,dm, but remain completely independent 
of the basic triangles, and in particular are not re- 
lated to the latter through any normal or colored 
symmetry relationships. The possibility of simul- 
taneous existence in a lattice of given symmetry of 
two analogous, interpenetrating, but independent 
mosaics, each with a symmetry equal to that of the 
whole, has already been pointed out more than once 
[3]. 

In Fig. IVb the four main colors are the same 
as in the original Fig. IV, while the four additional 
colors will be white, orange, gray, and lilac. 

The drawbacks noted for the tetragonal mo - 
saics also apply to some of the mosaics given in [1] 
to illustrate hexagonal colored groups. In order for 
a mosaic of symmetry P38, not to have six tiles con- 
verging on each three-fold axis, as occurred in the 
formally correct Fig. X of [1], the corresponding 
elementary tile must be an n-sided polygon with 
three angles of 120° (cf. the tetragonal symmetry 
case), and this is only possible when n is not less 
than 5. Two mosaics made up of pentagonal tiles, 
both of symmetry P3,, are shown in Figs. Xa and 
Xb, and we see that only three tiles converge at 
the points of emergence of the three kinds of three- 
fold axes. 
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In Fig. XI, which should illustrate group R38, 
curved edges were wrongly given to the hexagonal 
colored tiles, and the mosaic took on the higher 
symmetry of R3m. In Fig. Xla, the correctly as- 
signed curved tile edges destroy the plane of sym- 
metry, and the colored group will, in fact, be only 
R38. This group can also be illustrated without 
curved edges, if we return to the pattern of mosaic 
XIII in [1], but use here only three colors instead 
of six. Two versions of such a mosaic are shown in 
Fig. XIb and XIc. A mosaic like XIc but in six 
colors would be another illustration of group R38c. 

The same arrangement as that in mosaic XIII 
can be used to illustrate the group P6,, if the tile 
outlines do not contain curved lines and if it is not 
necessary for six tiles to converge at the points of 
emergence of both the six-fold axis and of the two 
kinds of three-fold axes, as was the case in mosaic 
XIV. The corresponding two mosaics are shown in 
Figs. XIVa and XIVb. 
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A method has been devised for constructing simple three-dimensional mosaics. 


1956 


The 


simplest solutions are obtained when starting from the plane colored mosaics corre- 


sponding to the 15 colored symmetry groups. 


colored crystallographic group. 


When working with two colors, three-dimen- 
sional mosaics are derived most simply from the 
plane mosaics corresponding to the 46 two-color 
groups [1] by adding layers of identical plane mo- 
saics (of a certain thickness) in the third direction, 
but with the colors rearranged so that a unit "tile" 
of one color lies above a tile of the other color. 
This operation is repeated for each new layer. In 
other words, we introduce a vertical colored trans- 
lation. The corresponding two-color lattices were 
listed in [2]. 

Colored translations may, however, also be in- 
clined to the vertical (the layer direction), and this 
occurs in several simple mosaics derived from 
Fedorov systems of parallelohedra. 

If we use the usual set of nine heptaparallelo- 
hedra (Fedorov cubo-octahedra with 14 faces) to de- 
pict a body-centered lattice (i.e., its unit cell),and 
we make the middle cubo-octahedron a different 
color, we will obtain a very simple three-dimen- 
sional mosaic with a two-color cubic Bravais lat- 
tice (one of the 50 such lattices, see [2]), with the 
oblique colored translation P[Im3m. The unit cell 
of this lattice will contain one cubo-octahedron of 
each of the two colors. 

If we start from the simplest system of tri- 
parallelohedra, namely, close-packed normal cubes, 
and color alternate cubes along each of the coordi- 
nate directions in two different colors, then we ob- 
tain a three-dimensional chessboard pattern, which 
is oftenused to describe the structure of fluorite 
(CaF,) by saying that occupied cubes alternate with 
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A Klein V-group is illustrated by a 


empty ones (Fig.1); here, however, we will prefer 
to speak of the alternation of cubes of two different 
colors (but identical in all other respects). The cor- 
responding Bravais lattice will be cubic face-cen- 
tered, with all edges color-centered according to 
Fgm3m (cf. [2]). 

It is easy to see that this packing is a special 
case of the mosaics mentioned at the beginning of 
this article,i.e., it can be broken up into layers of 
identical plane mosaics with the two colors in each 
plane mosaic interchanged on passing from one 
layer to the next. In this particular three-dimen- 
sional mosaic, however, this operation of breaking 
up into layers is possible along three mutually per- 
pendicular directions. 

By applying uniform Fedorov deformations to 
the original (and thus two-color) systems of paral- 
lelohedra, we arrive at the corresponding three- 
dimensional mosaics of lower symmetry. 

Starting from the system of hexaparallelohedra 
used for the cubic face-centered lattice (rhombic 
dodecahedra), we immediately find that we cannot 
apply the simplest technique to obtain a three-di- 
mensional mosaic with the same cubic symmetry, 
because the unit cell contains four rhombic dodeca- 
hedra. If these are colored in levels along a four- 
fold axis, the symmetry is lowered to tetragonal. 
There is then no reason to have a base-centered 
lattice, and we go over to a primitive tetragonal 
lattice with P}4/ mmm color body-centering and 
with one rhombic dodecahedron of each color in the 
unit cell. 
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If the rhombic dodecahedra are colored in levels 
along a three-fold axis, this leads to the rhombo- 
hedral symmetry Ry3m (cf. [2]), with two rhombic 
dodecahedra of different colors in the primitive 
rhombohedral unit cell. 

The packing of rhombic dodecahedra does, how- 
ever, retain cubic symmetry if all four of the do- 
decahedra in the unit cell are given differentcolors. 
This first representative of a multi-color (here, 
four-color) three-dimensional Bravais lattice mo- 
saic will be primitive cubic, with inclined colored 
translations which center the different faces of the 
unit cell: this can be conventionally written (cf.[2]) 
in the form PA'B!c'. 

Using ordinary descriptive terms, this three- 
dimensional mosaic with rhombic dodecahedra of 
four colors can be described as a set of four "inter- 
penetrating" primitive lattices, all geometrically 
identical, but of different colors. All the rhombic 
dodecahedra in each of the four lattices are related 
by a system of translations of the primitive cubic 
lattice, translations which are noncolored or, more 
precisely, which do not change the colors of the 
rhombic dodecahedra. This set of translations will 
be denoted by T=1. Any lattice will be related to 
the other three through three-color diagonal trans- 
lations (along faces) of the primitive cube, denoted 
by TA, TR, and Ty (or, more exactly, by the set 
of these), i.e., these translations not only displace 
the rhombic dodecahedra by half a face-diagonal, 
but they simultaneously alter their color. 

In three-dimensional colored mosaics it is not 
possible to compare the coloring with a level which 
is higher by a certain fraction of the full repeat dis- 
tance [1], and the relationships between colors have 
a general character. In our case, the compatibility 
of these relationships is easily established by ex- 
amining Figs. 2a and 2b, which depict the unit cell 
of our mosaic. If the transition from the zero 
(black) lattice to the first (yellow) lattice occurs by 
means of the translation Tj, that to second (blue) 
lattice by the translation T4,, and that to the third 
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(red) lattice by the translation T¢, then it is not dif- 
ficult to see that the transition 1 — 0 (yellow to black) 
can be described by the translation Tae the transi- 
tion 1 — 2 (yellow to blue) by the translation TG@, and 


the transition 1— 3 (yellow to red) by the translation 
1 
Tp: 
If we consider all points (of the rhombic do- 
decahedron) in the same system (i.e., of the same 


color) to be equivalent, we find that 
=—-T2=TéZ=T =1, 
[4-Te=To, Te-To=La, To-Ts =f. 


ie = ry ete., 
From this it follows that Th ° TR “Taal; what- 
ever the order of the factors. 

These relationships join together the four trans- 
lations T, TA, Tp, and Té (or, rather, theircom- 
binations) into a mathematical group; the basic 
property, according to which the product of two 
colored translations is equal to a third, shows that 
in this combination of translations Pa rp tons we 
have another example of a mathematical V-group 
(Klein group [3], or fourth-order equation group), 
for which only three representatives were previous- 
ly known in crystallography [4], namely, the point 
groups V = D, ='222, Coy = mm2, and Cop = 2/m.? 


V =Dy Coy Con Pararcr 
2 
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diy 2, tty A ee @ G2 dg = Ay 
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As in [2], the primes mean that the corresponding translations, or 


symmetry elements in general, are colored. 


>This is formally not quite true, because the Klein V-groups also 
include the corresponding black—white groups [2]; 2'2'2, m'm'g, 
mm'2', 2'/m, 2/m', 2"/m'. However, the properties of these 
groups are extremely close to those of the noncolored crystallo- 
graphic groups, and so we will only mention them here. 


THREE-DIMENSIONAL MOSAICS WITH COLORED SYMMETRY 


It will be obvious that the groups of translations 
PA'B'C! will also be retained as such after uniform 
deformation of the system of hexaparallelohedra or 
rhombic dodecahedra. A more important point is 
thatthe system formally loses its cubic symmetry 
even earlier, as a result of the imposition of the 
four colors. Although the colored figure in the cen- 
ter of face A can be superimposed in turn on the 
centers of faces B and C, and then on itself, where 
it regains its color, this is brought about through 
the sequence of different operations As chal TO = 
1, and not by repeating the same operation three 
times, e.g., a left-handed (or right-handed) rota- 
tion about a three-fold axis. These latter axes are 
absent from the system of colored rhombic dodeca- 
hedra, and it is these four three-fold axes which 
would determine that a system of figures belonged 
to the cubic system. It is easy to see that the spa- 
tial symmetry retained by our packing will only be 
the orthorhombic PA 'B'c mmm with the three unit- 
cellconstants the same only "by chance." Witha rhom- 
bically deformed cell the three-fold axis will evi- 
dently not be present, but the face centers will 
translate one into the other as before through the 
colored translations, so that a fourth translation 


will bring a face center back to its original position. 


We will arrive ata particularly simple pseudo- 
cubic system of figures, withthe same symmetry of 
Pa'rp'c mmm, if we start off from the fluorite struc- 
ture examined above when it is represented in the 
form of Pauling cubes (Fig.1), and apply three new 
colors to the "occupied" cubes at the centers of all 
the faces of the unit cube. In actual fact this willbe 
the same mosaic as that which we have just exam- 
ined, because a rhombic dodecahedron is the same 
as a cube with a pyramid added on each face, where 
the pyramid is obtained by cutting the cube in six 
along its diagonal symmetry planes (Fig.3). If then 
these pyramidal attachments are taken off, it is 
easy to see that spaces will be left between the 
cubes which will also be cubes of identical size,and, 
as we have repeatedly emphasized, the symmetry 
of this identical system of empty cubes will be the 
same as the symmetry of the system of occupied 
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cubes, and also the same as that of the combination 
cf both cube systems. 

In the symmetry groups of those mosaics which 
we constructed from parallelohedron systems, the 
basic colored operation was a second-order trans- 
lation. In order to obtain groups with coloredopera- 
tions of third and higher order, it would be sim- 
plest to revert to the 15 multicolored plane groups, 
and their corresponding mosaics, which were given 
in [1,5]. If we consider the latter to actually be 
plane mosaics of a certain thickness and possess - 
ing colored elements, in particular colored sym- 
metry axes, we can extend these upward by adding 
layers of similar mosaics but with the colors al- 
tered by a colored translation, in the same order 
as that which the colors follow around the colored 
axes (or along glide planes, particularly the d- 
planes). Thus, if we use the four tiles in mosaics 
I, Il, Ia, and Ib (from the color tables in [1,5]), 
which appear in the order yellow (y), green (g), 
red (r), and blue (b) around the four-fold colored 
axis, we will add on four tiles with the color order 
g-r-b-y, and then the four r-b-y-g and b-y-g-r, 
finishing up with the original order of y-g-r-b. The 
colored vertical translations are fourth order. 
Operated once, they raise a color by 90°, twice by 
180°, then by 270°, and so on: y.TIy=g, y.Tty = 
Tee, y-Thy = yT yh, =hb, y. Thy = We 8-Tiy — wets 

If we obtain a plane mosaic by starting with a 
4, axis and noting the change in color with each 
raise along this axis, we then find that our mosaic 
has a 4go9] axis, and because of the layering of the 
tiles through the action of the colored translations 
we restore the ordinary one-color 4, axis, whichis 
fixed on the four tile levels by the four differently 
colored spirals (Fig.4). 

The symmetry group of the colored mosaic 
formed may be denoted by either Pe'/44col1 or 
Pe1/44,4. In the first notation the derived symmetry 
element will be 4; = 4e9]°¢'/4. We also have it that 
4? = 2; = 4¢9]‘c'!/2. With the second notation the 
derived symmetry elements will be 


Aoi 41° (7/4) ao Zoo 21° (07/4)? 
Similarly, from the mosaics belonging to other 
colored symmetry groups depicted in [1], we must 
be able to find in these the multicolored symmetry 
elements 60]; Ope nite 3e01°2; 3e9]1, Ago] and the 
plane three-color lattice He9] which correspond to 
the lattice R. By piling up layers of identical mo- 
saics, using colored translations of the same order 
as that of the major symmetry element of the plane 
mosaic, we "reconstruct" the one-color symmetry 
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elements 6;, 6), 3,;, d, and R from which we started 
off in [1] to derive the colored plane symmetry 
groups, and which now will be realized in the form 
of the corresponding packings of each color. 

We note that the three-dimensional mosaic 
Re'/33m, derived from the plane mosaic R3m (XII 
in [1]), gives the spatial pattern of ilmenite, FeTiO;, 
if we take it that the empty octahedron in this struc- 
ture plays the same crystallographic role as the oc- 
tahedra with Fe and Ti [6]. 

Although the mosaic XIII has the symmetry 
group R8c, and six colors, the R(H) lattice in the 
plane pattern corresponds to a sequence of only 
three colors, and the corresponding vertical trans- 
lation will break down into three colored transla- 
tions 120° apart, i.e., along the vertical we will 
have either black (bk)—red (r)—blue (b) or white (w) 
—yellow (y)—green (g) alternating, in accordance 
with our accepted cyclic color sequence w—b—y-— 
bk—g—r-—w (Fig. 5). 

We can form more three-dimensional mosaics 
by making the colored translation 180° instead of 
90° for a four-fold colored axis, and by having 
translations of 120° and 180° as well as 60° witha 
six-fold axis. The resulting (one-color) groups will 
respectively replace P4,(43) by the monoclinic 
P2,/m, and P6,(6;) by P2,/m again (180°) and P3, 
(120°). Higher symmetry will only occur with four- 
fold and six-fold colored axes. 

We cannot use the above principle to obtain 
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three-dimensional mosaics from plane mosaics 
with a three-fold colored symmetry axis (in [1,5], 
the mosaics X—P3,, XI—R3, XII-R3m, XIII-R3c, 
XV-P6,). A horizontal mirror plane can be intro- 
duced into the appropriate three-dimensional mo- 
saic by other means, namely, by the use of further 
colors (by raising each of the three basic colors by 
180°). This is most simply achieved with mosaic 
XIII, where the additional colors are already pres- 
ent [1]. On the layer depicted in mosaic XIII we 
stack a similar layer with the colors interchanged 
on the complementary principle, and this is re- 
peated for each new layer. In the three-dimension- 
al mosaic the colored translation c{/, appears, but 
the lattice R will be absent. The retained glide 
planes plus the horizontal mirror plane generate a 
horizontal two-fold axis, i.e., the noncolored sym- 
metry will be P3/mc2 = P6c2 = Dp. 

To introduce a horizontal mirror plane into the 
three-dimensional mosaic constructed from mosaic 
XII, which has the symmetry group R3m, we alter- 
nate layers of the latter with layers of a similar 
mosaic made up with colors raised through 180°, 
according to mosaics XIII and XIV. The space (non- 
colored) group of the resulting three-dimensional 
mosaic will be P3/mm2 = P6m2 = Dip. 

In the same way, an m plane present with a 
perpendicular colored translation (180°) can be in- 
troduced into the three-dimensional mosaic based 
on mosaic XV, with colored group P6,. The re- 
sulting noncolored (space) group will, however, 
only be the monoclinic group P2/m. 
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It is now reasonably well established that the 
theoretical curves for electron scattering derived 
from the Born approximation agree well with ex- 
periment, at least for atoms of low and medium 
atomic number. This is confirmed by recent ex- 
periments [1,2], which also gave the deviations for 
the heaviest element used (gold). It has been pointed 
out [8] that the Born approximation is inadequate 
for heavy atoms. The establishment of the scatter- 
ing law applicable to a crystal in essence consists 
of two stages: 1) type identification (kinematic or 
dynamic), and 2) derivation of the atomic scatter- 
ing curve. The first aspect must be settled before 
the second can be discussed; but the scattering type 
(critical thickness) is itself dependent on the atom- 
ic amplitude [4], which makes the two aspects re- 
lated. This is not surprising, because the physical 
basis of the Born approximation is neglect of sec- 
ondary waves within the atom, the secondary waves 
throughout the crystal being neglected in the kine- 
matic approximation [5]. 

The Born approximation is reliable for light 
atoms, so their fe curves are of interest in rela- 
tion to bond effects. These fe curves are much 
more sensitive to changes in the outer parts of the 
electron shells than are the x-ray f curves, because 


Z—fe 
jeter (sin 97)? (1) 


at small angles, where f approaches Z (the nuclear 
charge), the electron scattering becomes essential- 
ly a differential effect; and it is a general property 
of the f curves that their behavior at small angles 
is governed mainly by the charge and potential of 
the outer parts of the atom. The peak heights and 
shapes for Si and F in (NH,),SiF, are [6] such as to 
show that Si has a positive charge and F a negative 
one. The shape of the fe curves may also be exam- 
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ined directly without Fourier synthesis of the poten- 
tial. 

Li,O is very convenient here, for it has aface- 
centered structure of antifluorite type; the struc- 
ture amplitude is finite for unmixed indices and 
takes the form 

feot2feri for 


h+k+l=4n (2a) 


feo for hh EES et eb) 


le o— 2fe, 11 for 


One of the species of atom (Li) thus takes no part 
in the scattering if the reflection indices are odd, 
and this gives fe © directly. Again, this curve may 
be deduced from half the sum of the curves drawn 
through the Pak that satisfy (2a) and (2c), the com- 
bination (2a)-(2c) then giving fe, Li- 

The LigO was made by distilling lithium metal 
onto collodion films at about:10+ mm Hg in the dif- 
fraction apparatus; this pressure is such as to pro- 
duce complete oxidation in nearly all cases (the 
lines of metallic Li are only occasionally seen 
among those-of Li,O). Multiple exposures and 
microphotometry were used as in [6,7] to deduce 
the intensities. The specimens were usually tex- 
tured, but the measurements were based only on 
polycrystalline Li,O showing no signs of texture,as 
deduced from oblique-incidence patterns [1]. Very 
weak reflections gave no peaks on the microphotom- 
eter trace, so these were evaluated visually. The 
intensities Ip,k7 were converted to the ®p,] via 


Dnt = V Inne / Cher p- 


Averages over four specimens were taken; the I 
were measured to about 10% so the @ were found to 
about 5%. 


Pari = 


htk+l=4n+2 (2c) 
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TABLE 1. Eaperimental &),) for Li,O (a = 4.60 A) 
hkl © hkl © hkl © 
111 10.6 420 Ao) 933 ned 
200 —1.3 422 EE) 444 2.6 
220 16.0 | 333. 541 3.2 711, 554 1.5 
Sule a? 440 3,7 64 2.0 
raped 1.4 531 Due d03. 731 1.4 
400 9.7 | 600, 442 0.5 
334 4,7 620 aul 


The lines were broad and indicated a small 
crystal size, while the fe curves fell even more 
rapidly than the theoretical ones (see Fig. 2; thisis 
due to the temperature factor); the two features to- 
gether point to kinematic scattering. Table 1 gives 
the relative values of the ®);, ; Fig. 1 gives the re- 
sulting curves. 
These f, curves must be normalized in some way 
for comparison with experiment; it is convenient to 
do this for sin¥’/A small, because the temperature 
factor here has only a small effect; but this is pre- 
cisely the region where the other effects on fe are 
pronounced. For this reason, we neglect the tem- 
perature factor temporarily and normalize the 


All except 99) are positive (Fig.1). 
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Relative values of @pp7 for Li,O. 


curves by equating f © to the theoretical value for 
sin 3/X = 0.3-108 cm, the theoretical value being 
taken from the tables of [8], which themselves are 
based on MceWeeny's f [9] for light elements. This 
sin ’/XA represents roughly the point at which the 
f (and hence the f,) found in different ways come 
into agreement. Figure 2 shows the experimental 
and theoretical fe curves calculated from (1) in ac- 
cordance with [8,9] and also with [10] for free O and 
Li atoms, and also for O?" and Li’. Table 2 gives 
the f, for the elements as normalized in this way 
(without allowance for the thermal motion). 

Figure 2 compares theory with experiment; 
small sin’/A are of special interest. The features 
of (1) are such as to produce very marked discrep- 
ancies between the theoretical curves; the experi- 
mental vie (curve 1) is clearly closer to that 
given by [8] (curve 2), which is based on the data of 
[9], than it is to that given by [10] (curve 3). This 
shows that McWeeny's basic assumption (the need 
to allow for bonding in the lattice) gives (as would 
be expected) closer agreement with experiment than 
does the calculation of [10], which was based on 
Hartree's method for free (unbonded) atoms. But 
curve 1 for fe,o runs below curve 3 at small 
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TABLE 2, Experimental f ¢* 


(sin $/A)-10-8 0.45 | 0.2 0.29 }o.3 | o- [0.0 0.6 0,7 Jos 
feo 0.55 | 0.51 | 0,47 | 0.425 | 0.315 | 0.215 | 0.135 | 0,085 | 0.065 
fe zi | 0.30 | 0.24 | 0.195 | 0.120 | 0,070 | 0.040 | 0.029 | 0.022 


* The factor required to convert to absolute fe is 2.393: 1078 cm. 


sin’/A , because Li,O does not have a purely co- 
valent lattice. Curve 4 gives fe, o- as calculated 
from (1) for the free anion from the data of [10]. 
The lower level of curve 1 is thus to be ascribed to 
a certain proportion of ionic bonding. Our results 
show that the common view (that the bonds in com- 
pounds of Li,O type are mainly ionic) is wrong, for 
here the proportion of ionic character is low. It is 
difficult to draw more detailed conclusions because 
Li,O has a = 4.60 A only, so the experimental f, 
can be found only from sin ¥/A © 0.15 + 10® cm. 
The theoretical results [9,10] for Li (curve 5) 
agree extremely closely; ionization effects appear 


only around sin ¥/d of about 0.1 - 10° cm (curve 6), 


0.50 
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Fig. 2. Theoretical and normalized experimental f¢ curves; 
the factor required to convert to absolute f¢ is 2.393 - 10° cm: 
1) oxygen, experiment; 2) theory [8,9] for bound oxygen; 3) free 
oxygen [10]; 4) free O?" [10]; 5) lithium [8-10]; 6) Li* [10]; 

7) lithium, experiment. 


which is inaccessible to measurement for Li,O.This 
makes it difficult to draw any conclusions for Li, 
although the experimental Fe, Li curve 7 does devi- 
ate somewhat from the theory, especially in that it 
does not show the predicted sharp upward trend. 
The experimental fe curves in the range 
sin ¥/A = 0.3-0.8 - 10° cm allow us to select suit- 
able temperature factors for the theoretical f ¢ 


curves, the Bin e-# einen" being Bo ~ 1.0 A? and 
ST > 80 A*. The former is the larger because, 
for remote regions of sin 8/2 we should have 

f e,0/f e, Li close to ZO/Zyj = %% = 2.66, if we 
neglect the thermal motion, whereas experiment 
gives it as about 2.9. 
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A direct method (statistical equations) is shown to be applicable to a compound charac- 


terized by 60 atomic parameters. 


Direct methods are much used in chemical 
crystallography, especially the statistical method 
for the signs, which is used in various forms. 

Zachariasen [1], in 1952, proposed a relation 
between the signs of the amplitudes Fy, Fx, and 


FH+K! 


Sy = S(Sx,; > Su+x;, ); (1) 


which he applied to metaboric acid; it gave the signs 


of 13% of the structure amplitudes (the largestones). 


Kitaigorodskii [2] gave the theoretical basis for this 
formula in terms of the structural product of the H 
and K indices: 

Xun = Fu FP x Fuse. (2) 


Belov and his colleagues have developed several di- 
rect methods for the signs and have applied these 
to minerals [3]. 

We have used another form of the method [4] 
via the formula 


Dificos (Hj + Kj) + Dif; 008 (Hj — K;) 
j 
= 2 D\f;cos H; cos K;, (3) 
) 
in which 
Hy = 20 (hx; + ky; + 12;), 

which relates the signs of Fy+K, Fy—K:; FH, and 
Fx. This formula becomes, in the particular case, 


cos Kj =1 


Srix:|Fusn| + Su—x|Fu—-x|=2-Sx:Sq:|Fu|. (4) 
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This shows that only the signs of two ampli- 
tudes appear here, as against the three of (1), and 
only two amplitudes, namely,Fy and Fx, so (4) may 
be put in general form as Sy 4K =SySx OF Sy_K = 
Sy‘ SK, provided that the moduli of the structure 
amplitudes satisfy 


|Fu—x| >|Fu+x|- 
Equation (4) can be used only via a Statistical survey, 
except where there is a heavy atom ina special posi- 
tion, because any deviation of the cos Kj from unity gives 
rise to errors in the signs; statistical averaging of 
the signs then involves as necessary and sufficient 
condition that the modulus of only one structure am- 
plitude in the product IF should be large. It is of 
interest to apply the method implied by (8) to com- 
plex structures; we first used it to advantage in joint 
work with G.S. Zhdanov and K.I. Tobelko on the 
structure of realgar, As,S, [5]. The method was sub- 
sequently applied to the coordination compound 

F3B «-NC;H,;, which contains only light atoms. The 
x-ray results on this are @ = 17.71, b = 5.89, c= 


|\Puysx| >| Pax) or 


Fig. 1. Patterson projection on the (010) plane. 
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Fig. 2. Electron-density projection on(010) plane ;signs of the F(hOZ ) 
determined directly from statistical equations. 


14.34 A; B = 118°40'; op = 1.48 g-cm™, oy =1.49 
g-cm°, z= 8; space group C3), = P2,/c. 

The general fourfold position 4(e) in group 
P2,/c is taken twice for z = 8. Figure 1 shows the 
Patterson projection first calculated; this revealed 
the layered structure, but the outline of the mole- 
cule was not seen. The signs of the F(h0/) were 
then deduced via the statistical equations on the 
basis of 112 reflections of h0l type, four equations 
being drawn up for each amplitude. The statistical 
averaging was based on structure amplitudes giving 
moduli [Fx] greater than 0.25. 

Figure 2 shows the resulting electron-density 
projection, which reveals all of the atoms. The mu- 
tual displacement of the molecules explains the 
largest peak in Fig. 1; the calculated intensities al- 
so confirmed the configuration. The signs were re- 
fined successively by reference to the coordinates 
of the atoms; Fig. 3 shows the resulting electron 
density projection. The final coordinates showed 
that the direct method gave errors in the coordi- 
nates less than 0.1 A, and that 75% of the signs had 
been determined directly; about 15% of the signs for 
weak reflections were undefined, while 10% of the 
F(h0/) had been given incorrect signs. This led to 
distortion of the shape and size of the electron- 
density peaks, so the heights of these at first did 
not correspond to the numbers of electrons on the 


Li-N« LVONKOVA 


Fig. 3. Electron-density projection on the (010) after successive ap- 
proximations. 


atoms. Electron-density series based on the cor- 
rect signs gave electron numbers approximately 
proportional to the peak heights (see [6]). A point 
here is that for NaCl, with Zy,+= 10 and Zq]- = 
18, the correction for deviation from proportion- 
ality is '/, of an electron for ZNat> Which lies at the 
limit of detection for x-ray methods. The coordi- 
nates were further refined from three-dimensional 
series. The statistical method is thus applicable to 
a coordination compound characterized by 60 atom- 
ic parameters. 
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Following x-ray diffraction studies of single crystals, the crystal structures of Th, Znj7 
and U,Znji7, constituting new structural types, have been completely determined. 


Introduction 


In 1946 there appeared a paper by Nowotny [1] 
indicating that the thorium—zine system contained 
a compound ThZny with hexagonal structure and lat- 
tice spacings 


a = 9.237 kX, c = 4.442 kX. 


As regards positioning of the atoms in ThZng, 
Nowotny found that the structure constituted a modi- 
fication of the CaZn; type, in which the place of 
calcium was occupied by a "statistical atom" 

(Tho, ¢Zno4), while the zinc sites were fully occu- 
pied by zinc atoms. Hence, the composition of 
ThZn, should be given the formula (Thy ¢Znp 4)Zn;. 

It should be noted that Nowotny prepared his al- 
loys in open crucibles, using flux, which apparent- 
ly led to severe oxidation and contamination. He ob- 
served that the rotation x-ray photographs contained 
not only layer lines corresponding to the spacing 
4.442 kX along the z axis, but also other layer lines 
which he ascribed (without adequate grounds, inour 
cpinion) to impurities of thorium dioxide. 

In 1953 we began an x-ray structural study of 
single crystals and powders of a uranium —zinc 
compound with composition initially taken as cor- 
responding to the formula UZng (in analogy with 
ThZn,). Further study, however, showed that the 
real compositions of these compounds corresponded 
to the formulas Th, Zn; and U;,Znj7, and their crys- 
tal structure, although quite akin to that of CaZn,;, 
was considerably more complex and belonged to a 
new structural type. 

In this paper we present results of our solution 
of the crystal structures of Th)Zn,7 and U)Znj7. In 
the initial stages of the work it became clear that 
the two structures were very similar and that so- 
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lution of U,Zn;7 should be preceded by solution of 
the simpler structure of Th)Zn,7. By combining 
harmonic analysis with the method of trial and er- 
ror, we were able to make a complete determina- 
tion of the atomic positions in crystals of Th) Zn,7 
and U,Zn47. 


Preparation of the Alloys 


For preparing the alloys we used metallic 
uranium of 98.95% purity produced in the Soviet 
Union, metallic thorium of 100% purity from 
Germany, and 99.95% pure zine from our own in- 
dustry. By roughly filing monolithic pieces of these 
metals, we obtained powders; after taking the ap- 
propriate proportions of the components these were 
pressed into 12-g tablets. The tablets were placed 
in a beryllium -oxide crucible provided withacover, 
which was then inserted into a thin-walled quartz 
ampoule. The ampoule was evacuated to the region 
of 10 mm Hg and after sealing placed in a verti- 
cal tubular electric furnace. The temperature was 
raised gradually, holding at 400°C for 1 h, 800°C 
for 1 h, and at the maximum for 10 min; this was 
1000°C for uranium alloys and 1250°C for thorium 
alloys. There was no bursting of the ampoule. In 
order to guard against sublimation of the zinc, the 
ampoules were made as short as possible and 
placed in a constant-temperature zone of the fur- 
nace. The original mixtures corresponded to the 
formulas ThZng and UZng. After fusion the furnace 
was cooled slowly, passing from 1200 to 800°C in 
3 lal 

The alloys had homogeneous micro- and macro- 
structure, consisting of fairly large single crystals 
grown together. In the uranium alloys the single 
crystals had a lamellar appearance and were formed 
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as stacks growing together along the base. In the 
thorium alloys the single crystals had the form of 
fragments showing various faces. 

We selected the finest fragments with well- 
formed faces and 0.1 to 0.3 mm in size for x-ray 
examination. 

The single crystals of the uranium compound 
were extremely stable in air and retained the luster 
of their faces for two years. The single crystals of 
the thorium compound were less stable and tar- 
nished in air after a few months. 


Chemical Composition of the Single 
Crystals 


For an exact determination of the chemical 
composition of the compounds, we made a chemical 
analysis of single crystals selected by means of a 
lens. The weights of the single crystals were 0.5 
tones 

The method of determining the uranium was 
based on its reduction from the hexavalent state in- 
to the tetravalent by means of zinc amalgam with 
subsequent titration of the solution with potassium 
permanganate. From two analyses with 0.53 and 
1.04 g of the specimen we determined uranium con- 
tents of 29.5 and 29.9 wt.%, The mean uranium 
content of the uranium—zine compound was thus 
29.7 wt.%, corresponding to 10.4 at.%, or the for- 
mula UZng.¢, i.e., on rounding off to whole num- 
bers, U,Zny7. 

The thorium was determined as ThO,. From 
three parallel analyses of 0.5 to 0.7 g we deter- 
mined percentage contents of 31.00, 31.15, and 
31.19 wt.% Th, giving a mean value of 31.11 wt.%. 
This gave the formula ThZn; 5, close to ThZng. 

Thus, chemical-analysis data show that the 
single crystals of the uranium compound had com- 
position U,Zn,7 and those of the thorium compound 
ThZng. However, as we shall show below, the com- 
position of the latter compound is really nearer to 
Th) Znj7 than ThZng. We note that the structural 
data also allow the compositions ThZng and UZng as 
a result of a variation in the number of zinc atoms 
in the unit cell of the variable-composition phases 
Thy Zny, 47 and UpZny,_.47. In thiscase, it is more 
or less immaterial which composition we take as 
fundamental for determining the structure, ABs or 
Ay Bi7- 

As initial composition of the compounds we took 
ThyZny7 and UyZny7, which was later confirmed by 
the number of atoms in the cell obtained from ex- 
perimental densities. 
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(a) Th,)Zn,7. Rotation x-ray photographs 
with a Th)Zn,7 single crystal rotating around the c 
and a axes were obtained in unfiltered cobalt radia- 
tion, using an RKV-86 camera. The c-rotation 
photograph gave six layer lines and the a photo- 
graph gave four. The symmetry of the Laue and 
oscillation photographs showed that the Th) Znj,7 
crystals belonged to the hexagonal system. Meas- 
urement of these photographs gave the following 
lattice spacings: 


a = 9.03 A, c = 13.20 A. 


Thus,our spacings are, respectively, V3 and3 
times larger than those proposedby Nowotny [1]. 

The experimental density of Th,Zn,7 single crys- 
tals determined pycnometrically in benzene was 
8.32 g/em®. Hence, the number of atoms per cell 
becomes 56.3. If we round this off to 57, the unit 
cell contains three formula units of Th)Zn;,;. The 
x-ray density of Thy Zn; will be equal to 8.39 g per 
cm?, 

(b) U,Zn,z7. Rotation x-ray photographs of 
U,Zn,; single crystals taken around the c and @ axes 
were obtained in filtered cobalt radiation, using the 
RKV-86 camera. 

From the symmetry of the Laue and oscillation 
photographs, hexagonal symmetry was assigned to 
the U,Zn,7 crystals. 

Measurement of the x-ray photographs gave 
the periods: 


a = 8.99 A, c = 26.35 A. 

We note that the a spacings of Th)Zn,7 and 
U,Zn,7 are approximately the same, while the c 
spacings of U,Zn;,7 are twice those of Thy Znj7. 

The experimental density value of the U,Znj7 
single crystals was 8.50 g/cm’. 

This gave 113.1 atoms in the unit cell. 

This number agrees closely with the figure of 
six formula units of U,Zn,7 in the cell, in all 114 
atoms. On this basis the x-ray density of U,Zn,, 
equals 8.57 g/cm’. 


Determination of the Space Group 


In order to determine the reflection intensities 
and to effect indexing, a series of oscillation photo- 
graphs was obtained in the RKV-86 camera around 
the a and c axes, the oscillation range being 15°. 
The relative intensities were determined by stand- 
ard blackening marks. 

(a) Th,Zn,7. On considering the indices 
of the spots on a series on oscillation photographs 
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of TheZn,7, it became clear that all satisfied the 
rhombohedral extinction conditions; this yields the 
following possible space groups: D3q = R3m, Ds = 
R32, C3y = R3m, Ci = R3, Ch= RB. 

‘From the symmetry of the Laue and oscilla- 
tion photographs, the last two groups drop out, 
Since the diffraction class is determined as Dsg. 
Since we had no experimental indications suggest- 
ing noncentrosymmetric groups, we chose the cen- 
trosymmetric D3qg—R3m. 

(b) U,Znj,7. It follows from the symmetry 
of the Laue and oscillation photographs that U,Zn4; 
belongs to the diffraction class Dg,. Among in- 
dices of the 00/ type, only those satisfying the con- 
dition 1 = 3n are present. From this it would be 
possible to select x-ray group No. 93 [3], includ- 
ing the space groups Dj = C6,2 and D2 = Cé,2. It is 
also possible, however, that spots with indices 001, 
where l/ ~ 3n, were absent because of low intensity. 
This supposition found support in later study which 
showed that it was more correct to state that there 
were no systematic extinctions in this case. Thus, 
we have the following space groups of class Dg}, to 
consider: Diy = C6/mmm, D} = C62, Ch, = Cémm, 
Di, = C6m2, and D’y = C62m. 

Subsequently, on the basis of the close struc- 
tural kinship of U,Zn,; and Th,Zn,;, we found that 
a satisfactory disposition of atoms in regular sys- 
tems of points could only be set up in the space 
group Din = Cém2. 


Determination of Atomic Coordinates 


The rotation x-ray photographs of Th) Zn,7 and 
U,Zny7 around corresponding axes are very similar 
to one another both in the dispositions and intensi- 
ties of the spots. The powder photographs of the 
two compounds are even more alike. These experi- 
mental facts indicate that the structures of the com- 
pounds are also very similar. 

In the unit cell of Th) Zn,7 we have to site 6 Th 
atoms and 51 Zn. In the unit cell of U,Znj; we must 
site 12 U atoms and 102 Zn. From the similar in- 
tensities of corresponding reflections in the two 
compounds we may suppose that the heavy metal 
atoms are analogously situated in both structures, 
and that the doubling of the c spacing of U,Zn,7 as 
compared with Th)Zn;7 is due to the appearance of 
a horizontal mirror plane on passing from the 
rhombohedral group Dq = R3m to the hexagonal 
group Din = C6m2. 

The procedure for determining the coordinates 
of the atoms in the two cases was in general re- 
spects as follows. First, we constructed Patterson- 
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Harker cross sections along the basal planes (using 
101 reflections) and diagonal plane (1120) (using 285 
reflections) for U,Zny7, with 7.5° and 3° strips. The 
interatomic vectors U—U could clearly be deter- 
mined on these sections; there were no clear indi- 
cations of the U—Zn and Zn—-Zn vectors. 

From the sections referred to the basal plane 
it was concluded that the horizontal coordinates of 
the heavy atoms had values 0,0,'4,°4, and %,%. 
From the diagonal section it was established that 
the z coordinates of the uranium atoms could have 
values’ 2 0F7, 45, eae ae 

From considerations of the analogous disposi- 
tion of heavy atoms in U,Zn,7 and Th) Zn;7, we as- 
sumed that the thorium atoms in Th)Zn,7 could only 
have the following values of z coordinates: 0,4/3,7/. 
Only the six-fold position 6(c) with parameter z = 
tf, satisfies this condition in the space group Dig = 
R3m. In this case we may suppose that the ar- 
rangement of uranium atoms in one half of the 
U,Zn,7 cell will be similar to that of the thorium 
atoms in the Thy)Zn,7, while the other half of the 
U,Zn,7 cell will be the mirror image of the first 
half. These relationships are illustrated in Fig.1, 
which shows the cross section of the two cells in 
the (1120) plane. 

The coordinates of the majority of the zinc 
atoms in the two structures were established by 
trial and error on the basis of a model with the 
CaZn; structure. 

As seen from the above relationships for the 
cell dimensions, the volume of our Th, Zn,7 cell is 
nine times that of the (Tho ¢Zny 4)Zn; cell of the 
CaZns type proposed by Nowotny [1]. 

Structure of the CaZns type belongs to the space 
group Di, = C6/mmm. The Ca and Zn atoms in this 
structure have the following coordinates: 1 Ca in 
1(a) :000; 2 Znz in 2(c): 4,74, 0, %%4, 44,0; 3 Znqy in 
3(g): ve 0,%, 0,45, h; We %, 


U,Zny, Th, Zn, 
Fig. 1. Cross section of the structure of ThpZny7 


and Up,Zm,7 along the diagonal plane (1120). 
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9Zn in 9(d), 
: 1 
18Znin 18(/), c= 3, 
18 Zncin A8i(h) eect eee 
nin 18(h), c= >,2=%, 
6 Zn in 6 (c)yzz, 2 = 0.097. 
(b) U,Zny7- Space group Di, = Cém2.Co- 


@-Thu O-2Zn, 
Projection of ThyZny7z and U,Znyy 
structure on the basal plane. 


@-Zn, 
Fig. 2. 


We started from the assumption that the zinc 
atoms in ThoZnj7 and U,Zny7 were arranged accord- 
ing to the CaZn; motif as in Fig.2, which gives a 
projection of the CaZn, structure (cell indicated by 
broken lines) on the basal plane of the ThyZn,7cell. 
The structural motif of the CaZn; type gives us po- 
sitions for 45 zine and 9 calcium atoms, six of 
which are occupied by thorium atoms in the Thy Znj,_ 
cell, as shown above. In the three remaining posi- 
tions not occupied by thorium we must have the re- 
maining six atoms of zinc. Geometrical analysis 
showed that these six zine atoms could only occupy 
positions on certain vertical axes containing thori- 
um atoms. A Harker section with respect to the 
hexagonal axis confirmed this conclusion and led to 
the value z = 0.097 for ThyZn47. 

From the proposed arrangement of the zinc 
atoms and the coordinates found for the thorium 
atoms, we determined the signs of the structure 
amplitudes for the 0k! and hk0 zones and con- 
structed electron-density projections on the basal 
plane and the diagonal plane (1120). The synthesis 
was made in an orthohexagonal cell from 56 inde- 
pendent Okl reflections and 11 reflections of the 
hk0 type. In view of the considerable truncation of 
the series, good agreement was only found for the 
thorium atoms in the electron-density synthesis. 
Very large maxima on both electron-density charts 
confirmed the coordinates found for the thorium 
atoms. 

Further study of intensities showed that our 
motif for the arrangement of the zinc atoms on the 
CaZn;-type structure with six zinc atoms in posi- 
tion 6(c) with parameter z = 0.097 was correct. 
Thus, all these experimental data lead us to the 
following arrangement of atoms in the Th,Zn,; and 
U,Zn47 cells: 

(a) Th» Z Ni7- 
ordinates of the atoms: 


Space group Dig = R38m.Co- 


1 


6Thin 6(c), z= 3, 


ordinates of the atoms: 


{Und ce, 2. ins2ie eee 


7 
1U in 1(d), 2U in 2(g)\n, 2= 5, 
1U in4(e), 2U in (hyn z=, 
1U int(f), 2U in 2(i, 2=4 

6 Zn in 6 (i), t=a, y =0, 

6 Znin 6(m), z=+, y=0, 
12Zn int2(o, 2=+, y=0, 2=5 


12Zn in 12(0)1, c=z, y=0, sas 


6 Zn in6(n), ©= 5, 2=5) 
6Zn in 6(n)n, c=, z=, 
6Zn in 6 (nym, =+, zane, 
6Znin 6(n)v, =, 2= a5 
6 Zn in 6 (ny, =, 2=4, 
6 Zn in 6(n)v1, t=, La, 
6 Zn in 6 (n)yu, c= 2, c= 
6Zn in6 (nym, c= , s=4, 
6Zn in 6 (nix, =>, t= pe 


2Zn in6(g);, z = 0.049, 
2 Zu in 2(g)y, 2 = 0.451, 
2Zn in 2(h)y1, 2 = 0.118, 
2Zn in 2(hk)yq, 2 = 0.784, 

2Zn in 2(i)yz, 2 = 0.284, 
2Zn in 2(i)yy, z= 0.618. 


Figure 3 gives a comparison between the theo- 
retically calculated and experimental values of F* 
for the zero layer line of oscillation x-ray photo- 
graphs of ThyZn,; around the c axis. Here we have 
used averaged values of the experimental intensi- 
ties for spots of a given hkl family. 

Figures 4 and 5 give a comparison between the 
experimental intensities of the leading lines on pow- 
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Fig. 3. Comparison between theoretical and ex- 
perimental values of FeO for the oscillation x- 
ray photographs of a Th,Zm, single crystal ro- 
tated around the c axis. 


der photographs of Th, Zn;7 and U,Zn,; estimated 
visually from blackening standards and those cal- 
culated for the above coordinates. The agreement 
between experimental and theoretical values is en- 
tirely satisfactory; this confirms the validity of our 
assumed structures for the two compounds. In view 
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of the difficulty of calculating the true absorption 
factor for the single crystals, we preferred to use 
intensities from powder photographs, where this 
factor varies monotonically with scattering angle. 

For the U,Zn,7 structure we calculated theo- 
retical intensities for several lines of the 001 type 
with J ~ 3n and the hkl type with 1 = 1 and 5. Re- 
sults showed that the intensities of these lines were 
not zero, but very small. Thus, for example, for 
the 0.0.10 reflection the intensity was 0.0007, while 
for the weakest observed reflection, the 322, the 
intensity was 0.5. These results confirm the valid- 
ity of the conclusion regarding the absence of any 
systematic extinctions and the choice of space group 
Di, = Cém2. 

We also calculated the theoretical intensities 
for composition ThZng, corresponding to 6 Th and 
48 Zn atoms in the cell, three atoms being in the 
3(a) positions (at the sites corresponding to the ab- 
sence of thorium atoms in the CaZn; model). The 
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Fig. 4. Comparison between theoretical and experimental intensities on powder 
photographs of ThgZn,7. 
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Fig. 5. Comparison between theoretical and experimental intensities on powder 


photographs of 


UgZ0y7 . 
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=i) KO= 27% 


C- Zi 


O-Zny 
Fig. 6. Structure of Th)Znq7. 
agreement between theoretical and experimental 


intensities in this case was not so good as forcom- 
position Th )Zn,7; with 6 atoms in the 6(c)y] position. 


Description of the Crystal Structure 
of Thy Znji7 and U,Zni17 


The unit cells of these two compounds are 
shown both in axonometry and in projection on the 
basal plane in Figs. 2, 6, and 7. 

The structure of U,Zn,7 is a particular kind of 
noncentrosymmetric superstructure of Tho Znj7. 
Both of these structures are based on the CaZns 
motif and can be considered as superstructures of 
CaZn,;. Matters become more complex owing to the 
presence of zinc atoms in the six-fold positions 
6(c)yjj for Th)Zn,7 and in the six two-fold positions 
2(g)im,1v> 2(h)qr, m1, and 2(i)q7, ny for Uz, Zny7, since 


these positions disrupt the principal motif of CaZn,. 


We may suppose that the occupation of these posi- 
tions by zinc atoms is statistical, thus providing 
for the formation of variable-composition phases 
Thy Zny5 +17 and U,Zny5 47. If these six-fold and 
two-fold positions are not occupied by zinc atoms, 
the crystals will have compositions Th) Zn,5 and 
U,Zn;5. When the positions are fully occupied as 
in our case the compositions will be Thy Zn,7 and 
U,Zn47.- 

By analyzing Figs. 6 and 7 we come to the fol- 
lowing conclusion regarding the coordination of the 
atoms in the two compounds. In both structures the 
heavy atoms stand away from each other, being 
evenly distributed over specific positions in the 
hep Zn lattice. The nearest neighbors of each 
heavy atom will be six zine atoms set in the same 
layer at distances of 3.00 to 3.01 A. In addition, 
there are 12 zinc atoms at 3.40 to 3.41 A, lying in 
the upper and lower layers, and also (for Th, Znj,7) 
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Fig. 7. Structure of U,Zny7. 


@-2n, @-2n, 


one zinc atom at 3.12 A. In U,Zny7, in addition to 
the 18 nearest Zn neighbors mentioned, individual 
U atoms have one or two extra Zn atoms at 3.10 A 
lying on the hexagonal axes. Thus, the total coor- 
dination number for thorium atoms in Th) Zn,7 
equals 19, and for individual U atoms in U,Znj;7itis 
iO, 19, Cie ZO. 

The numbers of nearest neighbors for thorium, 
uranium, and zinc atoms in various complexes and 
the values of interatomic spacings are given in 
Table 1 for ThoZnj7 and Table 2 for U,Znj7. 

The over-all coordination number of the zine 
atoms is, on the average,12. The minimum inter- 
atomic distances Zn—Zn are 2.56 to 2.58 and 2.60 
to 2.61 A; these are 7% smaller than the sum ofthe 
atomic radii of zine (2.78 A). On the whole, the 
analysis of the Zn—Zn, Th—Zn, and U—Zn inter- 
atomic distances given in Tables 1 and 2 shows 
good agreement between the interatomic distances 
found and the sums of the atomic radii of zine (1.39 
Dy thorium (1.79 A), and uranium (1.52 aye 


Conclusions 


1. Nowotny's [1] composition and crystal struc- 
ture for the compound ThZng are erroneous. In- 
stead of ThZny we find Th)Zny7;. In the U—Zn sys- 
tem, the analogous compound has composition 
U, ZN47- 


2. The compound The Zny7 has rhombohedral 4, The structural data obtained suggest the 
structure (space group D3q) with hexagonal-lattice possibility of a variable-composition range in these 
spacings: compounds; this may be represented by the chemi- 


r cal formulas Thy ZNy5—17 and U, Zny5— 17° 
(a= 9.03 A, ¢= 13.204. 
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TABLE 1. Number of Nearest Neighbors and Interatomic Spacings (A) for Th,Zn47 


Th | Zny Zn qty Zoqir 
No. of| interat. No.of |interat No.of| interat]No. of | interat. 


neigh, dist. lneigh. | dist. |meigh dist. | meigh,| dist. 
Th in 6(c), 4 4,40 6 3.014 | 12 | 3,41 1 3,12 
Zn; iN 18 (f) | 2 3,01 3 3.04 6 | 2,66 2 3.27 
Znyzz1" 18 (h) 3 3.44 4 2.66 4 | 2.61 1 2.77 
Zny, in 9 (d) 2 3.44 4 2.66 4 | 2.64 2 2,77 
ZnyzpN 6(c),,| 4 3,12 6 3,27 CO ezeng 1 2.56 


TABLE 2. Number of Nearest Neighbors and Interatomic Spacings (A) for UgZny7 


U Zny VAN i Z0qII 
No.of |interat.| No.of | interat] No.of |interat.| No. of | interat. 
neigh. | dist, | néigh.| dist. |neigh, , dist. | neigh.| dist. 
U in | 
I(e), I (d) 2 4.39 = ~ 
I(c), I(f) — — | 6 | 3.00 12 3.40 ye 3,10 
2 (8), 2 (2) 72 (h) yz 2 (i)y 1 4.39 i 3.10 
Znyz iN 6(Z), 6 (m) 2 3.00 Salo. 00) 6 2.66 2 eu 
12 (0); + 12 (o)zy 
Zn; in 
6(n),, 6(n)y77,6(n a 
mael Pass noniid | a | a bene ee ate 
6 (n)qy, 6 (n)zy, 6 (n)yx 2 3.40 ») PRRs 
in 2 
Herp ie RM 1 3.10 Gishiane 6 Deis 1 2.58 
2 (i)rz, 2 (i) yyy | | H 
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The existence of compounds UGa and UGay in the U—Ga system is established from 
chemical and x-ray examination and their crystal structures are determined. The struc- 
ture of UGa aetermined from single-crystal data constitutes a new type. The struc- 
ture of UGa, found by the powder method is of the Al1B, type. 


Introduction 


There is no published information on the phase 
diagram of the U—Ga system. It is merely known 
that the compound UGa3 has a fcc structure of the 
AuCusz type [1]. 

Between 1949 and 1955 we studied the U—Ga 
system by x rays and discovered compounds UGa 
and UGag. 

We here present results relating to the crys- 
tal structures of these compounds. 


Preparation of the Alloys 


Uranium—gallium alloys 10-17 g in weight 
were prepared from 99.86% pure uranium and 
chemically pure gallium in a high-frequency vacu- 
um furnace. After filling with the charge, the fur- 
nace was evacuated to 1074 mm Hg. In the course 
of heating and fusion, the vacuum fell to 107! mm 
Hg. On reaching the melting point, there was a 
short rest for 2-4 min in order to achieve better 
mixing of the components. The cooling rate be- 
tween 150 and 1000° was, on the average, never 
greater than 50 deg/min. The temperature was 
regulated visually. 

The shrinkage cavities of alloys with composi- 
tion close to 50 at.% contained fine single crystals. 


TABLE 1 
Determination Determination 
of gallium of uranium 
Percentages | 
5 %U= % Gas 
% Ga |= 100—sGa ~ feo u] %U 
Weight 21.45 79.09 21.95 78.05 
Atomic 48 25 oleio 49.99 51.04 
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Chemical analysis of these gave the figures shown 
in Table 1. 

The chemical analysis was conducted by A.I. 
Elovatskii under the direction of A.G. Karabash by 
the Hillebrandt method [2], based on separating the 
gallium from the uranium by extraction of gallium 
chloride with ether from a solution of the alloy in 
6% hydrochloric acid. 


Cieystales temo re Onur 


The single crystals selected for x-ray exam- 
ination had the form of interlinked rhombic scales 
with dimensions of the order of tenths of a milli- 
meter (Fig.1). 

The x-ray study was carried out with filtered co- 
balt radiation. The symmetry of the Laue photo- 
graphs showed that the crystals had rhombic struc- 
ture. The unit-cell parameters calculated from ro- 
tation x-ray photographs (RKV-86 camera) had the 
following values: a = 9.40 A b = 7.60 AS c= 9.42 
A. 

Rotation photographs taken around the cell di- 
agonal indicated that the ab face was centered. 

In order to determine the reflection intensities, 
we took three series of 10° oscillation photographs 


03mm Odom 


Fig. 1. Form and dimensions 
of a UGa single crystal, 
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around the three principal axes of the crystal. The 
intensities were estimated with the help of a scale 
of film blackening. The cobalt radiation enabled us 
to fix only 245 independent reflections. On chang- 
ing from intensities to structure amplitudes, we 
allowed for the Lorentz, polarization, and kine- 
matic factors. 

The results of measuring and calculating the 
oscillation photographs reveal the following extinc- 
tion laws: 


hkl: h +-k = 2n 
hOl: Linares 


This corresponds to x-ray group No. 29, consisting 
of three space groups Dif = Cmmce, Ge ='Cmc2, 
On =C2cm. The smallest multiplicity of points in 
each of these groups is 4. 

The number of atoms per cell was determined 
from the additive density, preserving the smallest 
multiplicity and the least deviation from the com- 
position given by chemical analysis of the single 
crystals. The cell contains 16 uranium and 16 gal- 
lium atoms. 

The additive density of U;gGajg equals 12.8 g 
per cm® (aU = 18.97 g/cm®; da = 6.087 g/cm’). 
The x-ray density of U,,Gajg is 12.1 g/cm’. 

Analysis of the cross sections of three-dimen- 
sional F* series based on space-group symmetry 
shows that, when the structure contains atoms in 


f (hot) 


Anil 


the general position, the (u0w) section gives an un- 
equivocal indication of the space group. The ab- 
sence of atoms in general positions reduces the re- 
liability of the determination and demands the con- 
struction of additional cross sections. After ana- 
lyzing the cross sections of the F* series in the 
plane (0vw) and parallel to the coordinate axes,we 
arrived at a single centrosymmetric group Dh = 
Cmcm. 

By analyzing the cross sections of a series of 
interatomic vectors, examining the intensities, and 
making a geometrical analysis by trial and error, 
we produced a structural model giving good agree- 
ment between the experimental and computed values 
of F(0k/) and F(h0l). The computed series of 
F(hk0) values was in poor agreement with experi- 
ment. Introduction of a correction for absorption 
in the sample [83] improved agreement to a satis- 
factory extent. This was apparently because the 
crystal was least "transparent" for the (hk0) zone 
(see Fig.1). 

Graphs comparing the structure amplitudes 
Fmeas 2nd Fealc for the zero zones are shown in 
Fig.2. Our electron-density projections on the co- 
ordinate planes (001), (010), and (100) appear in 
Fig.3. The uranium-atom peaks appear clearly in 
all projections [except the 4(a) group in the projec- 
tion on the (010) plane]. The gallium-atom peaks 
are comparable in magnitude with the effects of 


Fig. 2. Fmeas and Fcalc for UGa. Thick lines = Fmeas; thin lines = 
Fealc; circles = F(hk0)meas, without allowing for absorption. 
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Fig. 3, Electron-density projection on the face of the UGa cell. 


truncating the Fourier series, so that only those 
with double "weight" appear, and even their posi- 
tions are shifted. In a check Fourier series with 
computed values of F(h0/) the gallium peaks are 
also comparable with "ghost" peaks and consider- 
ably shifted from their original positions. 

Table 2 shows the coordinates of the atoms in 
the cell. 

Table 3 gives data on coordination and the 
shortest interatomic spacings. 


Figure 4 indicates the cell in projection on the 


coordinate planes, and Fig.5 gives a view of the 
cell in axonometry. 

It is not hard to see that the UGa structure 
preserves the @ U [4] motif of broken chains along 


TABLE 2 
Posi- 
AtoM |tion | xla y|b z\c 
ie NATO | 0 0 0 
4U 4 (c) 0 0,212] 0.250 
8U 8(e) | 0,300} 0 0) 
8Ga |8(f) | 0 0,689] 0,118 
8 Ga 8(g) | 0.260] 0,354] 0.250 


Fig. 4. Positions of atoms in projections on the 
faces of the UGa cell. 


the Z axis and two bonds (for positions 4a and 8e) 
parallel to the X axis (in @U the distance is 2.85 
A, and in UGa, 2.82 A). The gallium atoms in po- 
sition 8(f) are arranged as pairs analogous to the 
Ga, "molecules" in the gallium structures (dis- 
tance in Ga, 2.45 A, and in UGa, 2.48 A) [5]. 


CRYSTALS TRUCTURE 


TABLE 3 
; No. of 
Atoms Position Distance) nearest 
neighbors 
U—U | (a)—(c)| 2.85 9 
(a)—(e)| 2,82 9 
U—Ga | (a)—(f)| 2.64 5 
(e) —(f)] 2,605 9 
(c) —(g)] 2.66 2 
(e)—(.)} 2,66 9 
Ga—Ga |(f)—(f)| 2.48 { 
(g)—(f) | 2.86 ) 


Crystal Structure of UGa, 


On indexing the x-ray photographs of UGa, by 
means of Hull curves, we found that it had hexag- 
onal structure with lattice constants: a= 4.21 A 
G= 4.01 A, c/a =0.954, 

The x-ray photographs also showed two lines 
which could not be indexed: these were apparently 
due to impurities. 

From crystal-chemical considerations it ap- 
peared the UGa, had the crystal structure of the A1B, 
type (space group Dén = C6/ mmm) with atomic po- 
sitions: 


1U_ in 1(a), 
2Ga in 2(d). 


Figure 6 gives a comparison of the theoretical 
and experimental line intensities for the powder 
photograph of UGa,. The theoretical intensities 
were calculated without allowance for absorption. 
The experimental intensities were estimated visu- 
ally on a nine-point scale. The good agreement 
between theoretical and experimental intensities in- 
dicates that the structure determined for the UGa, 
is correct. The x-ray density of UGa, equals 10.3 


O -Ga 
Fig. 5. Unit cell of UGa. 
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Fig. 6. Graphs of Imeags and Icalc for UGag. 
Thick line = Imeas; thin line = Igajc. 

g/cm®. The closest interatomic distances in the 
UGa, structure are: Ga—Ga = 2.43 A; U-Ga = 
S15 AP Ue-U 401s 

We may conclude from these data that the 
uranium atoms are not in contact with each other; 
the gallium atoms form a plane hexagonal lattice 
parallel to the basal plane, each atom having three 
nearest neighbors at a distance of 2.43 A equal to 
the smallest distance in the structure of pure gal- 
lium. 

The uranium —gallium distance equals the sum 
of the maximum interatomic radii found in the 
structures of gallium and @-uranium. 


Conclusions 


X-ray study of UGa single crystals shows that 
this compound belongs to a new structural type of 
rhombic symmetry. The lattice constants are a= 
SA0AL DS 160A c Odo 

The rhombic cell contains 32 atoms, corre- 
sponding to 16 formula units of UGa. The x-ray 
density is 12.1 g/cm®. The space group is Ds = 
Cmcm. The coordinates of the atoms are: 


4U in4 (a) 

4U in4 (c):sy= 0212 

8U in8 (e):2= 0.300 

SiGarin"Sa(7) sf 08 9 z= 118 
8Ga in8 (g):2=0,260, y=0.354. 


The UGa structure preserves the @-uranium 
motif of broken chains along the Z axis and two 
bonds parallel to the X axis. The gallium atoms in 
the UGa structure are disposed in pairs, in anal- 
ogy with the Ga, "molecules" in the structure of 
gallium. 

The compound UGay, crystallizes in the Al 
structure. The lattice constants are: a = 4.21A, 
c=4.01A,c/a=0.954. The x-ray density is 10.3 


g/cm’, 
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Chemical, x-ray structural, and phase analysis was used to establish the existence 
of the following compounds in the Th—Zn system: Th)ixZn, ThZn), |ThZn, ThyZnz, 


and Thy Zny T: 


Introduction 


No phase diagram of the Th—Zn system has so 
far been published. There are certain indications 
that the system contains compounds ThZn, [1] and 
Th, Zn [2]. Data of Nowotny [1] on the structure of 
ThZny are erroneous (see Makarov and Vinogradov 
[3]). Data on the composition and structure of the 
Th,Zn phase are confirmed in the present paper. 


Preparation of Alloys 


The original materials were 100% pure thori- 
um and 99.95% pure zinc. Solid pieces of these 
metals were coarsely filed to give powders, which 
were mixed in due proportion and pressed into tab- 
lets in a steel mold. 

In order to obtain an alloy the tablet was placed 
in a beryllium-oxide crucible provided with a cover, 
and this was placed in a quartz ampoule. After 
pumping out the air to a pressure of 10°* mm Hg 
and sealing off, the quartz ampoule was heated in 
a vertical tubular resistance furnace. The heating 
was carried out in stages: at ~400°C there was a 
rest of Ihe at 300°C. Ihe and at 1250°C,, +10 min, 
after which the alloy was cooled with the furnace. 

As a rule, the ampoule did not burst. In the 
upper part of the ampoule there was some con- 
densed zinc and the amount of this was taken into 
account. The samples weighed 8-10 g for alloys 
Nos. 3,4, and 5, and 23 g for Nos. 1 and 2. 

The original composition of the alloys, the con- 
ditions for preparing these, and the phase composi- 
tion as given by microscopic and x-ray analysis, 
appear in Table 1. 

The chemical composition of individual phases 
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The crystal structure of the first three of these was determined. 


was determined from the microstructure and chem- 
ical analysis. 

The temperature was measured with a calib- 
rated Pt—PtRh thermocouple placed at the level of 
the tablet in the immediate neighborhood of the 
outer wall of the ampoule. 

The surface of the alloys showed thin plates of 
rectangular form, having a silver-white metallic 
color similar to hoar frost. These single-crystal 
plates were studied by the Laue and rotation meth- 
ods; they had tetragonal symmetry and the follow- 
ing lattice constants: a@= 4.16 A, c= 14.40 05 

Considering the way in which these single- 
crystal plates were formed, we came to the con- 
clusion that they developed by condensation from 
the gas phase and constituted the result of interac- 
tion between the thorium and silicon reduced from 
the quartz ampoule by zinc vapor. The single- 
crystal plates were in fact the compound OThSiy. 
According to published data [2], OThSi, has a tet- 
ragonal lattice with spacings a = 4.18 A, c= 14.35 
A which are in good agreement with the above re- 
sults. 


Chemical Analysis 


A weighed portion of alloy of the order of 0.5 g 
was dissolved in aqua regia by heating. The solu- 
tion was evaporated almost to dryness, 3-5 mlcon- 
centrated nitric acid was added, and the result 
again evaporated almost to dryness. The thorium 
and zinc nitrates obtained in the dry residue were 
dissolved in a small quantity of water with the ad- 
dition of 1 ml of 5% nitric acid solution. 

The resultant solution was brought up to avol- 
ume of 40 ml, 2 g ammonium nitrate was added, 
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TABLE 1 
Composition in 
No. original mixture Phase composition of alloy Remarks 
of alloys} wt.% Zn| at.% Zn 
1 21.98 50.0 Alloy consists of two parts; (a) | Very small a-ThSip crystals 
uniform silver-white,easily oxi- | formed on surface 
dized phase Thy; x with small 
inclusions of an unknown phase; 
(b) loose dark-gray mass of 
thorium dendrites and cubic 
phase X 
2 36,03 66.66 | Single-phase alloy, ThyZn, Zinc (1.2 g) sublimes. Crys- 
tals of a-ThSi, formed on 
alloy 
3 46.0 75.0 Alloy consists mainly of phase 
ThZn, and traces of ThZ ny 
4 52.99 80.0 Alloy consists mainly of phase | Crystals of wa-ThSi, formed 
ThZn, on surface 
5 Talal 90.0 Alloy consists of phase ThyZnyz | Large a-ThSi, crystals formed 
on surface 


and the result heated to 60-80°C. From the hot 
solution thus obtained, thorium was precipitated as 
the peroxide by adding a 15% solution of H,O,. The 
filtered and washed residue was roasted to constant 
weight in a porcelain crucible at about 1000°C and 
weighed as ThO,. 

In order to determine the zinc, the remaining 
filtrate was evaporated almost to dryness. The 
dry residue was dissolved in water with 1-2 mlhy- 
drochloric acid added and diluted with water to a 
volume of 100 ml, after which 1 g ammonium acet- 
ate and 5 ml concentrated acetic acid were added. 
The zinc was precipitated as ZnS by H,S. The re- 
sultant ZnS residue was washed and roasted so as 
to transform it to ZnO, from which the original 
zinc content was determined by weighing. 


Crystal Structure of Tho,,Zn 


From alloy No. 1 we were able to isolate a 
few single-crystal fragments around 0.1 mm in 
size and to study these by the Laue and rotation 
methods. The Laue x-ray photographs indicated 
tetragonal symmetry and the rotation (around the 
V 2a and c directions) photographs gave periods of 
10.68 and 5.62 AG respectively. Hence, the dimen- 
sions of the tetragonal cell are: a = 7.62 A, c= 
5.62 A, 

According to published data [2], the cell di- 
mensions of Th.Zn are: a= 7.95 A, c= 5.64 A. 

The slight difference in the a values may be 
due to the variable composition of the Tho+xZn 
phase. This is partly confirmed by the fact that, 
according to chemical analysis, the unbroken part 


of alloy No. 1, in which the Th)i,Zn single crys- 
tals were found had composition 8.19 to 8.94wt.% 
Zn, corresponding to 24.1 to 25.8 at.% Zn, in- 
stead of the theoretically expected value of 33.33 
at.% Zn for the composition Th,Zn. 

From crystal-chemical considerations and the 
celldimensions, it appeared that Tho; ,Zn had the 
CuAl,-type structure [4]. In this case, the theo- 
retically calculated density equals 10.7 g/cm’, in 
satisfactory agreement with the value 10.2 g/cm? 
for the additive density of Thy;xZn. It was very 
hard to determine the experimental density owing 
to the small weight of the single crystals, and it 
was in fact not determined. 

In the case of the CuAl, structure (space group 
Dies the positions of the atoms in Tho, xZn for 
composition Th, Zn will be: 


4Zn in 4(a), 8Thin 8(h). 
For the parameter u we took roughly '%. 
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Fig. 1. Comparison between calculated and experi- 
mental spot intensities for the zero-layer line of 
ThpZn., 
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Fig. 2. Comparison between calculated and experimental intensities for ThZn, 
spots: a) zero layer line; b) first layer line. 


Figure 1 gives a comparison between calcu- 
lated intensities for the ideal composition Th, Zn 
and the experimental values for spots on the zero 
layer line of rotation photographs taken around the 
c axis of a Tho; ,Zn single crystal. 

Remembering that in calculating the theoreti- 
cal intensities the absorption factor was not taken 
into account, we may consider the agreement be- 
tween theory and experiment as quite satisfactory. 
This leads to the conclusion that Thi ,Zn does in 
fact have the CuAl, structure. 

The interatomic distances calculated from the 
cell dimensions obtained together with the parameter 
u= Ms generally agree with the corresponding sums 
of atomic radii of thorium and zinc; in order to ob- 
tain the exact interatomic spacings, however, the 
exact value of u is required, and this we did not 
have. The interatomic distances in Tho;,Zn have 
the following values: 


Va 27n = 2.81 A. Zn —8Th= 348A. 


Th 470 = 3.18 A, 
Th—2Th = 3.33A, Th—1Th —3.59A. 


Crystal Structure of ThZn, 


As seen from Table 1, alloy No. 3 contained 
the phase ThZn,. We were able to select a few 
very fine acicular single crystals with clearly ex- 
pressed hexagonal-prism faces. The thickness of 
the crystallites was less than 0.1 mm and the length 
about 0.3 mm. In air the ThZn, crystals were 
quite stable and lasted for several days. 

Rotation and oscillation x-ray photographs gave 
coarse values for the identity periods of the hex- 
agonal ThZn, cell; these were used to index the 
powder photograph, which gave the following lat- 
tice spacings: a = 4.20 A; c= 4.17 A. 


By comparing the atomic radii, cell dimen- 
sions, and additive density, we considered that 
ThZn, had the AlB,-type structure with three atoms 
per unit cell. 

A comparison of the theoretically calculated 
and experimental spot intensities of the zero and 
first layer lines of the rotation x-ray photographs 
of ThZn, is given in Figs. 2a and 2b. 

In calculating the theoretical intensities we 
considered the Lorentz and recurrence factors.The 
agreement between the experimental and theoreti- 
cal intensities is entirely satisfactory, indicating 
the correctness of the hexagonal AlB,-type struc- 
ture for ThZny. 

The closest interatomic distances in ThZn, 
have the following values: 


Zn — 3Zn = 2,.42A, Th—12Zn = 3.18 A, 
Zn—6Th=3.18A, Th—2Th=4.17A. 


If we assume that the atomic radii of thorium 
and zinc are, respectively, 1.80 and 1.38 A, then 
the Zn—3Zn distance in the hexagonal layers is 
12.3% below the atomic diameter. The Th—Zn dis- 
tance found agrees exactly with the sum of the atom- 
ic radii of thorium and zinc. 


Crystal Structure of Th2Zn4 


According to microscope analysis, alloy No.3 
(see Table 1) had an almost single-phase micro- 
structure. A slight accumulation of acicular ThZn, 
crystals occurred owing to segregation in the lower 
part of the ingot obtained from alloy No.3. 

Chemical analysis of monolithic fragments of 
alloy No.3 gave a zine content of 52.22 wt.%, cor- 
responding to 79.5 at.% Zn, close to ThZn,. 
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Fig. 3. Comparison between theoretical and experimental intensities for ThZn, 
spots: a) zero layer line; b) first layer line. 


Single crystals of ThZn, could easily be iso- 
lated in alloy No.3; these had the form of lamellas 
with curved edges. From x-ray rotation photo- 
graphs we obtained the following parameters for the 
tetragonal cell of ThZn,: a= 4.26 1 c=10.4A. 

From the x-ray densities of thorium, zinc, and 
the intermediate phases Thoi;Zn and ThZn,, we 
constructed a curve for the x-ray densities of the 
Th—Zn alloy system. The value interpolated on 
this curve for ThZn, was 8.6 g/cm®, giving ten 
atoms per cell, or two ThZn, groups. 

From crystal-chemical considerations we ga- 
thered that ThZn, had the AlB, type of structure. 
This was confirmed on comparing theoretical and 
experimental intensities for the spots on the zero 
and first layer lines of the rotation photographs, as 
seen in Figs. 3a and 3b. 

The space group is Dit —14/mmm. The co- 
ordinates of the atoms are: 2Th in 2(a); 4Znj in 
4(d); 4Znjj in 4(e). We took the z parameter as 
0.38, 

The closest interatomic spacings in ThZn, have 
the values: 

Th — 8Zn,— 3.26A, Th — 8Zny = 3.36A, 
Zap —AZny = 2.48A,, 2ny — 470 = 2.024, 
Zny—4Th = 3.26A, Znoy—4Zn; = 2,52A, 


Znyy—4Th — 3.36A. 
Comparing these with the sums of the atomic 
radii, we see that the Zn—Zn distances are severe- 


ly shortened. This indicates the presence of a homo- 
polar or covalent bond. 


The Th,Zn, Phase 


As seen from Table 1, alloy No.2 had original 
composition 66.6 at.% Zn. In the course of fusion, 
part of the original zinc volatilized and condensed 
in the upper part of the ampoule. The condensed 


zinc was weighed and the composition of the alloy 
corrected for the actual zinc content, equal to 32.6 
wt.% (63.6 at.%). Chemical analysis of alloy No.2 
gave 34.38 wt.% Zn, or 65.2 at.%. Alloy No. 2, 
as shown by study in the microscope, was single- 
phase. Thus, the chemical composition indicates 
the formula ThyZn;, or possibly ThgZnj5. 

We were unable to find single crystals corre- 
sponding to this phase, since the alloy constitutesa 
continuous monolithic formation. Hence, powder 
photographs were obtained from alloy No.2, using 
copper radiation. The Th,Zn; phase was very brit- 
tle and susceptible to work-hardening, as indicated 
by the reduced line intensities. In Table 2 we there- 
fore present a calculation of two x-ray photographs 
somewhat differing from one another in the number 
and intensity of the lines observed and correspond- 
ing to a fragment of the alloy and Th,Zn7 powder, 
respectively, the latter obtained by grinding in an 
agate mortar. We could not index the x-ray photo- 
graphs, and the question of the crystal structure of 
Th,Zn; remains open. 


Phase X 


As seen from Table 1, the upper part of alloy 
No. 1 contains a loose dark-gray metallic mass of 
thorium dendrite and a cubic phase X. Since we 
were unable to obtain phase X in pure form so as 
to analyze it chemically, the chemical composition 
of this phase remains uncertain. 

As regards the ratios of the line intensities, 
the x-ray photograph of phase X is very similar to 
that of a Th, and also that of ThO,. The lattice 
constant of thorium dioxide (a = 5.58 A) is close to 
that of the phase in question (a= 5.684 iN which 
suggests that phase X may have an oxide nature. 
We can imagine that the increased lattice constant 
of phase X (by comparison with thorium dioxide) is 
due to the dissolution of zinc oxide in thorium di- 
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TABLE 2 
Powder photograph of X-ray photograph Of 
ThyZny ThyZn, fragment 
No. of No. of : 
lines 1 sin? lines 1 sind 
1 w-m 0.2014 1 S. 0.2862 
o w-m 0.2636 2 m. 0.3423 
3 S. 0.2904 3 w-m 0.4470 
4 Ww. 0.3121 4 w-m 0.4637 
5 m-s 0.3447 D m-s 0.5400 
6 m. 0.4057 6 m-s 0.5659 
tt m. 0.4184 i m*+s 0.7428—0.7250 
8 w-m 0.4480 8 m-s 0.7424 
9 W. 0.4651 Yq m. 0.7923—0.7997 
40 v.W 0.4909 10, m-s 0.8119—0.8180 
14 v.W 0.5038 14 m-s 0.8345 
12 Ww. 0.5356 dy m-s 0, 8494 
13 W 0.5400 13 m-s 0.8601 
14 w-m 0.5659 14 m-s 0.8839 
15 w-m 0.5955 1936 S. 0.9237—0 .9269 
16 Ww. 0.6298 16 m. 0.9350 
i Ww. 0.6576 17, Ss. 0 9446—0.9484 
18 Ww. 0.6693 18 m, 0.9613 
19, Ww. 0.7143—0.7288 19. M-s 0.9640—0. 9667 
20, w-m 0.7383—0.7441 20,, S. 0.9762 
21 w-m 0.8165 Za m. 0.9795 
22 w-m 0.8369 le Ss. 0.9849—0 . 9864 
23 w-m 0.8625 
24 m. 0.8840 
Ze m. 0.9484 


oxide. The true nature of phase X awaits further 
study. 


The Th ,Zn,7 Phase 


As already indicated in the introduction, 
Nowotny's data [1] on the chemical composition and 
crystal structure of ThZn, are erroneous. Makarov 
and Vinogradov showed [3] that instead of Nowotny's 
phase of composition ThZn,, the actual phase 
formed was Th,Znj;, having a structure with the 
following hexagonal lattice spacings: a = 9.03 A, 
e=13.20 A, 

Details of the Th, Zn,7 structure are given in 


[3]. 


Conclusions 


1. X-ray phase analysis of alloys in the Th— 
Zn system, taken together with chemical analysis 
and a study of microstructure, have established 
the existence of the following chemical compounds: 


Tho+Zn, ThZn,, ThZn,, Th,Zn,, Th.Zn,, 


and a phase X of unknown chemical composition. 
The crystal structures of some of these phases 
have been determined by x-ray diffraction study of 
powders and single crystals. 


2. The phase Tho: Zn has variable composi- 
tion shifted to the side of greater thorium content 
and a Structure of the CuAl, type with spacings a = 
7.62A, c=5.62A. 

3. The phase ThZn, has the structure of AlBy 
with spacings a = 4. 20A, c=4,17A. 

4. The phase ThZny has the structure of AlBg 
with spacings a = 4.26 A, c=10.4A. 

5. The phase ThyZn,7 has rhombohedral struc- 
ture with spacings a = 9. o3A, @.= 13:20 A, 

6. Phase X has an unknown chemical compo- 
sition, but its fcc lattice has constant a = 5.68,A. 
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INOwE Gn 


A wide range of Pb(Ti, Sn)O, solid solutions extending from 0 to 75 mol.% "PbSnO3" 
was found to exist in the PbTiO,—"PbSnO3" system. The phase diagram of this solid 
solution is similar to that of Pb(Ti, Zr)O3 but differs from that of Ba(Ti, Sn)O3. Con- 


clusions are drawn regarding differences 


in the mechanism of spontaneous electric 


polarization in the ferroelectrics BaTiO; and PbTiO3, which have hitherto been re- 


garded as completely analogous. 


Introduction 


We showed earlier [1,2] that, when the oxides 
PbO and SnO, were sintered in the temperature 
range 800 to 1050°C, no single compound PbSnO3 
resulted; samples of composition PbSnO3 contained 
two phases, Pb).SnO, and SnO,. The samples exam- 
ined in the present paper may thus be regarded as 
belonging to a linear section of the ternary system 
PbTiO; — PbySnO, —SnO,, which can only nominally 
be written in the form of a binary system PbTiO;— 
"PbSnO3." Despite the two-phase nature of its com- 
ponents, however, this system PbTiO3—"PbSnO," 
may be expected to contain a wide range of solid 
solutions, in view of the similarity between the ra- 
dii of the tetravalent ions of titanium (0.64 A) and 
tin (0.67 A), which belong to the same group of the 
Mendeleev periodic table. The likelihood of this 
being so is supported by the fact that there are con- 
tinuous series of solid solutions in systems similar 
to the one considered, namely, BaTiO3—BaSnOg3 [8, 
4] and PbTiO,— PbZrOs [5,6]. 

A study of the PbTiO3—"PbSnO3" and BaTiO; — 
BaSnO3 systems is necessary for a deeper under- 
standing of the properties of the ferroelectrics 
PbTiO; and BaTiOs. Those who have studied the 
ferroelectric properties of PbTiO; [7-11] and the 
atomic displacements in these [12], consider the 
two compounds to be completely analogous. It is 
supposed that in the cells of both compounds the 


displacement of the titanium ion (below the Curie 
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temperature) from the center of the oxygen octa- 
hedron along the four-fold axis leads to the develop- 
ment of electric polarization in the crystal in the 
same direction. If the properties of these two com- 
pounds are similar, then the phase diagrams of the 
solid solutions Ba(Ti,Sn)O3 and Pb(Ti, Sn)O3 should 
also be similar. A study of the solid solution 

Pb(Ti, Sn)Os3 is also interesting for elucidating the 
effects of the different sizes and electron configura- 
tions of the mutually interchangeable cations of ti- 
tanium and tin on the ferroelectric properties. 


Characteristics of the Samples 


Samples of the PbTiO;— "PbSnO3" system were 
prepared in the Physics Laboratory of the Institute 
of Silicate Chemistry, Academy of Sciences of the 
USSR, by the well-known technique for preparing 
ceramic ferroelectrics containing lead [9]. The 
original materials included titanium dioxide and tin 
("pure" type) and commercial red lead. The tem- 
perature of the final anneal fell from 1060 to 1000 
°C with increasing "PbSnO3" content. The dielec- 
tric properties of these samples were studied 
simultaneously in the same laboratory. The results 
of these dielectric studies are presented in [13]. 


Method 


X-ray photographs were taken at room tem- 
perature in the RKU-114 camera. For high- and 
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low-temperature work special cameras of diam- 
eters 120 and 114 mm, respectively, were used. 
Both the low-temperature camera and the special 
apparatus for blasting and glazing the samples were 
developed in the L. Ya. Karpov Physicochemical 
Institute [14,15] and enabled photographs to be taken 
from room temperature down to the temperature of 
liquid nitrogen. In selecting radiation for the pre- 
cision measurement of parameters, a specialnom- 
ogram was used [16]. All the samples were taken 
in Cu Ka radiation, and the sample 40% PbTiO;— 
60% "PbSnO;" in NiK and Cr K radiation as well. 
For work with Cu Ka radiation the 314 and 413 lines 
of lead titanate lay at angles @ of 76 and 80.5°, re- 
spectively. 

The films were placed in the RKU-114 camera 
in an asymmetric fashion. The positions of the 
lines used for calculating the cell parameters were 
measured to an accuracy of A(2L)=+ 0.1 mm, and 
the value of 7Deff determined to an accuracy of 
+0.01-0.02 mm. Under these conditions the accu- 
racy of determining the cell parameters should have 
been 0.001 A. However, owing to the large therm- 
al-expansion coefficient along the c axis of the lead- 
titanate cell (about 5 - 107>/ deg) [11], the error in 
determining the cell parameters reached +0.0015 
A for a temperature variation of 3 to 5°. Thus, the 
errors in determining the axial ratio c/a and vol- 
ume v of the cell will,respectively,be A (c/a) = 
+0.0007 and Av = +0.07 A®. For the rhombohedral 
cells the angle @y, was determined both by the 
usual method (solution of two equations with twoun- 
knowns) and by the method described by Megaw [17]. 
The values of @yp calculated for several groups of 
lines: 422 (Cu Ka@ radiation), 332 (NiK radiation), 
and 222 (CrK radiation) differed by +1'. This ac- 
curacy in determining c/a, v, and Oy gave a re- 
liable check on changes taking place in the cell of 
the Pb(Ti, Sn)O3 solid solution with varying chemi- 
cal composition. The wavelength values for the cal- 
culation were taken from [18]. 


Results of Examining Samples at 
Room Temperature 


Our x-ray diffraction photographs showed that 
the PbTiO;—"PbSnO3" system contained a wide 
range of solid solutions on the PbTiO; side, extend- 
ing from 0 to 75 mol.% "PbSnOQ3." The range from 
75 to 100% "lead stannate" was not single-phase. In 
this region, x-ray examination revealed the phases 
Pb(Ti, Sn)O3, X, Pb,SnOQ,, and (Sn, Ti)O,. Condi- 
tions for obtaining phase X in the PbO—SnO, sys- 
tem were indicated in [1,2]. The probable bounda- 
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Fig. 1. Phase composition of the samples 
in the PbTiO,;—"PbSnO," system studied. 


aries between these phases in the range 75 to 100 
mol.% "PbSnO3" appear in Fig.1. 

For pure lead titanate we confirmed a tetrago- 
nal cell with parameters close to those given in [17, 
19]. With increasing "lead stannate" content, the 
tetragonal distortion of the cell diminished and for 
a 55 mol.% content of "PbSnO 3" the tetragonal cell 
transformed into rhombohedral. The rhombohed- 
ral distortion fell with increasing "lead stannate" 
content but still existed up to the boundary of the 
solid solution Pb(Ti,Sn)O3. The phase transforma- 
tion at 55 mol.% "PbSnO3" was confirmed by the 
following. On the x-ray photographs of samples con- 
taining 0-50 mol.% "PbSnO3," as the concentration 
of the latter increased the components of all groups 
of lines gradually moved closer together; between 
50 and 60 mol.% there was a considerably greater 
change in the diffraction pictures than that corre- 
sponding to a 10% change in the lower range (0 to 
50%). Careful analysis of the way in which various 
groups of lines (especially the 332 group) behaved 
in the x-ray photographs of these samples showed 
that between 50 and 60 mol.% "PbSnO3" the cell of 
the Pb(Ti, Sn)O3 solid solution passed from tetrago- 
nal to rhombohedral. Figure 2 gives a schematic 
representation of the 332 group of lines found on 
the 60-mol.% photograph for Cu and Ni radiations 
and the 222 group for Cr radiation. In the same 
figure (for comparison) we have the same groups of 
lines for the tetragonal and rhombohedral distor- 
tions of the lattice, respectively. We see from 
Fig.2 that the observed 332 and 222 groups of lines 
on the x-ray diffraction picture of Pb(Tizo, Sngq)O3 
are indexed on the assumption of a rhombohedral 
distortion of the cell and not on that of a tetragonal 
distortion. 

Figure 3 shows the results of a precision meas- 
urement of the cell parameters of solid solution 
Pb(Ti, Sn)O3. According to [20,21] the transforma- 
tion of one form into another takes place over a 
certain range, and the nature of the change in pa- 
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Fig. 2. Groups of 332 and 222 lines for: a) rhombo- 
hedral and c) tetragonal distortion of the cell, and 
also b) for the experimental x-ray diffraction picture 
of the sample Pb(TiggSngp)Oz. 


rameters must be similar to that shown in Fig. 3. 
In the transitional region in which the tetragonal 
and rhombohedral forms coexist simultaneously, 
the parameters change sharply on passing from one 
form to another. The axial ratio c/a for the tet- 
ragonal cell at the beginning of the transitional 
range is,moreover,larger than at the end. 

Figure 4 shows the variation of the axial ratio 
c/a and angle @pp of the solid-solution cell on the 
"PbSnO3" content. Figure 5 shows the change in 
unit-cell volume for the solid solution as a function 
of composition. The unit-cell volume of the solid 
solution rises with increasing content of "PbSnO3". 
This is because the tetravalent tin ion has a larger 
radius (0.67 A) than the tetravalent titanium ion 
(0.64 A). In the range 0 to 20 mol% "PbSnOs," 
the cell volume rises nonlinearly, and in the range 
20 to 75 mol.%, linearly. On transforming from 
the tetragonal to the rhombohedral form, the unit- 
cell volume changes smoothly (within the limits of 
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Fig. 3. Unit-cell parameters qa and c of solid 
solution Pb(Ti,Sn)Og as a function of composi- 
tion, 


measuring error), although we might have expected 
a sharp change, since ferroelectric transforma- 
tions are transformations of the first kind [22]. 


Phase Diagram of the Solid Solution 
Pb(Ti, Sn)O3 


In order to construct the phase diagram of the 
solid solution Pb(Ti, Sn)O3 experiments at room 
temperature are insufficient. We therefore made 
some high- and low-temperature measurements of 
Pb(Tigy,Sn49)03 and low-temperature measurements 
of PbTiOs. 

The x-ray diffraction photographs of PbTiO; 
and the sample containing 40 mol.% "PbSnO3" taken 
at liquid-oxygen temperature (-183°C) showed an 
absence of phase transformations in these samples 
between —183°C and room temperature. 

The cell parameters of lead titanate at —183°C 
and room temperature are given in Table 1. We 
see from this table that the axial ratio at —183°C is 
slightly larger than at room temperature. The unit- 
cell volume is practically constant over this range. 
At the same time it is well known that between the 
Curie point (490°C) and room temperature the unit- 
cell volume of lead titanate rises [10,11]. The ob- 
served fact that the volume remains constant may 
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Fig. 4. Variation of c/a and a;p of the solid 
solution Pb(Ti,Sn)O, with composition, 
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Fig. 5. Variation of the unit-cell volume of 
solid solution Pb(Ti,Sn)O3 with composition. 
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TABLE 1. Cell Parameters of Lead Titanate 
Temp, afk | cA | cla | v 
Control 9, 900\)| 4.458 | 406409" 63.45 
—183°C Broo) Noe or We OF. |" 638.0, 


be explained if the increase in unit-cell volume due 
to the temperature variation of polarization is com- 
pensated by ordinary thermal compression. 

The absence of phase transformations in PbTiO, 
on cooling with liquid oxygen agrees with dielectric 
measurements given in [10,11]. The sample con- 
taining 40 mol.% "PbSnO3" was also photographed 
at high temperature. The main aim of the high- 
temperature work was to discover the Curie point. 
For the sample in question the Curie temperature 
as found by x-ray diffraction was 320 + 5°C. This 
shows that the Curie point of the solid solution 
Pb(Ti, Sn)O3 falls as the "PbSnO3" content rises, 
since, for pure PbTiO; the Curie temperature is 
490°C [9,11, 23]. 

The transformation of the tetragonal form of 
the solid solution into rhombohedral and the reduc- 
tion in the Curie point as the "PbSnO3" content 
rises indicates a similarity between variations in 
the properties of the solid solution in question, 
Pb(Ti, Sn)O3, and the well-known solid solution 
Pb(Ti, Zr)O3 [5,6,24-26], the phase diagram of 
which is shown in Fig. 6. On replacing the titani- 
um in PbTiO; by either tin or zirconium, the Curie 
temperature falls. Moreover, for a certain con- 
tent of either tin or zirconium there is a transfor- 
mation from the tetragonal modification to the 
rhombohedral. 

Let us take the variation in the Curie tempera- 
ture in the PbTiO;—"PbSnO3" system as linear and, 
in analogy with the PbTiO;— PbZrO; system, the 
boundary between the tetragonal and rhombohedral 
forms on the phase diagram as vertical.In this case 
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Fig. 6. Phase diagram of Pb( Ti,Sn)O, 
from data of [6]. 
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lutions Ba(Ti,Sn)O, from data of [27]. 


the diagram of the phase state of the solid solution 
Pb(Ti, Sn)O3 will appear as indicated in Fig. 7. 

As established in [22], the similarity between 
the ways in which zirconium and tin ions act on the 
ferroelectric properties of PbTiO 3 on introducing 
them into its lattice instead of titanium ions is due 
to the fact that the dimensions of the Zr** ion (0.82 
A) and Sn** ion (0.67 A) exceed those of the Ti** ion 
(0.64 A). 

The phase diagram of the solid solution 
Pb(Ti, Sn)O3 differs from the well-known phase dia- 
gram of Ba(Ti, Sn)O3 [27] shown in Fig. 8. 

While the present investigation was being car- 
ried out, two other papers [28,29] were published 
on the dielectric and structural properties of the 
PbTiO; —"PbSnO3" system. Our conclusions regard- 
ing the existence of a solid solution in the range 0 
to 75 mol.% "PbSnO3" agree with the results of 
these papers. At the same time the two-phase na- 
ture of the samples with composition PbSnO3 and 
the presence of phase X in the heterogeneous re- 
gion was not noted by the other authors, who, more- 
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over, erroneously assumed that the unit cell of the 
solid solution remained tetragonal over the whole 
solid-solution range, and that on increasing the 
"PbSnO3;" content from 50 to 60 mol.% the tetragon- 
al distortion of the cell fell sharply. Figures 9 and 
10, taken from [29], show the corresponding varia- 
tions in the cell parameters of the solid solution 
and the axial ratio c/a as functions of the "PbSnO3" 
content. 

The erroneous data on the structure of the solid 
solution Pb(Ti, Sn)O3 in [28,29] make it possible to 
understand a certain lack of coordination in these 
papers. Thus, existing experimental data on anum- 
ber of solid solutions show that the behavior of the 
Curie temperature in a system may be judged from 
that of the axial ratio. Assuming the validity of the 
x-ray data given in [28,29], we should then have ex- 
pected that in the range 50 to 60 mol.% "PbSnO3" 
the Curie temperature would have fallen sharply, 
whereas, in fact, this did not occur, the papers in 
question indicating a practically linear variation of 
the Curie point over the whole range of the solid 
solution and this result also being confirmed by di- 
electric measurements by other authors [13]. 

Since the Japanese authors of [28,29] lacked 
correct data on the phase diagram of the solid so- 
lution Pb(Ti, Sn)O3, they failed to note the simila- 
rity between the phase diagrams of Pb(Ti, Sn)O3 and 
Pb(Ti, Zr)O3 and the difference between these and 
that of Ba(Ti,Sn)O3. This is perhaps why these 
authors mistakenly assigned a leading role to the 
difference in the electron structures of the titani- 
um and substituted tin ions. As we have already 
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pointed out, our own results lead to a different con- 
clusion, namely,that the properties of the solid so- 
lution Pb(Ti,Sn)O3 depend substantially on the dif- 
ferent sizes of the substituent cations; the differ- 
ence between the electronic structures of titanium 
and tin is of subsidiary significance. 


Discussion of Results 


Our most important experimental result is the 
establishment of a difference between the phase 
diagrams of solid solutions Pb(Ti, Sn)O3 and 
Ba(Ti,Sn)O3, which indicates a lack of agreement 
between the properties of PbTiO; and BaTiO3. In 
this connection it is interesting to compare the prop- 
erties of these two solid solutions. From the re- 
sults of [6,27,30], which were published while the 
present work was in progress, we may conclude 
that the phase diagram of the solid solution 
Ba(Ti, Zr)O3 is similar to that of Ba(Ti, Sn)Os, but 
different from that of Pb(Ti, Zr)O3. Thus, the dif- 
ference between the phase diagrams of solid solu- 
tions Pb(Ti, Zr)O3 and Ba(Ti, Zr)O3 also indicates 
that PbTiOs and BaTiO; are not completely analogous. 

Let us consider the factors which may lead to 
a difference in the properties of these two ferro- 
electrics. 

According to existing estimates of the geo- 
metrical criterion t for different compounds of the 
perovskite type [7,17,31], the value of t for BaTiO, 
is greater than unity. In other words, geometrical 
analysis indicates that the titanium (B-type) cation 
is "free" in the BaTiO; cell. As regards the t fac- 
tor of lead titanate, a number of workers give a 
value slightly greater than unity [7,31], while 
Megaw gives one slightly smaller than unity [17]. 
In the case of antiferroelectric lead zirconate, 
everyone [6,17,31] agrees that t is considerably 
smaller than unity, i.e., the lead (A-type) cation 
is "free" in the cell of this compound. It is ap- 
parently for this reason that it is assumed in [32] 
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that ferroelectrics may include ABO; compounds of 
the perovskite type, in the cells of which the B ca- 
tions are free (t > 1), while antiferroelectrics in- 
clude only those compounds where this applies to 
the A cations (t < 1). 

There are some experimental data, however, 
which disagree with this last conclusion. Thus, in 
[6,33] it was shown by x-ray diffraction that, in the 
cells of both the ferroelectric rhombohedral and 
antiferroelectric pseudo-monoclinic form of the 
solid solution Pb(Ti, Zr)O3, cations of the same type 
(lead) were free, so that at room temperatures the 
lead cations in the cells of these solid solutions 
were considerably displaced relative to the com- 
pact BO, = (Zr, Ti)Og octahedra along three-fold and 
two-fold axes of the original cells, respectively. 
Thus not only compounds with free B cations, but 
also those with free A cations, may have ferroelec- 
tric properties. 

The difference between the properties of the 
ferroelectrics BaTiO3 and PbTiO; isprimarily con- 
nected with the fact that in the case of PbTiO; the 
lead (A) cation is free in the cell rather than the ti- 
tanium (B) cation as in BaTiOs. 

On replacing the titanium ions in PbTiOs by tin 
or zirconium ions, the geometrical possibilities 
for the displacement of the lead cation increase.For 
the solid solution Pb(Ti, Zr)O3 it may be shown that 
cations of the A type (Pb in the present case) may 
be displaced along the four-fold (tetragonal form), 
three-fold (rhombohedral form), and two-fold 
(pseudo-monoclinic form) axes, in the same way as 
B-type cations (for example, Ti in BaTiOs). 

The displacement of the titanium cation in 
BaTiO; (below the Curie point) from the center of 
the oxygen octahedron leads to a polar rearrange- 
ment of the cell. In this sense the titanium cation 
in BaTiO; is a ferroactive cation. From this point 
of view the ferroactive cation in PbTiO; is the lead. 

Workers studying atomic displacements in 
PbTiO; by x-ray and neutron diffraction [12] esta- 
blished a picture of atomic displacements in PbTiO; 
differing somewhat from that in BaTiO; and ex- 
plained this as being due simply to the different po- 
larizabilities of the barium and lead cations; they 
had no doubts regarding the similarity between the 
properties of BaTiO; and PbTiO3, but found con- 
firmation for this in the continuous variation of the 
Curie temperature and cell parameters of solid so- 
lution (Ba, Pb) TiO; with composition, as establishec 
in [34,35]. No attention was paid to the anomalous 
temperature dependence of the cell parameters of 
solid solution (Ba, Pb)TiO3 in the concentration 
range 10-25 mol.% PbTiO;, established in [36]. 
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Up to the present time, the different mecha- 
nisms of electric polarization in PbTiO; and BaTiO; 
have not been considered in any of the known pub- 
lications [9, 37-39] dealing with factors responsible 
for the different Curie temperatures of these com- 
pounds. Hence, the conclusions drawn in these 
papers regarding factors affecting the Curie tem- 
peratures of the BaTiO; group of ferroelectrics re- 
quire checking. The Curie temperature of PbTiO; 
is determined by the energy state of the lead ca-« 
tion, and that of BaTiO; by the energy state of the 
titanium cation. 

We showed in [22] that, for ferro- and antifer- 
roelectrics with t < 1 and cations A of the same 
valence, the Curie temperature was the higher, the 
greater the polarizability of the ferroactive A ca- 
tion and the smaller the cell parameters. The fact 
that PbTiO; has a comparatively high Curie point 
(490°C) and SrTiO; is not a ferroelectric is ex- 
plained by the different polarizabilities of the lead 
and strontium cations. Both the lead and strontium 
cations have geometrical capabilities of displace- 
ment, but in fact only the lead cation is ferroactive 
in view of its high polarizability compared with that 
of strontium. 

The results obtained in the present investiga- 
tion enable us to explain the difference in the prop- 
erties of solid solutions (Ba, Sr)TiO3 and 
(Ba, Pb) TiO3. 


Conclusions 


1. A range of solid solutions extending from 
0 to 75 mol.% "PbSnO3" has been established in the 
Pb TiO; — "PbSnO3" system. 

2. Fora "PbSnO3" content up to 55 mol.% the 
cell remains tetragonal; in the range 55 to 75mol.% 
it becomes rhombohedral. 

3. The phase diagram of solid solution 
Pb(Ti, Sn)O3 has been obtained for various tempera- 
tures from x-ray data; it is similar to that of 
Pb(Ti, Zr)O3, but differs from that of Ba(Ti, Sn)O3. 

4. The mechanisms of electric polarization in 
ferroelectrics BaTiOg and PbTiOs3, hitherto con- 
sidered completely analogous, are apparently dif- 
ferent. 

5. The different properties of BaTiO; and 
PbTiO; are mainly explained by the fact that in 
BaTiO; the ferroactive cation is titanium (B type), 
while in PbTiO; it is lead (A type). 

We consider it our duty to thank Doctor of 
Physicomathematical Sciences G.A. Smolenskiifor 
providing the samples for study and discussing 
some of the results. 
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Note Added in Proof 


Some of the experimental results in our earlier 


paper [1] are independently confirmed by B. Jaffe, 
R.S. Roth, and S.Marzullo (J. Res. Nat. Bur. Stds., 
55, 239-254, 1955); we only became aware of this 
after going to press. 
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The compound SiB, crystallizes with cubic symmetry andis isostructuralwithCaB,. The 
lattice parameters are a =4.142 + 0.002 kX, 7 pyen = 2.15, Ox ray = 2.18 g/cm; 


Z=1, space group O}. 


Moissan [1] was the first to report a compound 
of silicon with boron. He reported that the fusion of 
silicon with boron gave two compounds, SiB3 and 
SiBg, which were very hard and conducted electri- 
city. In a handbook [2] it is stated that SiB3 consists 
of black rhombic crystals with a density of 2.52 g 
per em®, while black crystals of SiB, have a den- 
sity of 2.47 g/cm’. 

The samples of SiBg which we have studied 
were prepared by G.V. Samsonov (Moscow Insti- 
tute of Nonferrous Metals and Gold) from element- 
al silicon and boron and also by reducing a mixture 
of boric oxide and quartz with magnesium. The 
synthesis from the elements was carried out by sin- 
tering a mixture of powdered boron and silicon ina 
hot press. The samples were black powders. 

The x-ray photographs were taken in a 114.4 
mm RKU camera with CuKa radiation. Initial cal- 
culations indicated the presence of a cubic phase 
with a lattice parameter of a = 4.14 kX, and this 
suggested the presence of a silicon hexaboride with 
a structure of the cubic hexaboride type. A com- 
parison of the x-ray photographs of polycrystalline 
CaBg and SiBg showed their similarity and con- 
firmed the suggestion that they were isomorphous. 

The intensities I and the positions of the lines 
were measured on x-ray photographs of samples of 
SiB,, and the value of d was subsequently calcu- 
lated (Table 1). 

From an exact determination the lattice period 
of SiBg was found to be a = 4.142 + 0.002 kX. The 
interplanar distances calculated on this basis (see 
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Table 1) were in excellent agreement with the meas- 
ured values. It should be noted that the period a for 
SiB, is very close to those of the known hexabor- 
ides (CaB,, LaBg, etc.) [3]. 

The density of the crystals of SiBg was deter- 
mined as Opyen = 2.15 g/ em?’ by hydrostatic weigh- 
ing [4] in carbon tetrachloride; with Z = 1, we ob- 
tained Ox-ray = 2.18 g/cm®. The density of SiB, 
is less than those of silicon and boron and this prob- 
ably has some effect on the conditions for obtaining 
the compound. In structures of the CaBg, type, the 
atoms take up the a and f positions in space group 
O}, — Pm3m [5] 


Si 1(a), 000, 
B 6), +(2+ 5). 


The structure amplitude is: 


Fax = fsi + 2f2 [cos 2m ( ha a = ) 


+ cos 2x (ha + “3 *) + cos2n (le + “3*)). 


A value of 0.21 was taken for x. The line in- 
tensities were calculated according to the formula 
I~Ly pF?, 


where Ly is the Lorentz polarization factor and p 
is the frequency factor. 
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TABLE 1, Data from Measurements and Calculations on the X-Ray Photograph 


of SiBg 
ee ee ee ee eee eee 
hkl | I I 
dale d exp | exp theor 
100 4.14 — = 0.1 
110 2.93 2.93 very strong 463 
111 2.39 2-40 strong 199 
200 2.07 DeO7 mod, strong Ae 
210 1.85 1.85 strong 166 
211 1.69 1.69 mod, weak 14,5 
re 1.46 1.46 mod. weak 15,8 
300 1.38 0, 
521 138 | 1.38 strong 16 "\ 216.7 
310 iret ioe strong 133 
341 120 d225 average 30 
222 1.19 — — 0.3 
320 TAS AS mod. weak 1250) 
321 1.108 1,107 average G20 
400 1,036 ANG35 mod, weak 30.9 
410 1,004 Oe) 
399 1° 004 \ 1.004 average 573} 60.8 
Alt 0.977 2255 
330 0.977 } 0.976 mod, strong 55.6} 78,1 
Bol 0.951 — — 1,4 
420 0,927 02927 mod, weak 42,7 
424 0.904 0.904 mod, strong &1 
332 0,884 0,884 mod, weak 36 
422, 0,845 0.845 very strong Oh? 
500 0.828 \ — | 2 ¢ 
430 0.828 ee 1.6 Be 
510 0.8123 101 \ 
431 0.8123 | 0.81245) strong 44g 7249 
Sylel 0.7972 5 134 
333 0.7972 | Dt0 72” | HAS SHORE 0.2} 134.2 


A comparison of the calculated and measured 
intensities (Table 1) showed excellent agreement. 

The samples studied contained, apart from the 
basic phase of SiBg, small amounts of impurities, 
probably SiO, and SiC, which produced additional 
weak lines on the x-ray photograph (2-5 lines). The 
difference in density obtained byus and in the litera- 
ture [2] probably resulted from the presence of dif- 
ferent quantities of impurities, denser than SiBg in 
the different samples. The isostructural hexabor- 
ides have a simple cubic lattice in which the boron 
atoms form a three-dimensional skeleton in the 
holes of which the metal atoms are found — silicon 
atoms in the case of SiBg,. The boron atoms are at 
the vertices of octahedra, and each boron atom is 
surrounded by five others at a distance of 1.72 kX. 
Four of these atoms are in the same octahedron 
while the fifth is in a neighboring octahedron. In ad- 
dition, each boron atom is surroundedby four metal 
atoms in such a way that the four metal atoms and 
the boron atom form a tetragonal pyramid with the 
metal atoms at the base. 


The metal atoms have a high coordination num- 


ber. They are surrounded by 24 boron atoms and 


six metal atoms. 
siderably greater than the sum of the atomic radii. 


This indicates the formation of a stable skeleton 


The metal—boron distance is con- 


from the boron atoms with the metal atoms distrib- 
uted in its holes. The metal—boron bond is weak- 
ened as the atomic radius of the metal decreases 
[6]. 

In the case of SiB,, the atomic radius of sili- 
con is considerably less than those of the other 
light metals which form hexaborides, and the sili- 
con—boron distance is very.different from the sum 
of the atomic radii (Si +B = 3.05 kx; Rsi,. + RB,» = 
1.34 +0.91 = 2.25). This peculiarity of the SiB, 
structure may explain the difficulty in making the 
compound SiB; [7]. 

The author thanks Professor G.S. Zhdanov 
(Head of the Department of Solid State Physics), 
and dozent M.M. Umanskii for discussion of the 
results, aspirant A.A. Stepanova for the gift of x- 
ray photographs of some hexaborides, and dozent 
G.V. Samsonov (Department of Rare Metal Metal- 
lurgy, Mintsvetmetzoloto) for presenting the 
samples. 
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Solution of the Structure 


The crystals of tetraphenylphosphonium iodide 
[(CgHs)4P]I were light yellow tetragonal prisms with 
dimensions of 4.0 x 0.3 x 0.3 mm.! Goniometry of 
the crystals showed the presence of faces of only 
one simple form {110} (relative to our chosen axis). 
The x-ray study was by the oscillation method and 
by photographing the reciprocal lattice in copper 
radiation. The lattice parameters were: 


a=b= 11.88 +-0.06A, 
V = 982 A3, 
M = 466.3, 


c = 6.96 +. 0.02 A (prism axis), 


p= 1.5 g/cm? Cray =svli93) ene" 21(4-89). 

By rotation about the c axis in a camera for 
photographing the reciprocal lattice, the zero-, 
first-, second-, and third-layer lines were de- 
veloped. The intensities of the reflections were es- 
timated visually and converted to a general scale 
via the oscillation photograph. The number of in- 
dependent reflections of which the intensities were 
measured was 231. In converting the intensities to 
structure factors, only the Lorentz and polariza- 
tion factors were taken into account. 

The symmetry of the reflections indicated that 
the crystals belonged to the Laue class 4/m, and 
the presence of systematic absences withh +k +l # 
2n indicated three possible space groups: 14/m, 
1/4, and I4. A distinct piezoelectric effect ruled 
out the first of these space groups; in the other two 
the I” and (CgHs)4P" ions must occupy special posi- 
tions on the symmetry axes when the multiplicity is 
two, but since the cation has a tetrahedral configu- 
ration, the choice of space group is unambiguous, 
namely, I4= St. This group is one of four which 


permit close packing of molecules with symmetry 
res oly 

Because of the absence of a center of sym - 
metry, the atomic coordinates were determined by 
calculation of a three-dimensional interatomic vec- 
tor series: 


P (ayz) = >) (Fin + Fiz) cos 2eha cos 2ky — 
— (Faxi — Fyx1) sin Zhe sin 2rky] cos 2xlz. 


In practice, a section was calculated along the c 
axis through the region of the maxima in a projec- 
tion on the ab face. (The side of the lattice was 
divided into 60 parts.) In this projection (independ- 
ent of the region a/4, b/2) the maxima next in 
height above the background were atx = 0, y= ee 
Consequently, the cation and’anion are on different 
4 axes; the iodine atoms occupy the 2 (a) positions 
:000, while the phosphorus atoms occupy the 2(c) 
positions: Le Ain Since the I — Pand I-C maxima 
are considerably larger than the P—P, P—C, and 
C—C maxima, only the maxima of the first two 
types are revealed in the interatomic function which 
appears as the union of two structures joined by a 
plane of symmetry through the iodine atoms and 
perpendicular to the 4 axis. Hence, the "untan- 
gling of the structure" is reduced to the choice 
of one of two values of the z coordinate for each 
maximuin (atom). The choice was made by tak- 
ing geometrical considerations into account (inter- 
atomic distances for P—C and C—C) and ap- 
peared to be unambiguous for all atoms. The coor- 
dinates found are given in Table 1. The coordi- 


1The authors thank G. A. Razuvaev for presenting the sample for 
study. 
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nates of the maxima were determined by interpola- 
tion, while those of the hydrogen atoms were cal- 
culated with the usual assumptions (C—H = 1.08 A, 
and the C—H bond is the bisector of the C -C —C 
angle). 

The structure amplitudes Fy,.9 were calculated 
from the coordinates of Table 1. They showed that 
all were positive and the corresponding centrosym- 
metric electron density projection (Fig.1) gave 
practically the same x and y coordinates 


Description of the Structure 


The structure of tetraphenylphosphonium iod- 
ide is constructed from I” anions and tetrahedral 
(C,H;)4P” cations which occupy special positions on 
the 4 axes of different sets (Fig. 2). 

We shall first discuss the structure of the ca- 
tion. The length of the P —C bonds is 1.80 A, EOS A 
0.07 A less than the sum of the single -bond coval- 
ent radii (rg +rp = 0.77 +1.10 =1.87 A). This 
shortening is not remarkable, since interatomic 
distances between an element and a carbon atom of 
a benzene ring are usually shorter than the sums of 
the corresponding single-bond covalent radii. It 
should be noted that there is only one recorded P—C 
bond length in the literature [2], and this was ob- 
tained from an electron-diffraction study of trim- 
ethylphosphine vapor, (CH3)3P, and is equal to 1.87 
A. The bond lengths in the benzene ring were: 


Chexi€ la Wee se 7= Cpodias 
Open Crver lett oC: Grud At 
Cea MABOE SC, bCe WAST 


Average 1.41-++0.02 A 


The equation for the plane of the benzene ring 
(nucleus IV of Fig.2) was found by the least-squares 
method and has the form 0.707x — 0.150y + 0.692z — 


Fig. 1. Electron-density projection on the ab face, The electron 
density was calculated in arbitrary units. The contour lines for 
iodine are about three times as far apart as in the remainder of the 
diagram. 
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—0.098 = 0. The distances of the carbon atoms from 
this plane are: 


G, 0.024, .C, O02 
Gl 038) Te, 420.03: 
Ce moni audi Ce panies 


Thus, the average divergence from the plane is 
+0.03 A and, consequently, the error in determin- 
ing the bond lengths is not greater than 0.05 A, 
Hence, the bond lengths in the benzene ring are the 
same within the limits of experimental error. For 
the error in determination of the bond angles, the 
divergence of the P—C, bond from the C;—C, diag- 
onal of the benzene ring should be taken; it did not 
exceed 2°, 

The bond angles at phosphorus were not tetra- 
hedral: ~ = C,(I)—P—©,(l) = 112.5°, w = C,(1)—P— 
C(I) = 103° (the numbering of the benzene rings is 
shown in Fig. 2). It appears that the divergence of 
the bond angles does not gainsay the 4 symmetry 
and it has a natural explanation, as is shown below. 
The angle between the planes of the benzene rings I 
and II is 118.5°, and that between the planes of rings 
I and II is 92.5°. Itis not difficult to show that the 
mutual orientation of the benzene rings of the cation 
in the crystal is real and corresponds to a minimum 
energy. 

We shall discuss the cation (CgH;),P* with sym- 
metry 4, Its configuration is described by the bond- 
ing angles ~ and w, and by the angle ~ between the 
normal to the plane of a benzene ring and the "ver- 
tical" plane which passes through the 4 axis and the 
P—C, bond (Fig.3). The presence of 4 symmetry 
requires that the angle ~ should be the same for all 
four benzene rings (I, Il, DI, IV) of the cation. 


A characteristic of the tetraphenylphosphonium 
ion is the steric hindrance between the carbon and 
hydrogen atoms of different benzene nuclei. For 
example, the distance between the "bound" atoms 
C, of two benzene nuclei with an ideal configuration 
about phosphorus is 2+1.80+sin(109°28'/2) = 2.04A 
(Fig.3), whereas twice the intermolecular radius 
for carbon is 3.4-3.6 A. For simplicity we shall 
assume that all the bond angles of the cation have 
the theoretical value, i.e., ¥ = w =109°28', and 
we shall find the relation of the steric hindrance in 
the cation to the angle of rotation ¢ of the benzene 
rings about the P—C bonds. This relation is shown 
in Fig.4. Here, A is the sum of the reductions in 
the distances between nonbonded atoms relative to 
the sum of the intermolecular radii: 


A= )\8n= >) Ri+ Re — dp, 


i,k isk 
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TABLE 1. Coordinates of the Symmetry-Independent Atoms 


A Relative 
Atom 
ae | y | be x | y | Zz 
I 0 (0) 0 0 0 0) 
P 0 5,94 1.74 0 Vo 1/4 
Cy 0.80 4.78 0.63 0.067 | 0,402 0.091 
(ep 0.48 3.40 0.59 0.040 | 0,286 0,085 
Cg 1.14 2.46 —0.22 0.096 0,207 —0,032 
Cg 2,06 2.98 —1.11 Ona: 0.251 —0,160 
G; eos 4.38 —1.07 0.200 0.369 —(),154 
Ce 170) Deo —0.26 0.143 0,447 —0.037 
H, =O 5 3.03 1.23 |—0.021 0,255 OBA 
Hg 0.92 1,36 —().20 0.077 0,115 —0,029 
Hy, 9.55 PAR PAY) —1,78 0.215 0.193 —-0,256 
H; 3.411 4,77 —1.71 0.262 0,402 —0,246 
Hs 1.89 6.44 —0.28 0.159 0,540 —0.040 
Center of grav - 153 0.94 —0,24 0,120 | 9.327 |—0.035 
ity of benzene 
ring 


Fig. 2. Projection of the structure on the face 
ab, A “skeletal” diagram is shown in the up- 
per half, together with the shortest distances 
between nonbonded atoms. The lower half 
shows the packing with the intermolecular ra - 
dii used in the calculations, 
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TABLE 2. Shortening of the Distance (A) Between Nonbonded Atoms in the Cation 


Se TG ee ee a >) ee 


"Real cation” 


"Ideal" cation "Ideal cation" 


== 119,5°, o = 103° = © = 109°28" » =o = 109°28', 
No. of ar se giacs mee ¥ > = 28° ouaae 
Distance distances 

dan ik diy 5% qi ein 
Ci) 1G ( 4 2.98 0.64 2.94 0.68 2.94 0.68 
Cie 4 BAT 4 | 0.454 0S oan 0-495 Hanesues = 
Cs) D) 4 3.29 0.33 3.13 0.49 3.19 0.43 
C5 (1) CD) 4 3.44 0.18 3,22 0.40 aay 0.14 
H,({) — C, (II) 4 2.81 0.06 2.60 0.27 2.70 0.17 
UDC.) 4 2.81 0,06 2.49 0.38 2.70 0.17 
C,(1) — C, (IL) > 2.8) 0.80 2.94 0.68 2.94 0.68 
CODES ONILY So 4 3.20 0.42 3.47 0.15 3227 0.35 
CD CAG) 2h 213.70 a eT eae 70) = 3.34 0.28 
Sum of the 
shortenings = — 8.96 — 11.60 — 9.56 


Fig. 3. Configuration of the cation, 
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Fig. 4. Relation between the steric hindrance 
and the angle of rotation of the benzene rings 
about the P—C bond. 


where Rj and Ry are the intermolecular radii of 
atoms i and k, and dj, is the actual distance be- 


tween them. 


Sums were taken over all pairs of atoms for 
which Rj + Rx > dj, (the values of the intermolecu- 
lar radii found in the given structure were used in 


the calculation; see below). 


It is shown in Fig.4 that minimum steric hind- 
rance corresponds to y = 52°. Consequently, this 
is the energetically most suitable angle of rotation 


of the benzene rings and should occur in the "free" 
cation with ideal tetrahedral bonding angles. The 
minimum will shift with changes in the bond angle. 
From the experimental results, ~ = 28° in the crys- 
tal. It is not difficult to see that additional decrease 
in the steric hindrance is caused by this change in 
gy if it is accompanied by a change in the bond 
angles. The shortest distances in the "ideal" and 
"real" cations are shown in Table 2. 

So the divergence of the bond angles in the real 
cation from the theoretical value is not fortuitous. 
While the minimum sum of the shortening is 9.56 A 
at ~ = 52° in the ideal cation, the sum is consider- 
ably less in the real cation: 8.96 A at gp = 28°, 

A.I. Kitaigorodskii [1] analyzed the configura- 
tion of tetra-aryl compound on the basis of less 
complete structural data, but his results are in 
qualitative agreement with ours. 

It is seen from Fig. 4 that rotation of the benz- 
ene nuclei in the cation (CgHs)4P* (and generally in 
tetra-aryl compounds of the general formula A,X) 
with retention of 4 symmetry (i.e., with "synchron- 
ous" rotation of all four nuclei) is not free. A con- 
siderable potential barrier must be surmounted. It 
can be shown that the potential barrier is retained 
if the rotation of the benzene nuclei is not synchron- 
ous (i.e., if the symmetry 4 is lost). In this con- 
nection, the introduction of a limited number of sub- 
stituents in the ortho- and meta-positions of the 
benzene rings might lead to optical isomerism. 

We shall now turn to a description of the pack- 
ing. The structure of tetraphenylphosphonium iod- 
ide must be described as close-packed layers of 
cations, with a coordination number of six, dis- 
tributed in the (110) planes [two cations are ob- 
tained from the initial one by the translation c, the 
other four by the translations +l, (a+b) + 6/2, or 
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TABLE 3. Shortest Intermolecular Distances 
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Dist, Py. | Diff He oh 
Atoms Bonds es adil. Atoms Bonds iy radii) 

A A 

A A 

(ly ee) ce | 3.63 | S628 0204 | Gy) a Hal e2r lh 2.00) 2870.08 
C;(I) — H, (II) @ | 2.84") 2-87 0.06 | Ha(Il)— Hay) 2, | 2.46 | 2.12 | 0.04 
Cie 6,10) e | 3.64 | 3,62 | 0.02 | G,(II) —Ca(iV’)| 2; | 3,68 | 3.62 | 0,06 
C4(1) — Hg(I1) ¢ | 2.83 | 2.87 | 0.04] H,(1l)—C,(IV’)| 2, | 2,94 | 2.87 | 0.04 

C;(I1) — H, (II) © 288" ees7 | ONO lec. =r = | Ae) = 


Fig. 5. Packing of the cations in a (110) layer and the adjacent anions. 


by the axes 2,;]. These layers (Fig.5) are related 
to one another by 4 and 2, axes. The layers are 
superimposed less densely than are the cations 
within the layers. The considerable holes left be- 
tween the layers are occupied by iodide ions. The 
coordination number of the cation is 14: six ca- 
tions in the same (110) layer, 2 cations from the 
layer underneath, 2 cations from the layer above, 
2 anions between the given and underlying layers, 


TABLE 4 
((C.Hs)PI]| (C,.H,)Si 
a 11.88.A|411,30A 
c 6.96A | 7,08 A 
Space group I4 P42,¢ 
n 2 2 
Symmetry = + 
of molecule 4 4 
(cation) : : 
X=C 1.80A]1.94A 
Packing 
coefficient] 0.76 0.74 


and 2 anions between the given and overlying layers, 
The coordination number of the anion is also form- 
ally 14 (10 anions and 4 cations) but it only touches 
the four cations. 

The shortest intermolecular distances in the 
structure, which correspond to atomic contacts are 
shown in Figs. 2 and 5 and cited in Table 3. 

The C—C distance leads to a normal carbon in- 
termolecular distance of 1.83 A and the H—H dis- 
tance gives a hydrogen radius of 1.08 A. Using one 
of these values to find the other from the C—H dis- 
tance we obtain 1.79 A for carbon and 1.04 A for 
hydrogen. Thus, the intermolecular radii in this 
structure are reasonably constant and are equal to 
1.81 A for carbon and 1.06 A for hydrogen. Accord- 
ing to Table 3, the average divergence from addi- 
tivity is 0.03 A. The value found for the intermole- 
cular hydrogen radius is 0.1 A less than the tabu- 
lated value. We believe that the tabulated value is 
too large because of the "compression" of hydro- 
gen in intermolecular contacts which has been fre- 
quently noted by us and by other workers (e.g. ,in 
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the structure of naphthalene [3]). The ionic radius 
of iodine in this structure is 2.22 A, which is some- 
what larger than the normally iceenien value (2.16 
A according to Pauling). Evidently, the value of 
PEIN may be taken as the radius of the free iodide 
ion since the polarizing effect of Pt can be neglected 
because it is so far from the I”. On the other hand, 
the values of the ionic radius given in the literature 
are determined from structures with monatomic 
cations — polarization of the iodide ion, which would 
shorten the interatomic distance, is not taken into 
account. 

The volume of the cation, calculated from the 
intermolecular carbon and hydrogen radii found 
above and the tabulated intermolecular radius for 
phosphorus (1.9 A), is 326 AS, the volume of the 
anion is 46 AS, so that the volume of the "mole- 
cule" is 372 A®, This gives a packing coefficient 
of 0.76. Thus, despite the "unsatisfactory" shape 
of the cation, the packing density is high enough. 

It is of interest to compare the structure de- 
scribed with the structure of other tetra-aryl com- 
pounds, e.g., with that of tetraphenylsilane 
(CgHs)4Si. The eryele of this compound belong to 
space group P42;c = S Var but the positions of the 
molecules are exactly the same as those of the ca- 
tions in tetraphenylphosphonium iodide [1]. 

The structures compared differ in the lengths 
of the X—C bonds and especially in the presence of 
large anions in the phosphorus compound. Never- 
theless, the similarity is very great.This indicates 
that there are considerable holes in the structures 
of tetra-aryl compounds, and these can be filled 
with anions on going to structures of the type 
[((CgHs),XTY~ without important changes in the 
character of the packing. The packing coefficient in 
[(CgHs5)4P]I is 0.66 without taking into account the 
anions. Thus, the "ready-made" holes between the 
filled (110)layers of molecules in the structures of 
tetra-aryl compounds are somewhat small for iod- 
ide ions. The introduction of the iodide ions causes the 
layers to shiftapart (the parameter a increases 
from11.30 to11.88 A), but this is compensated by 
filling of the holes (the packing coefficient increases 
from 0.71 to 0.76). Introduction of the anions can- 
not cause an increase in the parameter c. The in- 
crease in this from 6.96 to 7.08 A results from 
the increase in X —C bond length: 7.08— 6.96 = 
0.12 A *2(1.94-1. 80) sin (109°28'/2) = 0.15 A. 

We shall now discuss the hypothetical struc- 
ture (C,H;)4P with densest superposition of the 
layers. Shortening of the bond length from 1.94 
(Si-C) to 1.80 A (P—C) should lead to a decrease 
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in the parameter a= 11.30 A by 2(L. 94 —1.80) ° 
cos (109°28'/2)/V¥2=0.16 A. Thus, the parameter 
a in the hypothetical structure should be 11. ia Ae 
iiss Os (4 A less than the parameter a in the real 
structure of [(CgHs;)4P]I. Thus, the replacement of 
iodine with a small cation with a radius R= R; — 
0.74=1.48A requires no change in the structure 
of Ar,X with space group P42,c and it is to be ex- 
pected that the unchanged compounds (CgH;)4C 
(CgHs)4Si, etc., on the one hand, and tetraphenyl- 
phosphonium fluoride [(CgHs),P]F (rp- = 1.36 A) on 
the other, should be completely isostructural. The 
comparison carried out has justified the correct- 
ness of treating the structures of Ar,X and 
[Ar,X]*Y~ as the combination of densest layers of 
Ar,X particles with a coordination number of 6, the 
superposition of which is accompanied by a de- 
crease in density and the creation of holes which 
can be filled with anions or can remain vacant. 

The authors with to acknowledge the part taken 
in decipherment of the structure by L.A. Kamash- 
kina and K.F. Lugova. The authors express their 
thanks to Professor A.I. Kitaigorodskii for his at- 
tention and interest in the work. 


Summary 


1. The structure of tetraphenylphosphonium 
iodide has been determined by calculating the three- 
dimensional interatomic vector series. 

2. The following were determined: P—C = 
1.80 A, angle C-P—C = 103° and 112.5°, and the 
intermolecular radii of C, H, and Iare 1.81, 1.06, 
and 222 A, respectively. 

3. The configuration of tetra-aryl molecules 
and cations and steric hindrance and internal rota- 
tion in them have been discussed. 

4. The packing in the crystal of tetraphenyl- 
phosphonium iodide has been analyzed and com+ 
pared with that in structures of neutral tetra-aryl 
compounds. 


LITERATURE CITED 


1. A.I. Kitaigorodskii, Organic Crystal 
Chemistry [in Russian] (1955), p. 149. 

2. H.D. Springall and L.O. Brockway, "The 
molecular structures of the methyl derivatives 
of phosphorus and arsenic," J. Am. Chem. 
SoCrH, 60, 996 (1938). 

3. S.C. Abrahams, J.M. Robertson, and J. G. 
White, "The crystal and molecular structure 
of naphthalene," Acta Cryst., 2, 238 (1949). 


SOVIET PHYSICS ~ CRYSTALLOGRAPHY 


VOL li NO. NOV.-DEC., 1956 


THE STUDY OF POORLY-FORMED SINGLE CRYSTALS OF 
UNKNOWN SYMMETRY WITH AN X-RAY GONIOMETER 


E. Sandor 


Physics Institute, Budapest University 


Translated from Kristallografiya, Vol. 1, No. 6, 


pp. 677-680, November-December, 1956 
Original article submitted March 28, 1956 


A method is proposed for finding the direction of the axis of an important zone in poorly 
formed crystals of unknown symmetry, and of aligning this direction with the goniometer 
rotation direction. The zone axis is determined from a single Laue photograph prepared 


on a cylindrical film. 


All the unit-cell constants can be determined 
from oscillation photographs prepared on an x-ray 
goniometer, and from a few layer lines around one 
axis [1]. The problem arises here of determining 
the zone axis in poorly formed crystals and of 
aligning this direction with the goniometer rotation 
axis. In principle it makes no difference which zone 
axis is selected for the job, because any unit cell 
can be transformed by the appropriate methods in- 


to the Bravais cell characteristic of the lattice sym- 


metry of the particular crystal [2]. However, for 

practical reasons it is best to choose an important 
zone axis, if possible one corresponding to one of 
the crystallographic axes of the crystal. 


Determination of an Important Zone 
Axis 

The zone axis directions in poorly formed crys- 
tals are most simply found from Laue photographs. 
From a Laue photograph it is always possible to 
construct a stereographic projection of the planes 
giving the reflections, and this projection can then 
be used in the same way as the stereogram ob- 
tained from optical measurements of angles in a 
morphological study of a crystal. However, the 
treatment of an arbitrarily obtained Laue photo- 
graph is usually much more complex than that ofan 
"optical" stereogram, because of the following two 


factors: 


Fig. 1. Laue photographs obtained on a cylindrical film from a monoclinic crystal. The 
most important nodal points are circled and numbered. This group of points also in- 
cludes one which is forined by the intersection of zone curves outside the photograph. 
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Fig. 2. Stereographic projection of the 
points marked by circleson Fig.1. 


(a) If the x rays do not fall on the crystal along 
the direction of an axis of symmetry or normal to 
a plane of symmetry, then it is impossible to see 
the crystal symmetry on the stereographic projec- 
tion. 

(b) Quite a large proportion of the Laue spots 
are high-order reflections, which makes analysis 
of the stereographic projection even more difficult. 

A systematic method of overcoming these dif- 
ficulties has been worked out by Gross [3]. The 
basis of his method is that a series of flat Laue 
photographs, selected in a certain way, will give a 
unique result for the symmetry system (or rather 
the Laue class) of the crystal, which establishes 
the directions of the crystallographic axes. How- 
ever, even in the simplest case this means prepar- 
ing at least 16 photographs, and for crystals of 
higher symmetry even more are required. Kvitka 
and Umanskii have subsequently shown [4] that itis 
possible to choose important plane networks and 
zones in a Simple manner from the Laue photo- 
graphs and use these to establish the directions of 
the crystallographic axes and also the indices of the 
most important spots. This method requires only 
three suitably chosen Laue photographs. 

A method is described below, by means of 
which it is possible to determine the axis direction 
of an important zone from only one Laue photo- 
graph, obtained on a cylindrical film in an x-ray 
goniometer. The exposure time needs to be chosen 
so as to give clear images of the spots with the 
highest possible 3 (Fig.1). On the photograph, 10- 
15 spots are chosen which lie at the intersections 
of the largest possible numbers of zonal curves; 
these will be the spots with the largest empty areas 
around them. With the aid of the appropriate (S,®) 
diagram and a stereographic net, these "nodal 
points" can be transferred to a stereographic pro- 


Fig, 3. Screening cylinder with a wide 
aperture, used in obtaining four Laue 
photographs on an x-ray goniometer film. 


jection (Fig.2). 

Immediately after this, a check is made to see 
whether mirror planes can be detected on the 
stereogram or not. Since the normal to a mirror 
plane is simultaneously a zone axis, on the stereo- 
graphic projection a mirror plane can be recog- 
nized because, on the zone circle corresponding to 
it, there appear the images of points belonging to 
other zones. With a stereographic net it is easy to 
find the nodal points, which are the poles of zones 
passing simultaneously through images of points 
belonging to other zones. Each of these "nodes" is 
the stereographic projection of a mirror plane. 

The shaded dumbbell-shaped area in Fig. 2 
does not contain points because the corresponding 
spots were lacking on the photograph. This area is 
so small that it generally will not interfere with the 
treatment of the Laue photograph. In this particu- 
lar case, for example, it is easy to see that the 
zonal curve, which has its pole at the No. 1 point 
which did not appear on the photograph, is a mir- 
ror plane, with the circle corresponding to it on 
the stereographic projection passing through the 
points 10, 12, and 138, respectively. Since the 
photograph was obtained without a beam catcher for 
the primary beam, there will actually be a narrow 
annular region around the edge of the stereographic 
projection where no spots will fall. This region 
has not been drawn in on the stereographic projec- 
tion, because, by using a suitable beam catcher it 
could be made very small. 

A search for mirror planes can, in principle, 
lead to the following four alternatives: 

(1) no mirror planes observed on the projec- 
tion; 

(2) one mirror plane found; 

(3) several mirror planes found, all lying at 
similar angles to one another (30, 45, 60, or 90°); 
or 

(4) several mirror planes found, with angles 
of 30 or 45° between them, and, in addition, at 
least one more mirror plane perpendicular to these 
planes. 

In the first case it can be assumed that the crys- 
tal belongs to one of the two Laue classes withouta 
mirror plane (L or 3). 


POORLY-FORMED SINGLE CRYSTALS OF UNKNOWN SYMMETRY 


i ie 7 | WT 


Fig. 4. Rule for measurement of quadruple Laue photographs. 


In accordance with this, an important zone is 
chosen on the Laue photograph (i.e., the line with 
the highest number of spots on it and with the wid- 
est empty space on either side of it), and the axis 
of this zone is lined up with the axis of rotation of 
the goniometer. This approach will usually lead to 
a triclinic unit cell, which can always be converted 
into a unit cell characteristic of the particular sym- 
metry system by using the well-known Delaunay 
methods. If the crystal should happen to belong to 
the trigonal Laue class 3, this will be shown up dur- 
ing the conversion. 

In the second case, the crystal obviously be- 
longs to the monoclinic Laue class 2/m, when it 
is best to line up the goniometer rotation axis with 
the normal to the single mirror plane. 

In the third case, according to the interplanar 
angle shown, it may be assumed that the crystal 
belongs to the Laue class 6/m, 4/m, or 3/m, re- 
spectively, or,when the angles are equal to 90°, to 
the class mmm or ms. In this case, it is best to 
choose the goniometer rotation axis to coincide 
with the line of intersection of any two planes. 

Finally, if, in addition to the planes lying at 
angles of 30 or 45° to each other, there is a fur- 
ther plane perpendicular to these, the crystal can ~ 
only belong to the Laue class 6/mmm, if the angle 
is 30°, or the classes 4/mmm or m38n, if the 
angle is 45°. In this case it is also best to have the 
goniometer rotation axis lined up with the line of 
intersection of any two planes making an angle of 
30 or 45° with each other. 

Where the above method is concerned, the fol- 
lowing additional points are worth noting: 

(1) Since the overwhelming majority of crys- 
tals belong to Laue classes possessing at least one 
plane of symmetry, in most cases a check for mir- 
ror planes will give a positive result. 

(2) Even though Laue photographs obtained on 
a crylindrical film include a spot-recording area 
extending over almost all of the stereographic pro- 
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jection, it may nevertheless happen that the reflec- 
tion corresponding to a mirror plane will be absent 
from the Laue photograph. An example of this can 
be seen in Fig. 1, where the reflection correspond- 
ing to the only mirror plane of the monoclinic crys- 
tal (the point of intersection of zonal curves lying 
outside the photograph) is missing from the Laue 
photograph. 

However, on the stereographic projection it is 
easy to establish that this particular missing spot 
corresponds to a mirror plane. 


Aligning the Chosen Axis with the 
Axis of Rotation of the X-ray 


Goniometer 


After the position of an important zone axis has 
been established, it can be lined up with the axis of 
rotation of the goniometer using the goniometer head 
arcs. However, since the graphical treatment of the 
cylindrical Laue photograph cannot be expected to 
be accurate to more than 1-2°, the accuracy of the 
resulting setting will generally not be great enough 


Fig. 5. Quadruple Laue photograph, to refine the setting of 
a monoclinic.(NayS,O3) crystal along the normal to the mir- 
ror plane. a) Photograph in an imaccurate setting; b) control 


photograph from an accurately positioned crystal. 
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to satisfy the requirements for a goniometer photo- 
graph. In practice, the positioning is complicated by 
the fact that the required angle of rotation very often 
cannot be achieved through the use of the goniom- 
eter head arcs alone, so that the crystal can be 
brought to the desired position only through apply- 
ing other auxiliary methods (with a small universal 
joint, etc.), or by regluing, which reduces the ac- 
curacy of the setting even further. The rough set- 
ting obtained from the stereographic projection must, 
therefore ,be further refined. 

On an x-ray goniometer this refinement is most 
simply brought about by the four-photograph Laue 
method (see Winchell [5]), in which four Laue photo- 
graphs are obtained side by side on a cylindrical 
film. The first photograph is taken in a position 
such that the direction of the x-ray beam coincides 
with the plane of one of the goniometer arcs; the 
other three photographs are taken after rotating the 
goniometer axis through 180, 90, and 270°. If the 
zone axis coincides exactly here with the goniom- 
eter axis, the corresponding zone on the four photo- 
graphs will be parallel to the equatorial line, andif 
the chosen zone axis is a normal to a mirror plane, 
then all four photographs will also show a mirror 
symmetry plane (Fig. 5b). Two difficulties arise 
in the use of this method: 

(a) the four Laue photographs partially over- 
lap one another, which greatly complicates their 
treatment; and, 

(b) on the equatorial-zone curves the spots are 
sometimes so far apart that it is quite difficult to 
measure the distances between curves when # = 90°. 

The first difficulty is most simply overcome 
by using a screening cylinder with the film, like 
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that used for obtaining layer-line photographs, but 
with a wider aperture; this cylinder will leave free 
only a quarter of the film at a time (Fig. 3). 

It is possible to speed up the measurement of 
spots with ¥ = 90° by using a suitable millimeter 
rule, drawn on glass or Plexiglas, which is laid on 
the quadruple Laue photograph (Fig. 4). 

Figure 5a shows one of these Laue photographs, 
obtained from a monoclinic crystal (anhydrous so- 
dium hyposulfate, Na,S,O3) using a wide-aperture 
screening cylinder. The zone axis chosen is the 
normal to the mirror plane. From the control 
photograph (Fig. 5b) it can be seen that even in the 
unfavorable case where the original position was 
quite inaccurate, and the equatorial zone spots lay 
some distance away from each other, the method 
described here leads to asatisfactory result. 
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A study has been made of changes in morphological symmetry for all classes of crystal 
during possible ferroelectric phase changes. 


if Introduction 


At certain temperatures, the Curie points, 
ferroelectric crystals undergo a sharp change ina 
number of their physical properties. These changes 
in properties are due to the fact that a phase change 
of the first or second kind occurs at the Curie 
points. Some ferroelectric crystals have several 
ferroelectric phases, and thus have several phase 
transition points. The Curie point usually refers to 
the temperature at which the crystal goes over 
from the nonferroelectric to the ferroelectric state 
(or vice versa). The temperature of transition from 
one ferroelectric phase to another is usually called 
simply a phase transition point. In what follows, by 
a ferroelectric phase change we will mean either 
type of transition. Incidentally, it should be noted 
that ferroelectric phase changes are reversible. 
For simplicity in what follows, unless otherwise 
stated, we will consider changes in one direction 
only, but this does not reduce the generality of the 
arguments used. 

Changes in the physical properties of ferro- 
electric crystals during ferroelectric phase changes 


are always accompanied by microstructural changes. 


as structural investigations show, so that the space 
symmetry group of the crystal alters. When an 
area in a ferroelectric crystal changes from a non- 
ferroelectric to a ferroelectric phase, a domain 
structure arises in that area. In cases where a 
crystal alters from one ferroelectric state to an- 
other, the nature of the domain orientation changes. 
It is important to note that experiments show 
that the morphological (point) symmetry of a do- 
main differs from the symmetry which the crystal. 
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(or domain) had before the phase change. In other 
words, the structural changes occurring during a 
ferroelectric phase change in a ferroelectric crys- 
tal are accompanied by a change in the point sym- 
metry. 

Thus, for example, a crystal of Rochelle salt 
belongs to the orthorhombic class 2:2 above a 
temperature ot 24.5°C (the upper Curie point); in 
the temperature interval between the Curie points 
(-18.5°C to +24.5°C) the crystal domains belong to 
the monoclinic symmetry class 2. Similar sym- 
metry changes can also be observed in other ferro- 
electric crystals. 

A careful scrutiny of the fundamental changes 
taking place during phase transitions in ferroelec- 
tric crystals leads to the conclusion that the 
changes in crystal or domain symmetry are gov- 
erned by certain definite laws, from which it is 
possible, for example, to accurately foretell how 
the point symmetry of a crystal belonging to a given 
class will alter if the crystal passes through a fer- 
roelectric transition, how many transitions can oc- 
cur in crystals of a given class, etc. The present 
article is devoted to a discussion of these ques- 
tions. 


2. Symmetry Changes during Ferro- 
electric Phase Transitions 


The most important characteristic feature of 
a phase change in a ferroelectric is the creation, 
or disappearance, or change in direction, of spon- 
taneous polarization in the crystal. This feature al- 
so governs changes in other properties of ferro- 
electrics, i.e., electrical, optical, mechanical 
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properties, etc. Spontaneous polarization in acrys- 
tal is defined by a polar vector Pg, which has the 
symmetry ©-m. The individual domains in a fer- 
roelectric crystal can,therefore,only belong to one 
of the pyroelectric classes which are subclasses of 
the class © +m, which means only the following: 
Tere 


? 


3d, 4, 6, m, 2-m, 3-m, £m, 6-m. 


The change in point symmetry which takes 
place during a ferroelectric phase change can be 
found in two ways. The first of these uses the 
Curie symmetry principle [1]. According to this 
principle, the symmetry which a crystal attains 
when spontaneous polarization arises in it is given 
by the highest subgroup of the symmetry group of 
the crystallographic class which is common both to 
the crystal under consideration, when in the non- 
ferroelectric phase, and to the symmetry group of 
the polar vector .. for a given arrangement of the 
symmetry elements of both groups. 

In the general case, spontaneous polarization 
may be assumed to arise in a crystal in any direc- 
tion. Since there is an infinite number of different 
directions in a crystal, then formally each crystal 
can have an infinite number of physically different 
ferroelectric phases. However, we will limit our- 
selves to considering phases which differ in their 
symmetry. It is easy to see that the number of 
phases which differ in symmetry is not great and 
may be calculated simply. 

As an example, let us consider a crystal be- 
longing to the class m-2:m. If spontaneous po- 
larization arises along one of the two-fold axes in 
this crystal, its symmetry drops to 2-m; if the 
vector P, lies in a plane of symmetry and does not 
coincide with any of the symmetry axes, the crys- 
tal symmetry is lowered to the class m; and, final- 
ly, if the vector Pg does not lie in any of the sym- 
metry planes, the symmetry of the crystal is re- 
duced to the class 1. 

The second method of determining a change in 
symmetry is more graphical, and comes down to a 
direct examination of the uniform spontaneous de- 
formation. This deformation is piezoelectric and 
electrostrictional in nature, but in crystals which 
do not have piezoelectric properties in their initial 
phase (whether in general or along the direction of 
the vector Pg), it is only electrostrictional. Here, 
the character of this deformation can be found from 
an examination of the electrostrictional tensor or 
the piezoelectric constant tensor. It is easy to see 
that in a crystal belonging to, say, the class 2:2, 
spontaneous polarization arising along the X axis 
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will lead to shear deformation about the X axis (in 
accord with the form of the piezoelectric constant 
tensor), which will result in the crystal becoming 
monoclinic, with the symmetry class 2. (This 
example corresponds to the ferroelectric phase 
change in Rochelle salt.) If the direction of the 
vector P, does not coincide with any of the two- 
fold symmetry axes, the crystal will undergo shear 
deformation about all three axes, which will result 
in the symmetry being lowered tosymmetry class 1. 

For crystals not possessing piezoelectric prop- 
erties in their initial nonferroelectric phase, the 
change in symmetry can be found in a similar man- 
ner from the electrostrictional tensor. There isno 
need to examine electrostrictional deformation in 
crystals possessing piezoelectric properties in 
their initial nonferroelectric phase, because the 
electrostrictional deformation leads to the same 
symmetry changes as the corresponding piezoelec- 
tric deformation. 

Thus, by using either the first or the second 
method, it is possible to find the symmetry change 
occurring during ferroelectric transitions in any of 
the 32 crystal classes, i.e., when the vector Py 
arises or changes its direction. Tables 1-4 give 
the results of such an examination for crystals of 
all classes.} 

Tables 1-4 show symmetry changes for all 
classes during possible ferroelectric phase changes. 
In other words, it does not follow from the tables 
that a crystal of a given class will necessarily 
undergo the transformation indicated when the tem- 
perature is altered, because it may either not 
undergo any transformation at all, or it may be 
transformed without spontaneous polarization aris- 
ing (a nonferroelectric transformation). However, 
it does follow from the table that if the crystal does 
undergo a ferroelectric phase change when the tem- 
perature is altered, then this transition can only 
lead to the symmetry change shown in the table. It 
may be observed here that all known ferroelectric 
transitions conform completely with these tables. 


3. Discussion of Some Possible 
Ferroelectric Transitions. Orienta- 
tion of Domains 


1. For crystals belonging in their initial 
phases to one of the pyroelectric classes, two al- 


1A1l these examinations were essentially based on one assumption, 
that the change in crystal symmetry during a phase change was due 
only to the formation of spontaneous polarization, This assumption 
is quite natural and, as can be seen from examples, is supported by 
all known ferroelectric transitions. 
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ternatives are formally possible. The first alter- 
native relates to those crystals which, although be- 
longing to the classes in question, are not spon- 
taneously polarized in their initial states. In this 
case the formation of spontaneous polarization along 
the existing special polar direction cannot be con- 
sidered as a ferroelectric transition, because, in 
view of the uniqueness of this special direction, the 
spontaneous polarization cannot arise with equal 
probability in any other direction, and,consequent- 
ly, the crystal cannot undergo a change in symmetry 
and cannot be split into domains; splitting into do- 
mains is a necessary condition for a ferroelectric 
transition. 

There may also occur the other alternative, 
when the crystal has spontaneous polarization inits 
initial state, but is not split up into domains, i.e., 
it is a simple pyroelectric. It could theoretically be 
assumed that such a crystal would, at a certain 
temperature, undergo an abrupt change in the ab- 
solute value of its spontaneous polarization (without 
changing its direction). It is easy to see that this 
would not be a ferroelectric transition, since the 
direction of the vector Pe would still be unique, and 
the crystal could not split up into domains (and the 
crystal symmetry also could not change). 

Thus, the transitions of the 6-m~>6-m,4:m 
~<+4-.m type listed in Tables 1-4 cannot be con- 
sidered as ferroelectric transitions, and they are 
only included in these tables because they resulted 
from the complete and rigorous derivation proced- 
ure used. 

As well as these transitions, the pyroelectric 
classes listed in Tables 1-4 may formally be trans- 
fofmed into the classes m and1. Such transitions 
are not met with in practice and, moreover, in our 
view, would be unrealistic. It can,therefore,be 
considered that if a crystal from a pyroelectric 
class does not split up into domains, it cannot 
undergo ferroelectric transformations. 

It should be noted here that the above observa- 
tion regarding unrealistic transitions to classes m 
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and 1 does not relate solely to pyroelectric classes, 
but to all classes. 

2. The conclusion arrived at in the preceding 
section, that none of the pyroelectric classes could 
be the initial class for a series of ferroelectric 
phase changes, in no way contradicts the known 
fact that some ferroelectric crystals can pass from 
one ferroelectric phase into another, because, in 
the latter case, not only does the magnitude of the 
spontaneous polarization change, but also its direc- 
tion. A change in polarization like this is always 
bound up with a change in symmetry of the crystal. 
Examples of such transitions are those of BaTiO; 
at 45°C, from a tetragonal to an orthorhombic 
phase, and at —70°C, from an orthorhombic to a 
rhombohedrai phase. 

These transitions are generally possible be- 
cause, in its initial nonferroelectric phase, the 
crystal has a symmetry such that for differentdi- 
rections of the vector Ps (along simple crystallo- 
graphic directions) the initial symmetry can alter 
in different ways. It thus turns out to be possible 
to have several adjoining ferroelectric phases for a 
single initial nonferroelectric phase. (In BaTiOs, 
for example, the possible transitions are 6/4— 
4-m, 4-m+2-m, 2-m~+3-m, etc.) 

Because of this feature, it may be asserted 
that the character of symmetry changes in any fer- 
roelectric transition is determined, not by the 
symmetry possessed by the domains before the fer- 
roelectric transition occurs, but by the symmetry 
which the whole crystal has in its initial nonferro- 
electric state; moreover, the spontaneous deforma- 
tion in any transition is described by the electro- 


*The order in which the ferroelectric phases follow one another is 
not strictly constant within a particular trausformation series (for 
example, the perovskite solid-solution ferroelectrics may show the 
transitions 6/m-* 4: m= 2->m™ 3: mor6/m=™ 4: mor 
6/4 > 3: m, etc., for different compositions), Because of this, 
the order cannot be forecast in the tables of possible transitions for 
auy particular series, 


TABLE 1, Transitions from Cubic Classes 
Set of P, Initial phase 
directions ; 9 ie 
6i4 ele 3i4 6/2 3/2 
Possible ferroelectric phases 
<100> 4-m(b) 4(6) 2-m(6) 2-m(6) 2(6) 
<i 3-m/(8) 3(8) 3-m(4) 3(8) 3(4) 
<110S 2-m(12) 2(12) m(12) Lm(12) 
<hk0> 1(24) i 
chkk> |) -m(24) see m(12) 1(12) 
<hhk> (24) |p 1(24) | 
<hkl> 1(48) 1(24) 


TABLE 2, Transitions from Hexagonal and Rhombohedral Classes 


Initial phase 
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Note. 


Symbols are given for Bravais conditions (see, for example, [4]). 


strictional and (where applicable) the piezoelectric 
deformations of the initial, nonferroelectric phase. 

3. If spontaneous polarization arises along a 
particular direction in a ferroelectric crystal, then 
it could arise with equal probability along other 
crystallographically equivalent directions. Ina 
multidomain crystal, therefore, the domain polari- 
zation directions are found to be equally distributed 
among the complete set of these crystallographical- 
ly equivalent directions. 

The number of spontaneous polarization direc- 
tions n for each transition may be found from the 
ratio of the order of the symmetry group of the ini- 
tial nonferroelectric phase N;, to the order of the 
symmetry group of one of the domains Np, arising 
on transition: 


(i ING SING: 


Thus, for example, when BaTiO; transforms 
from class 6/4 (symmetry group order 48) toclass 
4+m (group order 8) the number of orientations of 
the vector P. will be 6, while for transition to 
class 3-m (group order 6), the number of Pg orien- 
tations will be 8, etc. In Tables 1-4 the numbers of 
orientations are given in brackets for each transi- 
tion, after the symbol for the class into which the 
crystal is transformed. 

Thus, the general character of domain orienta- 
tion in a multidomain ferroelectric crystal is de- 
termined by the symmetry of the initial nonferro- 
electric phase and the Ps direction in the given 
ferroelectric phase. 

4. The ordering of domain orientations along 
crySstallographically different directions of spon- 
taneous polarization leads to the position where the 
multidomain crystal is in essence a polysynthetic 
twin (or a special texture specimen with axes of 
finite order), and the domains are the components 
of this twin. The twinning elements are formed by 
the loss of the symmetry elements present in the 
initial nonferroelectric phase when it goes over in- 
to the ferroelectric phase. (This holds for any fer- 
roelectric phase. ) 

Here ,the multidomaincrystal, with respect to 
statistically averaged macrosymmetry, consequent- 
ly,reverts to the symmetry of its initial nonferro- 
electric phase.? This applies not only to its morph- 
ological symmetry, but also to all its physical 
properties when taken as a whole (its macroscopic 


3This fact, noted for Rochelle salt by Chernysheva [2], has recent- 


ly been examined in detail by us [8]. 
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TABLE 3. Transitions from Tetragonal Classes 
oe 


Set of P, Initial class 
directions 
m-4:m 4-m | 4:m 4:2 | 4 am | ve 
Possible ferroelectric phases 
<001> 4-m(2) | 4-m(1) 4(2) 4(2) 4(1) 2-m/(2) 2(2) 
100 
et oS. \ 2-m(4) \ m(4) Il m(4) \ 2(4) \ 2(4) | 
<hkO> 1(8 
<h0lS | m/(8) is as 1(A) \ 1(8) La) 
<hhlS ( | 1(8) (8) m4) 
<hki> 1(16) 1(8) j 1(8) | 
TABLE 4. Transitions from Orthorhombic, Monoclinic, and Triclinic Classes 
Set of Ps Initial class 
directions 3 ee 
me2:m 2-m | 232 | 2m | m 2 | 2 | of 
Possible ferroelectric phases 
rie hoe 2-m/(1) 2(2) m/(1) 2(4) 
m 2(2 de KP, 
<100S ( } m(2) (2) ; m(2) mit) 
<hk0> |  7(A) 1(4) 1(2) ; 12 1(1 
<hOl> m/(4) ny m(2) 1(4) mi) m/(1) Az) eC : 
<Okl> f \ 1(4 \ 12 
<pk> | 18) | 14), pee Pe: 
| | : 


properties); this has been rigorously proved by us 
for the example of Rochelle salt [3]. 

In particular, if a crystal belongs to a non- 
piezoelectric class in its initial nonferroelectric 
state, then the multidomain crystal (e.g. , BaTiOs) 
will not possess piezoelectric properties in any 
ferroelectric phase. If the initial nonferroelectric 
phase belongs to a piezoelectric class, then in any 
ferroelectric phase the multidomain crystal will 
show piezoelectric properties (e.g., Rochelle salt). 

It is natural that the conclusion concerning re- 
version to the symmetry of the initial phase is only 
statistical, as an average. If for some reason or 
other the domains in a multidomain crystal had a 
certain preferred orientation, this conclusion would 
not be valid. This would apply, in particular, if we 
were dealing with a specimen cut from a particular 
growth pyramid, rather than with a whole crystal. 

5. Analysis of Tables 1-4* reveals that not all 
the formally possible ferroelectric phase changes 
shown in these tables are equally possible. Some 
have specific features which have not been met with 
up to now, and these changes appear unlikely. (We 


noted above the formal nature of the series of trans- 


formations with initial pyroelectric classes, and of 
transitions into classes m and1.) In fact, several 
of the transitions included in Tables 1-4 (6/23, 
6:m <6; 6-m +2; 3:m +38, etc.) are such that 


the reversion of a multidomain crystal to the sym- 
metry of the initial class (see section 4) can only 
occur under conditions where the crystal in a fer- 
roelectric phase contains equal numbers (twin com- 
ponents) of left-handed and right-handed enantio- 
morphs. Such transitions do not seem to be very 
probable. 

In the transitions m:3:m <2°m, 3;2 ~2, 
3/2 +3, the vector Ps is directed along the polar 
axis of the initial phase, so that in this case an 
antiparallel arrangement of domains is not possible. 
In the transitions 6 +2) <=6) ms3'sim 3) imy 3 +23, 
etc., although the initial phases belong to piezoelec- 
tric classes, the spontaneous deformation is elec- 
trostrictional because, in these cases, the vector 
Pg is directed along the Z axis of the cell, a direc- 
tion in which the initial phase does not possess 
piezoelectric properties. In spite of this, itappears 
to us that in crystals with averaged symmetry a 
ferroelectric transition with a Bs vector directed 
along a high-order axis is more favorable than 
transitions with other Pe directions. 


4It is pointed out that these tables may also be used to determine 
the symmetry changes in a crystal when any influence described by 
a polar vector is applied to it, Thus, with the aid of the tables, it . 
can be easily seen how the crystal symmetry will change if an elec- 
tret is obtained from it, or from what crystal symmetries and 
through what action a Kerr cell may be obtained, etc. 
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4. Conclusions 


From Tables 1-4, which are in complete agree- 
ment with all known ferroelectric phase transitions, 
it is possible to make a series of predictions, which 
are quite unambiguous when the low probability (or 
even impossibility) of some of the formally possible 
transitions is taken into account. 

From Tables 1-4 it is possible to make an un- 
ambiguous forecast, for a particular ferroelectric 
belonging to a given series of transformations, of 
the absence or presence of new ferroelectric transi- 
tions, i.e., the number of possible ferroelectric 
phases differing in symmetry can be predicted. In 
order to discover whether these possibilities are 
realized in practice, it is, of course, necessary to 
carry out the appropriate studies in each particu- 
lar case. 


ZHELUDEV AND L. A. SHUVALOV 


The authors are deeply grateful to Academician 
A.V. Shubnikov for much valuable advice and for 
discussing the manuscript. 
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The characteristics of the behavior of Rochelle salt domains in alternating electric 
fields of a strength of about 50 V/cm and at frequencies of 32 and 1000 cps are de- 
scribed with reference to microcinematographic frames. Hysteresis of the reorienta- 
tion of the domains in alternating fields is followed by means of microcinematograph- 


ic frames. 


It is generally accepted to associate the polari- 
zation processes of ferroelectrics with reorienta- 
tion of the domains. Such a connection may be based 
on the observed variations of domain structure when 
constant electric fields are applied to ferroelec- 
trics. Considerable interest, however, attaches to 
the behavior of the domain structure in alternating 
electric fields. For this purpose, we have made a 
cinematographic study of the domain structure of 
Rochelle salt when electric fields of a frequency of 
32 and 1000 cps were applied to the crystal. 

For the cinematography of even comparatively 
slow processes, high-speed filming is necessary. 
However, we used another method, the ordinary 
filming method (about 32 frames per second), but 
using stroboscopic illumination. By the use of 
stroboscopic illumination, the observation of rapid- 
ly ensuing processes can be extended in time. The 
possibility of such an observation is evidence of 
the fact that for given field strength and given fre- 
quency, the process of movement of the domains 
takes place reproducibly from cycle to cycle. 

In filming at the frequency of 32 cps, the 
stroboscopic effect was obtained from the differ- 
ence between the frequency with which the camera 
objective was closed by the shutter and the frequen- 
cy of the electric field. The resulting slowing down 
of about 40 times made it possible to examine the 
polarization revetsal process as if occurring in 
one cycle of the electric field. For the frequency of 
1000 cps, a mechanical stroboscope was used, by 


543 


means of which filming took place with a slowing 
down of about 4000 times. Filming was carried out 
at the different frequencies in a microcinemato- 
graph camera especially constructed for the pur- 
pose. The specimens used were plates of X-cut 
Rochelle salt, 20 x 20 x 1mm in size. The elec- 
trodes on the specimens were Silver leaf, but in the 
observation area, instead of the silver coating, a 
transparent conducting layer about 1 mm? in area 
was applied. The specimen was thermostatically 
controlled during observation. Observation was car- 
ried out at a temperature of 11.5°C. 

Figure 1 shows the variation in domain struc- 
ture at a frequency of the electric field of 32 cps 
and a field strength of 36 V/cm. Concurrently with 
the observation and cinematography of the domain 
structure, observations were made on the charac- 
ter of the polarization of the specimens using an 
oscillograph and the usual circuit for recording 
hysteresis loops. In accordance with this, the 
frames selected to illustrate the variations in do- 
main structure were arranged on a curve corre- 
sponding qualitatively to the curve observed on the 
oscillograph. 

As will be seen from Fig.1, the polarization 
curve is far from saturation for the given field 
strength, but the dependence of polarization hys- 
teresis on the field strength is fairly clearly shown 
on the curve. In addition to the polarization curve, 
the figures also show the sinusoidal variation of the 
electric field with time (in the present case, the 
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Fig. 1. Variation in domain structure of Rochelle salt at a frequency of 32 cps. The Roman numerals on the field variation 
curve correspond to the Arabic numerals (frame numbers) on the polarization variation curve. The amplitude value of the 
electric field strength, 36 V/cm. 


Fig. 2, Variation in domain structure at a frequency of 1000 cps. The Roman numerals on the field variation curve cor- 


respond to the Arabic numerals (frame numbers) on the polarization variation curve. Amplitude value of the electric field 
strength, 56 V/cm. 
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time axis is the axis of ordinates). The points 
shown in Roman numerals on the sinusoid corre- 
spond to the frames shown on the polarization curve 
and marked with the corresponding Arabic numer- 
als (the point on the sinusoid corresponding to 
frame 3 is underneath the frame situated in the cen- 
ter of the figure). 

In the center of the figure is the frame corre- 
sponding to the domain structure of the portion 
filmed in the absence of a field. As will be seen 
from this frame, the domains of opposite orienta- 
tions in the limits of the filmed portion have un- 
equal areas. The part corresponding to the light 
component occupies the predominant area, while 
the part corresponding to the dark component oc- 
cupies an insignificant area. Thus, the domain 
structure in this portion is unipolar. It should be 
noted that unipolarity of the domain structure in 
the limits of the observed region is not evidence of 
unipolarity of the specimen as a whole. 

The frames located on the polarization curve 
correspond to field strengths differing from each 
other in phase by 45° (0°, 45°, 90°, 135°, etc.). 
Frame 1 corresponds to the amplitude value of the 
field strength. It may be seen from frame 1 that 
the light component of the domains occupies prac~ 
tically the whole field, while the dark component, 
antiparallel to the field, is practically absent. Re- 
duction in the field strength from the amplitude 
value to zero (frames 2 and 8) results in the forma- 
tion of the dark component, but its proportion is 
still insignificant. 

Increase in the field strength of opposite sign 
already leads to a sharper increase in the propor- 
tion of dark component (frames 4 and 5). The fact 
that at the amplitude value of the field (frame 5) 
the dark component does not occupy the entire field 
of view (in contradistinction to what occurs at the 
amplitude value of the field strength of reverse po- 
larity — frame 1), is understandable, since, in the 
initial state (frame in the center of the figure) in 
the field of view the dark component was quite weak- 
ly pronounced. The reduction in the field (from 
frame 5 to frame 7) does not result in any substan- 
tial change in domain structure; it is as if thecrys- 
tals "remembered" the state of domain structure 
which existed at the amplitude value of the field 
strength (frame 5). Increase in the field strength, 
illustrated by frames 7-1, shows the more pro- 
nounced variation in the domain structure in this 
portion. In this respect, the portion marked by 
frames 7-1 is similar to the portion marked by 
frames 3-5. 

In addition to the foregoing features character- 
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izing the variations in domain structure in an alter- 
nating field, careful examination of these varia- 
tions brings to light two facts. The first and more 
important is that,for the same field strength and 
the same direction of the field, the domain struc- 
ture has a different form, depending on whether 
these strengths are attained on increasing or de- 
creasing the field. This is quite clearly shown by 
a comparison of frames 2 and 8, 3 and 7. This in- 
deed constitutes the hysteresis of the variations in 
domain structure on the application of electric \fields. 

The second interesting fact is that the domain 
structure in the absence of a field (frame in the 
center of the figure) in its arrangement does not oc- 
cupy an intermediate position compared with the 
domain structure obtained when the field passes 
through zero value (frames 3 and 7). This fact is 
readily understandable if it is borne in mind thatin 
the initial state the domain structure in the field 
of view was unipolar. Evidently, under the action 
of the alternating field, in the first polarization re- 
versal cycles, this unipolarity diminishes some- 
what, there is some averaging out of the propor- 
tions of the tvo components, and the specimen as- 
sumes, in the field of view, so to speak,a new ini- 
tial state which, however, after removal of the fieldis 
lost fairly soon, and the domain structure again as- 
sumes the form shown on the frame in the center. 

Figure 2 shows the variation of the same por- 
tion of domain structure shown in Fig.1. In this 
case, filming was conducted at an electric field 
strength of 56 V/cm and frequency 1000 cps. The 
characteristics of the variation in domain structure 
that were pointedoutin the foregoing discussion of 
Fig.1, apply wholly in this case. It should, how- 
ever, be noted that at a frequency of 1000 cps, uni- 
polarity of the domain structure is less marked; in 
both limit positions (frames 1 and 5), the specimen 
is approximately equally remote from saturation. 
As our observations show, this circumstance is 
due precisely to the change in frequency, and not 
to the fact that the field strength at the frequency of 
1000 cps was rather higher, since, at the frequen- 
cy of 32 cps, unipolarity was quite definite for the 
same field strength. 

Observations at different field strengths but at 
the same frequencies show that with increase in 
field strength, the variations in domain structure 
occur not only by the growth of some domains at 
the expense of others, but also by the formation of 
fresh small domains in the limits of large domains 
previously existing (at smaller field strengths). 

The authors thank V.F. Parvov for assistance 
in the work. 
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The depolarization rate of the photoelectret state in sulfur single crystals is examined 
as a function of temperature; the thermophotoelectret state is discussed. 


We have previously [1] examined the charge on 
single crystals in relation to photopolarization con- 
ditions, the charge being measured by depolariza- 
tion by repeated illumination. Here we report de- 
polarization tests made at various temperatures. 

As previously [1], the specimens were plates 
of area 1 cm? and thickness about 1 mm cut perpen- 
dicular to the acute bisector of the optic axes from 
natural sulfur crystals. The depolarization current 
was measured with the vacuum-tube electrometer 
previously described [2]. A stabilized power supply 
was used in the polarization. 


Fig. 1. Apparatus for examining 
depolarization: 1) aluminum 
electrode; 2) specimen; 3) grid; 
4) springs; 5) quartz tubes; 6,7) 
thermocouple leads; 8) pumping 
tube. 
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The specimen was mounted in the evacuable 
quartz vessel shown in Fig. 1. The electrodes were 
the aluminum plate 1 and the brass grid 3, to which 
the polarizing voltage was applied. The specimen 2 
was gripped by the pressure of the springs 4 sup- 
ported by the quartz tubes 5. Electrode 1 was 
joined to the electrometer via a special lead (not 
shown). Grid 3 acted as a transparent electrode and 
was at the high potential. The temperature of the 
specimen was measured with a copper-constantan 
differential thermocouple attached to grid 3, with 
leads via tubes 5. Figure 1 shows light incident 
parallel to the polarizing field. 

The lower part of the apparatus of Fig.1 was 
fitted with a heater; it could also be inserted in a 
Dewar containing liquid nitrogen. The effects could 
be examined at a variety of temperatures, includ- 
ing continuously changing ones. 

The polarizing voltage was 300 V, which was 
applied for 10 min with the specimen exposed to the 
unfiltered light of a quartz—mercury lamp. The 
light was then cut off, the apparatus was placed in 
a Dewar containing liquid nitrogen, the required 
temperature was attained, the polarizing voltage 
was cut off, and the electrodes were short-cir- 
cuited. The depolarization was begun 5 min later, 
by exposing the specimen to the above light, whose 
intensity was as during the polarization. 

Figure 2 shows the depolarization of a single 
crystal; the initial falling part of the curve is for 
103°C and corresponds to only partial depolariza- 
tion. It is clear that this initial process was very 
largely complete within 10 min; at that time, the 
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Fig. 2. Discharge-current curve for a photoelectret 
heated and illuminated, with the corresponding tem - 
perature curve. 


heater was switched on, the specimen still being il- 
luminated. 

The heating produced a characteristic peak in 
the discharge current, which fell to zero shortly 
after the specimen had reached room temperature. 
The shape of the peak is governed by the heating 
rate; the curve of Fig. 2 is that appropriate to the 
initial temperature and rate of change shown. Simi- 
lar curves were obtained with other initial tempera- 
tures. 

These depolarization curves resemble those re- 
ported (e.g., in [2]) for thermoelectrets in the 
course of total-charge measurement by heating and 
integration of the discharge current; here the heat- 
ing was combined with continuous illumination, of 
course. 

This photoelectret state is thus also a thermo- 
electret one, because illumination at a low tempera- 
ture removes only part of the charge; total depola- 
rization requires heating to the temperature used in 
the polarization (here room temperature). Also, the 
state is not of purely thermoelectret type, because 
the heating must be accompanied by illumination; 
lack of the latter gives no measurable discharge 
current, although there is a slight loss of total 
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Fig. 3, Relation of In(Q* 10!) to 103/T; 
Q is charge released (coulombs) and T is 
absolute temperature. 


charge. The state may be called the thermophoto- 
electret one. 

This state arises by electron trapping at lo- 
calized levels and on account of a thermal activa- 
tion energy (gap between the excitation levels and 
the conduction band). The light excites the elec- 
trons to the latter levels, from which they can pass 
to the conduction band, but the transition probability 
is dependent on the temperature. This is the me- 
chanism that has been proposed [4] to explain photo- 
conduction in alkali halides; it gives an explanation 
of the depolarization rate as a function of tempera- 
ture and,hence,of the peak in the discharge during 
depolarization by simultaneous illumination and 
heating. 

The activation energy may be deduced as in [2] 
for specimens polarized as above at room tempera- 
ture. The depolarization curves were recorded at 
a variety of temperatures; the discharge current 
was integrated for each temperature, the charge be- 
ing proportional to the area under the falling branch 
of the curve in Fig.2. The logarithmic plot of Fig. 
3 is a straight line, and it gives the activation ener- 
gy as U= 0.09 eV (1.4 10: erg), which is very 
small, as is to be expected [3]. The result for the 
Carnauba wax thermoelectrets of [2] was rather 
larger (1.3 - 10°" erg). Polycrystalline sulfur 
gives photoelectrets of much lower activation ener- 
gy, which may be measured as described here. 

We are indebted for assistance to Yu. N. 
Martyshev and A.I. Delova, and we are indebted 
for discussions to Academician A.V. Shubnikov, 

G. Nadzhakov (Bulgarian Academy of Sciences) ,and 
I.S. Zheludev. Advice on the method from A. I. 
Froiman is also gratefully acknowledged. 
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The specific susceptibility and magnetization have been measured as functions of 
temperature up to 600°C for natural pyrrhotine; thermal hysteresis occurs in these 
properties below the A peak. The transition around the A peak is shown to be of 
isothermal diffusion type. Rapid cooling suppresses the A peak in the susceptibility 
curve; the magnetic parameters above the Curie point (about 300°C) are also very 
much affected by the presence of oxygen. This is discussed in relation to a second- 
ary ferromagnetic phase of Fe3;O, type with a Curie point of 570°C. 


Introduction 


Compounds of the transition metals (especially 
Fe and Cr) with sulfur, selenium, and tellurium 
are of considerable interest [1,6,8], because these 
belong to the class of substances that can become 
ferromagnetic under appropriate conditions. 

Natural pyrrhotine is hexagonal and has the 
composition FeS,; ;—FeS; 2, and is already one such 
substance; it has some very unusual magnetic prop- 
erties. 

There are two views on the cause of the ferro- 
magnetism. Néel's view [2] is that the lattice has 
vacant sites that are not filled by iron, whose dis- 
position is such that the mineral has two sublat- 
tices, whose spontaneous magnetizations are oppo- 
site in direction and unequal in magnitude. This 
gives rise to uncompensated antiferromagnetism, 
which is responsible for the strong magnetization. 
Néel considers that the vacant sites may become 
ordered at a certain temperature, which gives rise 
to even greater uncompensated antiferromagnetism. 
The second view is that taken by Japanese workers 
[3,4]. They neglect this vacant-site structure and 
consider that the uncompensated antiferromagne- 
tism arises from inequality of the two sublattices, 
one of which contains only Fe? ions, while the other 
contains Fe®’ and Fe", No final decision has been 
reached on the nature of the ferromagnetism; new 
experimental evidence is needed. 


We have measured the magnetic parameters as 
functions of temperature to obtain more detailed 
evidence on the behavior near the Curie point. 


Methods 


A specimen in the form of a disc, ellipsoid, or 
cylinder a few millimeters in size was suspended 
on a quartz fiber in a furnace with a bifilar wind- 
ing, which was placed between the poles of an elec- 
tromagnet. The susceptibility x was deduced from 
the force of attraction F via 

F 

X~ nH dH [de ’ 
where m is the mass of the specimen and dH/ dx is 
the field gradient. H(dH/dx) was measured with 
substances of known susceptibility. The force was 
measured with a VR-NV-20 direct-reading bal- 
ance. The field was measured with a small probe 
coil whose size was comparable to that of the 
specimen. The demagnetization factor for the spe- 
cimen was derived from tables. The specific mag- 
netization is given by 0 = XH. The measurements 
were made with fields of 1300 to 7500 Oe in air and 
with the specimen sealed in a bulb... The order of er- 
ror in the susceptibility was about 5%. 


Measurements Near the Curie Point 


The susceptibility decreases somewhat be- 
tween room temperature and 200°C, but it rises 
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Fig. 1, Susceptibility over the range 20-300°C in a field 
of 1340 Oe. Arrows show the sense of the temperature 
change, the numbers being the stages of heating and cool- 
ing. 


sharply at 220°C to give the A peak [5], which rep- 
resents very unusual behavior. The usual Weiss 
curve is followed between this peak and the Curie 
point (about 300°C). Temperature hysteresis occurs 
in the susceptibility curve below 200°C; the position 
and height of the A peak are dependent on the rates 
of heating and cooling [6]. 

Figure 1 shows the susceptibility for one spe- 
cimen at 1340 Oe at several rates of heating and 
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Fig. 2. Susceptibility as a function of time at a fixed 
temperature in a field of 1340 Oe. a) 221°C, halt 
during heating; b) 203°C, halt during cooling. 


cooling used in sequence. Curve 1 is the initial 
heating; cooling (curve 2) at the same rate gave a 
higher susceptibility below the A point. Curve 3 
lies below the previous cooling curve, which shows 
that the transition had not gone to completion dur- 
ing the cooling. Curve 4 of Fig. 1 shows that no A’ 
peak is produced when the cooling is very rapid, 
the entire curve from the Curie point down to room 
temperature being of the usual Weiss form. 

Néel's explanation for the A peak and the hys- 
teresis is that the vacancies become ordered at that 
point (superlattice I, in the terminology used by the 
Japanese workers [6]). Superlattice I is converted 
to another type of vacancy ordering below 220°C 
(superlattice II), which reduces the deviation from 
equivalence in the sublattices. The kinetic behavior 
suggests that the transition involves diffusion; Fig. 
2 provides strong evidence for this by reference to 
the susceptibility as a function of time at fixed tem- 
peratures of 221°C (this corresponds to the A peak 
after heating, curve a) and 208°C after cooling 
(curve b). The relaxation time for the transition 
from I to IJ is finite, so structure II is not per- 
fected during rapid cooling, structure I persisting 
to a certain extent. This explains why there is no 
X peak on rapid cooling (curve 4 of Fig. 1), be- 
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Fig. 3. Susceptibility over the range 20-600°C in a 
field of 1340 Oe. Arrows show the sense of the tem- 
perature change, the numbers being the stages of 
heating and cooling. The rates of change of tempera- 
ture are also shown, 
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Fig. 4. Susceptibility over the range 20-640°C in a field of 1340 Oe with- 
Out access to air. Arrows show the sense of the temperature change, the 
numbers being the stages of heating and cooling. 


cause structure Il scarcely forms at all, the ma- 
terial consisting almost entirely of the ferromag- 
netic phase I. Subsequent slower heating (curve 5 
of Fig. 1) allows diffusion and accelerates the con- 
version of structure I to structure II. 


Measurements Above the Curie Point 


The only paper on this region is that by Benoit 
[7], who found that the susceptibility fell rapidly 
near 570°C. It has been supposed [2,7] that struc- 
ture I is destroyed at this temperature, but our re- 
sults show that the effect at 570°C is of secondary 
origin. 

Measurements in air (Fig. 3) show a marked 
sensitivity to oxygen. The initial heating above 
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Fig. 5. Susceptibility of pyrrohotine at various field strengths 
over the range 300-600°C. The arrows indicate the direction 
of temperature change. 
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Fig. 6. Magnetization isotherms in the range 280-600°C. 
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300°C (broken line of Fig. 3) sometimes causes the 
susceptibility first to rise and then to fall at 570°C. 
Cooling curve 2 (after a certain time at 600°C) runs 
above the curve for the previous heating and has 

the typical shape of a curve for a two-phase state 
with two Curie points. This indicates that the min- 
eral takes up much oxygen at high temperatures 

and produces a second ferromagnetic phase whose 
Curie point is 570°C; this may be Fe3Q, (magnetite), 
whose Curie point is 570°C. 

A similar effect from a magnetite-type phase 
has been observed for FeTe, which falls in the FeS 
class of compounds [8]. 

We have recorded x(T) curves over the range 
20 to 600°C for specimens heated and cooled in 
tubes without access to air. Figure 4 shows that 
there is then no characteristic point at 570°C. Fig- 
ure 5 shows the range 300-600°C on a larger scale. 
The susceptibility shows a slight dependence on the 
field in this range, which almost vanishes at 570°C 
(the material becomes paramagnetic). Natural pyr- 
rhotine crystals evidently contain a small propor- 
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Fig. 7. Spontaneous magnetization in the range 
280-600°C, 


tion of the magnetite-type phase. Benoit's data give 
a curve lying above ours, which indicates a higher 
proportion of magnetite; the two sets of results 
give very similar susceptibilities above 570°C. 

The magnetite phase is indicated also by the 
magnetization isotherms from the Curie point to 
600°C (Fig. 6), which are straight lines above 2000 
Oe. Extrapolation to H = 0 gives the spontaneous 
magnetization 0, above the Curie point and up to 
570°C (Fig. 7). It is clear that og is due to atrace 
of a ferromagnetic phase of Fe3O0, type with 570°C 
as Curie point. 


Conclusions 


1. Thetransition near the A point in pyrrho- 
tine has been shown to be of isothermal diffusion 
type; no A peak occurs on very rapid heating or 
cooling, and the usual Weiss curve is obtained. 

2. Oxygen has a marked effect on the magnet- 
ic properties of pyrrhotine above the Curie point. 
The fall in susceptibility at 570°C is due to a phase 
of magnetite type formed on heating the mineral in 
air. 

We are indebted to Professor G. P. Barsanov 
for providing the specimens of natural pyrrhotine. 
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Structure and dielectric properties are reported for Ba(Ti, Fe)O3 with BaFeO,; contents 
from 0 to 5%; a discontinuity in properties occurs at about 3% BaFeO3;. A previous dis- 
cussion [1] of the properties of cubic barium ferrate (perovskite structure) is extended 
to the titanates and to isomorphous barium titanate—ferrate mixtures, whose proper- 


ties are reported here. 


1. Synthesis 


The cubic titanate and ferrate of barium were 
used as initial components for the Ba(Ti, Fe)O3. 
The titanate was made [3] by firing Ba(OH), -8H,O 
with TiO, in stoichiometric proportions for 20 h at 
400°C. 

The components were ground together in a ball 
mill for 5 h, and mixtures were pressed into discs, 
which were fired for 2 h at various temperatures in 
asilite oven. Table 1 gives the conditions used for 
some of the specimens. 


Dike Structure 


X-ray patterns were recorded with copper ra- 
diation. Specimens with up to 3% ferrate were used 
in an RKE camera, while those with 3-5% were al- 
so examined in a powder camera 86 mm in diam- 
eter. Figures 1 and 2 illustrate the relation of lat- 
tice parameters to ferrate content; Fig. 1 shows 
photometric curves for the group of lines having 
Zhi = 26 for various contents. The tetragonal ratio 
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here changes, and this is illustrated quantitatively 
by the c and a of Fig. 2, which omits the param- 
eters of specimen No. 4 (3% ferrate), because the 
resolution (Fig. 1) was too poor to allow reliable 
calculation. This photometer curve also has the ad- 
ditional weak peak denoted by H, which relates to 
the hexagonal form of BaTiO; [2]. This line be- 
comes stronger as the ferrate content is increased, 
while the lines corresponding to the perovskite 
structure become weaker. 


7) 


Fig. 1, Photometer curves for reflections hay - | 
ing rhy = 26 for BaTiO, containing BaFeO, 
(shown as mol.% on the curves). 
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Fig. 2. Relation of cell parameters 
to BaFeO, content for BaTiO. 


3. Dielectric Measurements 


The above change at 3% ferrate is revealed al- 
so by dielectric measurements. 

a. High Frequencies. Figures 3 and 
4 show the dielectric constant at 5 - 10° eps and E= 
100 V/cm as a function of temperature; the substi- 
tution method was used in conjunction with a hetero- 
dyne measurement of the resonance frequency. 
Figure 3 shows that the peaks in € corresponding 
to the first and second phase transitions (120°C and 
about 0°C, respectively) shift to lower tempera- 
tures as the ferrate content increases. Figure 4 
shows that ferrate contents of 3-5% reduce € (espe- 
cially at the peak), as well as the Curie tempera- 
ture. The last is linearly related to concentration 
(Fig. 5) up to 2% ferrate, but not above this. 


= 160 -50 0 50 


100 1,C 


Fig. 3. Relation of € to temperature at 5 - 10° cps and E = 
100 V/cm; curves 1, 2, and 3 correspond to specimens 1, 
2, and 8 of Table 1. 
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-150 100 7,°C 
Fig. 4. Relation of € to temperature at 5° 10° cps and E = 


100 V/cm; curves 4 and 5 correspond to specimens 4 and 5 
of Table 1. 


b. Low Frequencies. Figure 6 shows 
€ =f (T) as recorded with an MDP bridge at 50 cps 
with E = 550 V/cm; as in Figs. 3 and 4, there isa 
rise in € up to 2%, this being especially prominent 
around the high-temperature transition, but there 
is only a slight fall in the transition temperature. 
There is a sharp rise in € and a displacement of 
the weak high-temperature peak toward higher tem- 
peratures at 3-5%. The tan6=/ (T) curves at this 
frequency for 0-2% ferrate resemble the curve for 
pure BaTiO3, but the curves shift toward lower 
temperatures as the ferrate content increases, the 
minima in these curves occurring at the tempera- 
tures corresponding to maxima ine = f(T). The 
tan 6= f (T) curves for 3 and 5% ferrate have 
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Fig. 5. Relation of Curie point to BaFeO, 
content for BaTiOg. 
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Fig. 6. Relation of & to temperature at 50 cps and E = 
550 V/cm; curves 1-5 correspond to specimens 1-5 of 
Table 1. 


minima, but they differ greatly in shape from the 
others. 

Figure 7 shows € =f (E); € in all cases in- 
creases with E, but the gradient is inversely re- 
lated to the ferrate content. 


4. Discussion and Conclusions 


The structure and dielectric properties show 
that isomorphous replacement of Ti by Fe has 
marked effects. The material remains ferroelec- 
tric at low iron contents, but the phase transitions 
shift to lower temperatures. We have previously 
assumed [1] that ideal BaFeOs; with the perovskite 
structure would show a phase transition above the 
temperature found for BaTiO;; on this basis, and 
from the properties of other isomorphous series 


such as (Ba, Pb) TiO, [4], we would expect the Curie 


point to increase with the iron content, but this is 
not found. The reason appears to be that iron of 
valency four is unstable; the iron enters the lattice 
of BaTiO; in trivalent form, which reduces the 
Curie point. 
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Fig. 7. Relation of € to field strength at 50 cps; curves 1-5 cor- 
respond to specimens 1-5 of Table 1. 


The structure becomes increasingly irregular 
as the proportion of iron is raised; the collective 
interaction of the ions is reduced, and a hexagonal 
phase appears when the ferrate content attains 3%, 
this being detectable in specimens cooled slowly 
from 1380°C. The reduced stability range for the 
spontaneous polarization is easily explained if 
quadrivalent titanium is replaced by trivalent iron, 
as are the increase in the relaxation-type polariza- 
tion and the other effects. 
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These stresses are calculated for anisotropic elliptic cylinders and plates, and also 
for spheres of cubic crystals; constant rates of heating or cooling are assumed. De- 
tailed formulas are given for some particular cases. 
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These thermal stresses must be known inorder 
to choose the best conditions for growing or anneal- 
ing single crystals; they are also relevant to thede- 
sign of piezoelectric and pyroelectric devices for 
use under conditions of rapidly changing tempera- 
ture. It is general knowledge [1-4] that these 
stresses are relevant to all forms of heat treat- 
ment of materials, and the basic aspects of thermo- 
elasticity for isotropic bodies are considered in the 
papers cited. Anisotropic bodies have begun to be 
considered only recently. Shubnikov [5] has ex- 
amined the geometry of thermal expansion for crys- 
tals, and calculations have been reported [6] for 
the stresses in cylinders of cubic and hexagonal 
crystals for the cases of the cylinder axis along 
four-fold and six-fold crystallographic axes. In 
collaboration with others, we [7] have deduced the 
stresses in anisotropic plates with arbitrary linear 
anisotropy subject to the condition that the only 
temperature gradient is along the direction of the 
thickness. Results have been given [8] for the cy- 
lindrically orthotropic thick-walled cylinders. Here 
we give results for elliptic cylinders, plates with 
very general anisotropy, and spheres of cubic 


crystals. 
The first requirement is to find the tempera- 


ture distribution produced by the heating or cool- 
ing. We assume that the temperature at the boun- 
daries of the crystal changes at a constant rate h 
(deg/ sec), which is taken as positive for heating 
and negative for cooling. The temperature distribu- 
tion is then given by 
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eT aT 
dx" dx" aor 
subject to the boundary conditions Tp = Ty + ht, 
T(x, y,z;t) < ©, and to the additional condition 
T(x,y,z:0)= Ty. Here, x is the thermal-diffu- 
sivity tensor. The superscripts represent contra- 
variance with respect to both, i.e., that on trans- 
forming the coordinates from x! to x1' we have 


ik 


(k= 12,3) (1) 


iu ke 
ith! — OF : Ox” ik 


dx’ dak 


’ 


as distinct from covariant tensors, for which 


- dx? axl 

ik —- Oa. 

: dx” da" * 

Summation with respect to the various repeated in- 
dices is assumed. 


An approximate solution to (1) is 


T (x, y, 2; t) =T,+ht+ O(a, y, 2), (2) 
where © satisfies 


PD 


yik eee 
dx" Ox" 


=h (i, k= 2 3)\2 Op (3) 
The relative error arising from this approxi- 
mation does not exceed € provided that the time 


since the start is not less than 


R2 
Us = SS Ili. 
x € 


where R is the largest of the characteristic dimen- 
sions of the crystal, and ~ is the least of the prin- 


1See [3], for example, 
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cipal values of the thermal-diffusivity tensor. It is 
easy to solve (8) in all the cases of interest here; 
the general solution of the thermal-conduction equa- 
tion for anisotropic bodies has been discussed by 
Mitskevich [9]. 

Let the X and Y Cartesian axes lie along the 
principal axes of the ellipse in an elliptic cylinder 
defined by 


the Z axis lying along the axis of the cylinder. The 
XYZ system is used in conjunction with an £7 Zone, 
with = x/R,, n =y/R,. Then, 


h D } 
OS ey Ct ed (4) 


For an unbounded layer —R< z< R, heatedat 
the surfaces z = R and z = —R, we have 


h 


9472 


Oe (2? — ie) (5) 
For a sphere (radius R) cut from a cubic crys- 
tal, we have 


D=— (y+ 22 RY, (6) 


To find the stresses we have to solve a system 
of equations that in the general case consists of 
three Cauchy equations for elastic equilibrium and 
six equations for consistency. In the latter the 
components €j, of the deformation tensor must be 
expressed in terms of the generalized form of 
Hooke's law 


Si == STRIPS UMe Te is La Lek, abet == Hen) 


via the components 0 7m of the stress tensor, the 
components ©; of the thermal-expansion tensor, 

those of the elasticity SjxkjJm, and the temperature 
T. The boundary conditions are deduced from the 
absence of forces applied to the boundaries of the 

crystal. 


2. Stresses in an Anisotropic Ellip- 
tic Cylinder 


These stresses are completely determined by 
the temperature distribution of (4), so they also 
will be functions only of x and y or — and”. 
Lekhnitskii ({10], sections 17 and 21, pp. 87-97 
and 109-110) has used stress functions to solve this 
class of problem and has derived a system of equa- 
tions that these functions must satisfy; they were 
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derived on the assumption that there are no therm- 
al stresses, but they are readily extended to that 
case. 

The Z axis of the XYZ Cartesian system here 
lies along the axis of the cylinder, the X and Y 
axes lying along the principal axes of inertia of the 
cross section. We assume that mass forces are ab- 
sent and that no forces are applied to the bounda- 
ries of the crystal; the present solution may be 
superimposed on Lekhnitskii's if either class of 
force is present. The boundary conditions at the 
ends of the cylinder are satisfied only on average, 
zenerally speaking: the resultant force is zero, as 
are the torsional and bending moments of the forces 
applied to the end, the forces themselves being fin- 
ite. We foliow Lekhnitskii and introduce two stress 
functions F (x,y) and ¥ (x,y), with 


eee ee (7) 


The stresses given by (7) satisfy the equations 
of equilibrium identically; these stresses are re- 
lated to the deformations and the temperature via 
the generalized form of Hooke's law, and they must 
also satisfy the equations of compatibility and the 
conditions at the ends. For this to be so, it is 
necessary and sufficient that the stress functions 
are solutions of the following system of equations’: 


IPE BY pk Pe meets ad 


Tx 
Ly F + 2L,¥ = — Ly T+ (eye — ty) = 


Ty 
ye 


a) (Ce ae Vaz) 


Saat) 4 (8) 


subject to the boundary conditions® 


he. gale 


ie ae D 


OF 
On 


= (i) 


Here, the symbols are 
a a a 
L,= Nyzyz ron Nssayp Ondy + Nexzxz Oy?» 


08 08 
Ly = -— Nyyyz age tb (Myyxe + 2Nzxyyz) ax®ay 


03 08 
+3 (asyye ae 2Nxyxz ) Gc dy> + Nxxez Oye? 


’See[10], pp. 91-92, formulas (17.17) and (17.18). 


’These are applicable only to a singly coupled body; see [10],p.93, 
for extension to more general cases, 


Youlonn las 


a aa a 
LE, = Nyyyy Ca ee Any yxy Bay +2 (Mexuu ae ANxyxy ) Baay? 


o4 o* 
— ANxxxy CEO + Nexxxex0 aye (9) 
a) a) 
Ly=— ers a [x2 Gy 5 


0? a 
“teu Guay + Ts Gye’ 


e 
Ly = uy On? 


Nikim = Siklm —_~- 


(Clhy i SA, hy OE 


Vik = Ck 
ZZZZ 

The bar denotes averaging over a cross section 
perpendicular to the axis: 


7=-5\\/(, y) aedy. 


Also, # is a constant defined by the condition 


Y = 0, (10) 


which expresses the fact that there is no torque. 

The above system may be solved by noting that 
Tx = Ty = 0 in this particular case, and that njkJm 
and Yj, are covariant tensors in respect of coordi- 
nate transformations in planes z = const; then we 
put (7) and (8) in the 7 Z system as 


Ie Te io) 4 


oS8 = on = —_ —__ oS => oe 
me On Rt ob? R* d&dn ° 
ee! or prea has (ota 
= Re Oy? ie ee mt) 
LF + 2?L,¥ = — R*Ly®, 
L,F +2R? 1,8 = — R*L,O—R89, (12) 


in which the operators L' are derived from the cor- 
responding operators L by replacing x by é and y 
by 7, with R= vR,R,, and ® defined by (4). 

The solution to (12) is found in the form 


1 » : 
B= —— hRE ga (Va o* 77)", 


3 1 e ¢ 2 
WY = — + hRiq (bi — en) (1 —° — 7p). (13) 
Condition (10) gives ¥ = 0, while a, b, c, and qare 
given by 

A (vg + Yan) + (KQy = Ke Qe) vex + (KOs — Kn Qn) rns 
AD + 2KQ_Q, — Ke @ — KG | 


Dee KO: — Ke, O,\d = Freee Ketel 


bls & 


OSS [((KQn— Ke Qe) a — KY2 + Ke yez], 


1 
9 (65 + 7) 


q= (14) 
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Here, 


K = Merz, 
Ky = BM yenz + Nez2e2, 
IK = 3Nexe7 + Nnznz + 
Qe = 8ngce2 + 2Meqnz + Nonkee » 
Qn = 3Nynnz + 2enez + Neenz » 
A == K; K »— K?, 
D = 3nggeg + 3Mgnnn + 2NZenn + 4Nnenen - (15) 
The stress distribution is then 
3% = hga (1 — 2% — 37), 
am = hga (1 — 36? — 7). 
361 = 2hqazy, 
Che ee hq (c(1 — #2 — 34?) + 2bé7), 
om = > hag [b (1 — 382 — 4?) + ct], 


Q 
% 
t 


hq - x 1 j 
akark iat {ee (# sa =) + [a (Szezz + Snnzz) + CSe222 


+ OSnzz2] (4 ee 7?) —2 (ASnzz ++ OSyiz) &2 


ae (AS¢z > ia CS¢ zzz) ae +2 (2aSenzz + DSt22z + CSyz22) enh > (16) 


where a, b, c, andd are given by (14) and (15), 
o2Z being given by 


4 a: ( rl Co Ty aries) 
s " 2 y? 


52% = 


¥ H 3° nx s F i 
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(17) 


which is a generalization of the formula given in 
[10] (formula (17.9), p. 90]. 
in general, in what follows. 
The stresses near the boundaries may best be 
considered by introducing a local Cartesian system 
for each point on the boundary ellipse, the basis 
vectors coinciding with the unit vectors tangential 


(7) and normal (v) to the ellipse at that point. This 
system is defined by 


This o 24 is not given, 


t= f, écos p) — R, ysind, 
v= Rk, Esind+ Ry cos v, 
R& 
V Ree Rey” 
Ry7 


V REO + Rie 


sin) = 


cos) = 


(18) 
From (18) we find that on the cylindrical surface 
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zt = hg 


3Risin*6 + (R4+4R°R3 + R§) sin? 6 cos? + 3R}cost9 — R4sin?6— R$ cos?0 


Re sin? @ + RE cos? 8 : 


4 , ( b[3R3 cos? 6 + (R3 + 2R?) sin? 6 — R53] cos 6 


-5- hq ———— 
2 V R? sin? 6+ RS cos? 0 


¢ [3R*sin? 0 + (R? -+ 2R3) cos? © — R?] sin 0 


V R? sin? 8 + R3 cos? ie 
where 6 = arctan(n/é), whileo’”, o”T, and Ss5 = Seg = 2(Sq1 — S42), S14 = So5 = 0, Oy = Ay = Ag = 


oY’ are zero, as the boundary conditions require. 

It is convenient to put (16) in terms of a cy- 
lindrical coordinate system if the crystal is a cir- 
cular cylinder: 


rr = 3'T = hgg dy (R? —r”), 
¢=roe = 0, 


eo == 7? oY = figy dy (fh? — 372), 


1 : 
rZ = oT? == > hg (bo Sin g + Co cos ¢) (R? — 7), 


— 1 
oZ =ro?? = 


= = NG (b) cos © — cy sin ©) (R? — 37”). — (20) 


Here, @), bo, Co, and dy are given by (14) and (15) 


if € is replaced by x and 7 by y in the superscripts. 


The crystal may have a plane of elastic symmetry 
perpendicular to the axis of the cylinder, in which 
case the solution simplifies further if we introduce 
a coefficient A: 


Woe 1 ; Neos a Yuy (21) 
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and takes the form 


F=—+Ah(R?—7), 


; ae 0 (22) 


a 


rr = Ah(R?—r*), 99 = Ah(R? — 37%), 
(23) 


Table 1 gives A for cylinders cut from various 
crystals to give this plane of elastic symmetry, in 
which all the formulas are in terms of the principal 
systems of crystal—physics coordinates, * the four- 
and two-place superscripts being replaced in ac- 
cordance with the usual rules ([{12], pp. 389-390). 

As we would expect (see [2]), we have for an 
isotropic body that 


oo 4 af 


A= Te S11 + S12 saa 16x) (24) 


in which E is Young's modulus and v is Poisson's 
ratio. All the formulas in the table become of this 


form if we put Sg = S33 = Si1+ S13 = S23 = Sig» S44 = 


1 B= 1 a OZ a. 


Q, Wy = Xe = Wg =X, 


3. Stresses in Crystal Plates 


The temperature distribution in an anisotropic 
plate is defined by (5). The stress distribution cor- 
responding to a general temperature distribution 
along the thickness of the plate has been derived 
[7], and for the particular temperature distribution 
of (5) we have 


h 
6x22 


opz — ()) 


oik = 


Bik (R2— 32%), (i, k= 1, 2) 


(pene oe). 
(25) 


Here the Bik are given by the following formulas if 
the superscripts are replaced in accordance with 
the usual rules: 


nr 
or cil, 2, @ 
m A ’ (m,n 7 &) 
(26) 
where 
$11 Si2 Sig 
A= Si2 Soo Sag | ; 
Sig S26 Soe 


andA jy is the algebraic complement of the corre- 
sponding element in the determinantA. Detailed 
formulas have been tabulated [7] for the Bik of crys- 
tal plates of various orientations and various 
classes of elastic symmetry. 

One case of interest is a thin circular disc of 
radius R subject to a planar temperature distribu- 
tion in the squares of the coordinates: 


T=T,+k(a? + y?—R?). (27) 


The stress function may be put as follows by 
analogy with (22): 


4See [11], pp. 271-273, 326, and 327, or [12], pp. 18-21, 44-45, 


and 50-58 for the orientation of the axes of these systems, 
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F=—7 B(R—r), (28) 
where 
Be te ty 
88 sexe 3Syyyy + 8axyy + 4Seyxy (29) 


and the stress distribution becomes 
rr=B(R®—r%), o9 = B(R?—3r), 79=0. (30) 


4. Stresses in Spheres of Cubic 
Crystals 


Here the temperature distribution is defined by 
(6), and this produces the following stress distribu - 
tion: 


a — Ch(R®— x? — 2y? — 22”), oY? = Chyz, 
ow = Ch(R? — 22? — y? — 227), o** = Chazz, 
a7? = Ch(R? — 2x? — 2y? -- 27),  o*¥ = Chay, (31) 
or in polar coordinates, 
rr=Ch(@—r), $9 =Ch(R?—2r), ro=0. (82) 
Here 
me 20 é l 
ox 88) + 12sj2-F Sag (33) 


As we would expect [8], for an isotropic sphere: 


ak 
15x(4—v) ~ 


o 1 


Susy 


(34) 


We are indebted to V.R. Regel' for suggesting 
the topic and for advice on the work, and to B. N. 
Grechushnikov and V. L. Indenbom for valuable dis- 
cussions. 
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It has been found that there are specific directions of propagation of elastic waves in 
crystals belonging to different symmetry groups of the tensor of elastic constants. 
The specific directions are divided into several forms according to the character of 
the propagation of ultrasonic vibrations. It is shown that in the preparation of elas- 
tic shear waves in a crystal, the phenomena observed are in many respects similar 
to familiar optical phenomena. By means of quarter-wave crystal plates, it is pos- 
sible to excite a circularly polarized ultrasonic shear wave in an isotropic body. 


The classic theory of the elasticity of aniso- 
tropic media was developed in the 19th century in 
the works of Green, Christoffel, Kelvin, and other 
investigators [1] to explain optical phenomena in 
crystals on the basis of the mechanical theory of 
light. Direct confirmation of this theory could not 
then be made for elastic waves because of the in- 
adequacy of experimental technique. 

Currently, because of the development of meth- 
ods of ultrasonic research, it is possible not only 
to confirm the position of the classic theory of elas- 
ticity, but, in a number of cases, to develop and 
amplify it. This question is dealt with in many re- 
cent articles [2-5,8] by foreign authors. Until now, 
however, as far aS we are aware, there are no 
publications dealing directly with the analogy be- 
tween the propagation of elastic and electromagnet- 
ic vibrations in crystals except slight references to 
this analogy in a few articles [2,8,10]. 

The purpose of the present work was to investi- 
gate the propagation of ultrasonic vibrations in spe- 
cific directions in a crystal, and to draw an analogy 
between these phenomena and similar phenomena in 
optics. 


Specific Directions of the Propagation 
of Elastic Waves in Crystals 


It is sufficiently well known that three plane 
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elastic waves may be propagated in an arbitrary di- 
rection in a crystal: one quasi-longitudinal wave 
and two quasi-transverse waves. The vectors of 
the particles displaced in the waves are perpen- 
dicular to each other and, in the general case, some 
of them may coincide with the wave normal. 

Depending on the form of symmetry, many 
crystals may have specific directions for the prop- 
agation of plane elastic waves, directions in which 
purely transverse and longitudinal waves are prop- 
agated. Recently, Borgnis [2] has given a method 
of determining such specific directions for crystals 
of any symmetry. 

Following this method, we shall determine the 
specific directions for crystals of the rhombic sys- 
tem, and derive from the expressions obtained spe- 
cific directions also for crystals of symmetrical 
systems of higher symmetry (tetragonal, hexagonal, 
and cubic). 

We shall denote by n the unit vector, determin- 
ing the direction of propagation of the wave and de- 
termined by direction cosines 1;, ly, 13, while the 
vector of displacement of the particles in the wave 
connected with this direction, having the direction 
cosines mj, My, m3, will be denoted by r. As 
shown previously [2,3], for any n, the direction of 
the vector r and the velocity of propagation of the 
waves vj may be found from the system of equations 
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my, + mel’ 32 + mg Lg = my pv®, 
my Ix, + mg Loo + mz Pog = m2 pv", (1) 
m, I'y3 + mz D235 + mz I'33 = mz pv”, 


where Ij; = Phy for crystals of the rhombic system 
they have the form 


y= it Ci + i Cog + I C55? 

P= Lites “4 3 cos ls U5 cass 

P33 = Ui 55 + bb Cag 1 I C33; 

T'93 = Ly Ls (Cog + C44) 

Dy3 = 1, ls (Cis + 55), 

DP yo = Uy Ls (C12 + €g6)- (2) 


Here, p is the density of the crystaland the cjk are 
its elastic constants. The velocities of propagation of 
the three waves, connected with n, may be ob- 
tained from the determinant of the system of equa- 
tions (1) 


Li, — pv" Ih, I'h3 
L192 pdisa— pv”. Las = 0. (3) 
Ih3 Po, I'33 — pv? 


To each value of the velocity vj, obtained from 
(3), there correspond in equations (1) its permitted 
values of the direction cosines of the vector of dis- 
placement of the particles in the wave mj4j, Mg, M3j. 
For any one of the three values of the velocity, the 
quantities on the right-hand side of Eq. (1) may be 
regarded as components of a certain vector 
pitm4j Pv}; Mp j PV}; m3; vi} coinciding in direction 
with the vector r. The condition for the existence 
of specific directions of propagation of elastic 
waves in the crystal may then be written as follows: 


T,X HO o Py x D=sO, (4) 


Condition (4) means that one of the displacement 
vectors is parallel to the vector of the wave nor- 
mal. Actually, if one of the three displacement 
vectors is parallel to the wave normal, and the dis- 
placement vectors are at right angles to each other, 
all three waves are pure, while the direction of the 
normal found is the specific direction. Condition 

(4) may be written otherwise: 


lP2—Ilsp,=0. (4') 


Simple substitutions and replacement in (1) of 
mk by l,, in agreement with condition (4), for a 
crystal of the rhombic system, leads to Eq. (5), the 
solutions of which with respect to 11, 1), 13 give 
the directions sought: 


ly Pg —l3p2=0, 1,p3 — Isp, = 0, 


Il, (C (645 — 45 — 2C¢6) + I (— Coz + C13 + 2Cge) 
+ U3 (cig + 2C55 — Cos — 2C44)) = 0, 


ly ly (Ui (— err + C13 + 2C55) I3 (Cog + 2C44 — C12 — 2C66) 
+ 13 (¢s3 — Cis — 2C55)] = 9, 
I, Lg [U4 (Cra + 2Cgg — C13 — 255) + EB (Big Cxg — 20g) 
+ U3 (— C33 + Cog + 2C4a)] = 9. (5) 


Equations (5) are satisfied if 


bp le) ee be (6) 


Thus, the directions of all three crystallographic 
axes in crystals of the rhombic system are spe- 
cific propagation directions, independent of the 
relationships between the elastic constants of the 
crystal. 

Furthermore, equations (5) are satisfied 


C22 — Cog — 2Caa 


2 
3 

for = 6) — = 5 . 
I C33 — Cog — 4C4q 
2 . 

for i =0 Er C11 — C13 — 2055 

< 2 33 — Cg — 2655 ” 

2 

for ls = 0 5 Gi "Cis — 266 


#Z C22 — Cig — 2056 © (7) 


Specific directions in planes perpendicular to 
one of the crystallographic axes exist only for de- 
finite relationships between the elastic constants of 
the crystal, when the numerators and denominators 
in expressions (7) have the same signs. 

For the case 1; ~ 1, = 13 # 0, there may still 
be one more specific direction in a rhombic crys- 
tal, given by the direction cosines determined by 
the relationships 


2 
2 bd — kb —ad Ut df —kf+ ck 
2 ac—bce—af ’ 2 ~ kb bd + ad’ (8) 
where 
@ = C11 — Crp — 2Cyq, d = Cy, — Cig — 2Cpp, 
b = C33 — C13 — 255, f = Cgg — Cog — 2Cya, 


a 


= Co2 — Cig — 2WCog, K = Ca. — Cog — 2eyy. 


(9) 
It is easy to find from expressions (6)-(8) spe- 
cific directions also for crystals of the tetragonal, 
hexagonal, and cubic systems. 
For example, for the cubic system, putting 
C11 = Coz = C335 Cag = C55 = Cag, Cyg = Cy3 = Co3, We ob- 
tain the specific directions 
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L=1 (@=1, 2, 3) (6") 
l, =0, L=+h, l, =0, ly= +h, 1,=0, l=+1; 
(7") 
ES ee ee ee eee (8') 

Equations (9) reduce to the form 
Q=b=...=k=Cy— Cy. — 2Cyy, (9") 


i.e., they characterize the so-called anisotropic 
crystal. 

In the already cited article [2], the author, 
while finding the specific directions in crystals of 
the cubic system, failed to discover the directions 
expressed by formula (7') which will be discussed 
later. 

Before proceeding to the investigation of the 
specific directions, two remarks must be made. 

1. It is by no means necessary for all the 
specific directions determined from formulas (6)- 
(8) to exist in a crystal. Only formulas (6) are in- 
dependent of the relationships between the elastic 
constants of a crystal, and the directions defined 
by formulas (6) always exist in crystals of the sys- 
tems enumerated. With regard to the directions 
defined by formulas (7), they exist necessarily in 
crystals of the cubic system in planes perpendicu- 
lar to the principal axes, but they may not be found 
in the other systems. The directions defined by 


formulas (8) certainly exist in crystals of the cubic 
system, agree with the directions defined by for- 
mulas (7) in crystals of the hexagonal system, and 
are very rarely found in crystals of the tetragonal 
and rhombic systems. 

Thus, for example, for a Rochelle salt crystal, 
the directions defined by expression (8) do notexist, 
while of the three possible directions, defined by 
expression (7), only one exists in the XZ plane at 
an angle ~ = 15°20' to the X axis (for values of the 
elastic constants taken from [6]). 

2. When speaking of a specific direction, it 
must not be forgotten that it may occur several 
times in a crystal, corresponding to the symmetry 
of the elastic properties of the crystal. The com- 
plete set of specific directions in crystals of all 
systems is shown in Fig.1, where the light circles 
denote directions necessarily existing in the crys- 
tal, while the dots represent directions which are 
not to be found in some crystals of the given elas- 
ticity group. The specific directions are now shown 
in the bottom hemisphere. 


Investigation of the Propagation of 


Ultrasonic Vibrations along Specific 


Directions in a Crystal 


Figure 1 shows that the majority of specific 
directions of the propagation of elastic waves in a 
crystal are associated with the elements of sym- 


Fig. 1. Specific directions of propagation of elastic waves in crystals belonging to different 
symmetry groups of the tensor of elastic constants. 
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metry of the medium (coincide with the axes of 
symmetry or are situated in planes of symmetry). 
Knowing from equations (3) the velocities of prop- 
agation of the waves y; in a specific direction 

{1 1, l., 13}, if we find from equations (1) the di- 
rection cosines of the vectors of displacement in 
the waves, and carry out this investigation for all 
crystals, the specific directions may be divided in- 
to several groups. 

1. In any specific direction coinciding with 
the axes of symmetry 3, 4, 6 of the elasticity 
group, it is possible for the purely longitudinal 
wave to be propagated at one velocity and the two 
shear waves at another velocity. With regard to the 
directions of displacement in the shear waves, 
equations (1) fail to give a definite answer in this 
case. 

To confirm this statement, experiments were 
made on specimens cut normally to the axis of 
symmety of the elasticity group from NaCl and 
NH,H, PO, crystals, using an ultrasonic pulse ap- 
paratus [7] and two plates of Y-cut quartz, as trans- 
mitter and receiver. A short pulse of ultrasonic 
Shear vibrations of frequency 1.67 Mc was injected 
into the crystal. The quartz plate at the other end 
of the crystal received the pulse passing through 
the crystal; after amplification, this pulse was ap- 
plied to the plates of a cathode ray tube. A plate of 
Y-cut quartz transmits and receives ultrasonic vi- 
brations with a fixed direction of displacement 
along the X axis, just as the polarizers and ana- 
lyzers used in optics pass only light vibrations of 
a definite polarization. 

If an elastic shear wave is injected into aglass 
specimen and the relative orientation of the trans- 
mitting and receiving plates is varied, then, if the 
directions of displacement of the particles in both 
plates are parallel, the received signal will be a 
maximum. If the plates are crossed, the pulse on 
the screen of the instrument disappears. In all the 
intermediate positions, a law is obeyed, known in 
optics as the Malus law, 


A, = A) cos 9, 


where A, is the amplitude of the ultrasonic pulse 
with parallel position of the plates, and A, is the 
amplitude of the pulse when the directions in the 
plates are rotated in relation to each other through 
an angle ¢. 

In what follows, we shall call the transmitter 
and receiver the polarizer and analyzer, and en- 
deavor to use the more familiar terminology 
adopted in optics, emphasizing the close analogy 
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between the propagation of electromagnetic and 
elastic vibrations in crystals. 

The experiments made showed that with any 
arbitrary direction of polarization of the trans- 
mitted ultrasonic wave and with crossed positions 
of polarizer and analyzer, the ultrasonic pulse was 
not received by the analyzer, and the wave was 
propagated in the crystal with given direction of po- 
larization. 

This case is to some extent analogous to the 
case of the propagation of light along the optical 
axis of a uniaxial inactive crystal. A plate, cut 
perpendicular to the optical axis and placed be- 
tween crossed Nicols, remains dark on rotation of 
the microscope stage about the direction of propa- 
gation of the light. For this reason, we shall call 
the specific directions of propagation of elastic 
waves, coinciding with the 3,4,6 axes, axes of iso- 
tropy of the first kind. 

It should be noted that in crystals there may be 
more than one such axis. In cubic crystals there 
are three such axes, in tetragonal, hexagonal, and 
trigonal crystals there is one, and in crystals of 
the remaining classes there are no such direc- 
tions (see Fig.1). 

2. An exception to the case considered in 1 is 
provided by the [111] directions, coinciding with the 
axis of elasticity group 6 of crystals of the cubic 
system. As shown by de Klerk and Musgrave [8], in 
the propagation of an ultrasonic wave in the [111] 
direction, conditions are provided for the internal 
conical refraction of the elastic waves. Just as in 
the case of the axis of anisotropy of the first kind, 
two shear waves may be propagated at different 
velocities in this direction, and one longitudinal 
wave. The direction of displacement of the shear 
waves is also arbitrary. The specific nature of the 
[111] direction is the deviation of the ultrasonic ray 
from the normal to the wave front and the forma- 
tion by rays having different polarizations of a 
circular cone of rays whose aperture angle, A;,, de- 
pends on the anisotropy of the crystal and is de- 
fined by the expression 


tan; a2 | a ee 


V2! C1 — ere + Cag (10) 


It will be recalled that in the propagation of 
light along the optical axis of a biaxial crystal, the 
phenomenon of internal conical refraction is also 
observed; when the crystal is sufficiently thick and 
the beam of light sufficiently narrow, it consists in 
the formation of a ring of light (see Fig. 2a), the 
diameter of which depends on the thickness of the 
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Axis of 
of the second kind 


Fig. 2. a) Path of rays in internal conical re- 
fraction of light in biaxial crystals; b) path of 
rays in internal conical refraction of elastic 
waves in cubic crystals. vjy, VL2, Vi3) Velo- 
cities along the ray; v\y, VN2, YN3) Velocities 
along the normal corresponding to the two 
shear waves and one longitudinal wave. 


crystal and the relationship between its optical con- 
stants [9]. If the phenomenon of internal conical 
refraction is examined in monochromatic polarized 
light, a luminous point will be seen, moving along 
a circumference during rotation of the polarizer. 
In exactly the same way, in the case of the 
propagation of elastic waves, on rotation of the 
transmitting plate about the direction of propaga- 
tion, the beam passing through the crystal des- 
cribes a circle about this direction (Fig. 2b). The 
phenomenon of the internal conical refraction of 
elastic waves is characterized by large angles of 
refraction (for example, about 20° for KI) and also 
by the fact that the direction of propagation is here 
the axis of the refraction cone, and not its genera- 
trix, as in the optical case. We shall call such di- 


rections of propagation of elastic waves axes ofiso- 


tropy of the second kind, 

Axes of the second kind have so far been dis- 
covered only in crystals of the cubic system [8]. It 
is interesting to note that the cubic crystal has 
three axes of the first kind and four axes of the se- 
cond kind; speaking in the language of optics, it is 
a heptaxial crystal, perhaps the most interesting 
with regard to its elastic properties. 

3. Specific directions, coinciding with two- 
fold axes of symmetry, are widespread in crystals 
(Fig.1). Analysis of Eqs. (2) and (8) in this case 
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gives three different values of the velocity and three 
directions of displacement perpendicular to each 
other. For example, for a cubic crystal with prop- 
agation of an elastic wave in the [110] direction, the 
velocities of the waves v are equal, respectively, 
to 
M110) [(C11 + C12 + 2Ca4)/ 2p)", 

Biol. [(Cr — Cy2)/29]?, Moo) = [Cas / el", 
where the subscripts against the velocities denote 
the directions of displacement in each wave. 

Investigations of the propagation of ultrasonic 
pulses in the directions considered have been car- 
ried out on specimens cut from cubic (rocksalt, po- 
tassium iodide), tetragonal (ammonium dihydrogen 
phosphate), trigonal quartz), and rhombic (Roch- 
elle salt) crystals. 

As previously reported [10], analogy is ob- 
served in this case with optical phenomena occur- 
ring in the propagation of light perpendicular to the 
optical axis of a uniaxial inactive crystal. Ifa 
monochromatic beam of linearly polarized light is 
passed into a plate cut parallel to the optical axis 
of the crystal, and is regarded through an analyzer 
crossed with the polarizer, then on rotating the 
plate about the axis of the light beam, darkening of 
the plate will be seen every 90°, when the plane of 
vibrations of the incident light coincides with the 
plane of vibrations of one of the light waves prop- 
agated in the crystal. 

As shown in [10], in the example of the prop- 
agation of ultrasonic pulses along the X axis of am- 
monium dihydrogen phosphate, a pulse injected in- 
to the crystal with any direction of polarization in 
the YZ plane, is divided into two pulses in the 
crystal. If the dimensions of the crystal are ade- 
quate, these pulses may be observed separately on 
the oscillograph screen. The time of lag of one 
pulse in relation to the other is determined by the 
dimensions of the crystal and the magnitude of the 
double refraction, depending on the value of the 
elastic constants of the crystal. If the direction of 
displacement of the wave injected into the crystal 
is oriented along the bisectrix of the angle between 
the permitted directions of polarization of shear 
waves, then, if the analyzer is rotated about the 
direction of propagation, there will be seen on the 
screen either one pulse, when the analyzer andone 
of the waves have parallel directions of polariza- 
tion, or two pulses in all intermediate positions. In 
parallel and crossed positions of polarizer and 
analyzer, both pulses are approximately equal in 
amplitude, and the images on the screen of the 
cathode ray tube will be qualitatively identical [10]. 
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Fig. 3. Specific direction of propagation 
in the YZ plane in a beryllium crystal. 
VL1» VL2: VL3) Velocities along the ray; 
VN1> YN2> YNg) Velocities along the nor- 
mal, of the two shear waves and one longi- 
tudinal wave, respectively. 


In what follows, we shall call the directions co- 
inciding with the axes of isotropy of the second kind 
in the crystal the directions of maximum double re- 
fraction. There are three such directions in a rhomb- 
ic crystal, four in a tetragonal crystal, six ina 
cubic crystal, and an infinite number in hexagonal 
crystals (any direction perpendicular to the 6 axis). 

4. As an example of the specific directions 
situated in planes of symmetry, we shall examine 
the direction which may exist in crystals of the 
hexagonal system; it is defined by the formula 


€33 — C13 — 2Caq 
Soak ie srasemmraneons 
11 — 613 — 4€q4 


(11) 


where ¢ is the angle reckoned from the 6 axis. In 
Fig. 1, the aggregate of such directions is shown 
by a dot-dash line. 

For the example, we shall consider crystals of 
zine and beryllium and use the nomenclature of 
their elastic constants given in [4]. Equation (11) 
shows that in a plane containing the principalangle, 
zinc has no specific directions differing from ¢ = 
0° and y = 90°, which have already been considered. 
A beryllium crystal in this plane has a new specific 
direction defined by the cosines [0, 0.594, 0.805], 
Mss 5 (O= Blais) a 

Calculations give for this direction [0,15,13] 
the following values for the speeds (vNj) and direc- 
tions of displacement of particles in the waves, 
5.19 - 10° cm/sec and 5.80 - 10° em/ sec for shear 
waves; the directions of displacement give [0, 1,13]; 
[1,0,0], and [0, -l3,l.], respectively. 


K. S. ALEKSANDROV 


All three rays are situated in a plane contain- 
ing the direction [0,/2,/3] and the 6 axis, and only 
the longitudinal wave deviates from the normal. The 
transmission of energy in the shear waves occurs 
along rays deviating from the normal by angles A,= 
4°20' and A, = 9°50' (see Fig. 3). As in the case of 
propagation along the axis of isotropy of the second 
kind, the velocities of sound along the ray for the 
two shear waves differ from the velocities along 
the normal. Here, however, the angles of devia- 
tions of the rays from the normal are different for 
waves having different directions of polarization; 
there is double refraction, and the shear wave with 
arbitrary direction of displacement is resolved in- 
to two waves, emerging at different points on the 
surface of the crystal. In other words, the shear 
waves are propagated in this direction in the crys- 
tal in the same way as light incident perpendicular 
to an arbitrarily cut plate of a biaxial crystal. In 
both cases, both waves deviate from the normal by 
different angles and have different normal and ray 
velocities. 

5. With regard to the specific directions 
which do not coincide with the elements of sym- 
metry of the medium and have not been especially 
investigated by us, it is possible to say merely 
that, according to the definition itself of specific 
directions (4), the longitudinal wave will be prop- 
agated along the normal. The points of emergence 
of rays of shear waves will be situated arbitrarily 
in space and the waves themselves will have differ- 
ent ray and normal velocities. In optics, this case 
is analogous to the foregoing. 


Fig. 4. Quarter-wave orienta- 
tion of plate cut from quartz 
perpendicular to the X axis. 
1-1 and 2-2) Directions of dis- 
placement of particles in waves 
propagated at velocities v, and 
V- 
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Circularly Polarized Elastic Shear 
Waves 


Just as in optics, circularly polarized elastic 
waves may be excited in an isotropic body by pass- 
ing a linearly polarized ultrasonic wave through a 
linearly double-refracting plate, displacing the 
phase difference of the two waves propagated in it 
by 7/2. 

It is not a difficult matter to make a quarter- 
wave plate for elastic waves. As will be seenfrom 
the foregoing, the most suitable are the directions 
which we have termed the directions of maximum 
double refraction. We shall calculate such a plate 
of quartz, selected for this purpose from considera- 
tions of mechanical strength. In the direction of the 
X axis in quartz (Fig.1), two purely shear waves 
are propagated at velocities vy = 3.33 ° 10° cm/sec 
and v, = 5.09 - 10° cm/sec. Figure 4 shows the 
polarization directions of these waves. In thesame 
figure, the dashed line shows the method of cutting 
the plate. We determine the thickness of the plate, 
taking into consideration that the available appa- 
ratus [7] gives ultrasonic pulses of a frequency of 
1.67 Mc. The phase difference G between the waves 
expressed in wavelengths is [11]: 


Gd —s>), (12) 


where d is the plate thickness, A, and A, are the 
wavelengths, whence, for our case, d= 1.445 mm. 

To reveal a circularly polarized elastic wave, 
the quarter-wave plate was stuck to a glass block 
with the dimensions 25 < 25 xX 100 mm, and pulses 
of shear waves polarized in the direction of one of 
the edges of the plate (along the diagonal between 
the permitted polarization directions) were injected 
into the block through the plate. The pulse was re- 
ceived by the analyzer at the other end of the block 
and had approximately the same amplitude for any 
position of the analyzer. If a second-quarter plate 
was inserted between analyzer and block in crossed 
position with the former, compensation of the phase 
difference occurred, and on rotation of the analyzer 
the pulse reached a maximum value twice per re- 
volution and disappeared twice. 

With a quartz plate of double the thickness (d= 
2.89 mm), the plane of polarization of the elastic 
wave could be rotated through 90°, which in optics 
is a characteristic feature of a half-wave plate. 


Summary 


On the basis of the method used by Borgnis [2], 
relationships have been derived for defining the 
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specific directions of propagation of elastic waves 
in crystals of the rhombic system. From these 
equations, by means of various simplifications, 
corresponding expressions may be obtained for 
crystals of the tetragonal, hexagonal, and cubic 
systems. 

It is shown that in crystals of the cubic sys- 
tem there are specific directions defined by for- 
mula (7') and not mentioned in [2]. A stereograph- 
ic projection (Fig.1) is given showing all the spe- 
cific directions which may exist in the different 
groups of elastic symmetry (symmetry of the ten- 
sor of elastic constants). 

An experimental investigation has been made 
of the propagation of ultrasonic vibrations along 
specific directions in various crystals. In allcases, 
an analogy has been found between the propagation 
of ultrasonic and light vibrations in crystals. On 
the basis of this analogy, all the specific directions 
in crystals have been divided into several groups, 
depending on the conditions of propagation of elas- 
tic waves in each direction. 

A direct continuation of the analogy made is 
the production of quarter-wave plates for elastic 
waves, and the experimental production of circu- 
larly polarized elastic waves. 

The aforementioned analogy is true only of 
shear waves and only in specific directions in the 
crystal. We have not considered the longitudinal 
wave, Since, in any specific direction in the crys- 
tal, the longitudinal wave is propagated along anor- 
mal. In the experiment,transmitters of Y-cutquartz 
were used, giving a pure shear wave with sufficient 
approximation for qualitative experiments. 

The case of the propagation of elastic waves in 
a crystal is to some extent a generalization with re- 
spect to the propagation of electremagnetic vibra- 
tions in crystals. 

If later a general theory of the propagation and 
refraction of elastic waves in crystals were to be 
formulated and the Fresnel surfaces and Huygens 
structures, which have become classic for optics, 
were to be determined for elastic waves, the sur- 
faces used in optics would form a special case of 
these surfaces. A considerable advance in this di- 
rection has been made in recent articles by Mus- 
grave [3,4,5,6], although they are not free from 
individual errors and shortcomings. 

The author is sincerely grateful to Academi- 
cian A.V. Shubnikov for his general supervision of 
this work. 
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An apparatus is described for forming spheres a few tenths of a millimeter in diam- 


eter. 


It is often difficult to calculate structure fac- 
tors from measured intensities because it is very 
complicated to compute the absorption as a function 
of direction for a single crystal of irregular shape. 
This difficulty is very largely overcome if the crys- 
tal is cut to a shape such that the absorption is 
easily computed for all directions. The most con- 
venient form is a sphere, and the absorption cor- 
rections for spherical crystals have been tabulated 
[1], so the calculations become quite simple. 

Here we describe a device that enables one to 
cut a sphere of diameter 0.2 mm from a single 
crystal of any shape. 


ie ranciples sand Desien 


Several devices have been described [2-5] for 
this purpose, the basic principle being as follows. 
The crystal is enclosed in a chamber coated with 
abrasive material and is caused to rotate freely; 
the crystal tends to rotate around its axis of larg- 
est moment of inertia, so the parts that strike the 
walls are the ends of the largest linear dimension. 
The result is that the crystal gradually acquires a 
shape having no largest dimension, i.e., a sphere. 
The crystal may be driven by centrifugal forces or 
by a jet of compressed air; the abrasive may be 
emery paper or grindstone material. The following 
conclusions may be drawn on methods previously 
described. 

A. Compressed air is much preferable tocen- 
trifugal devices, because the crystal is then kept 
continuously on the wall, and so it becomes rounded 
much more rapidly. Moreover, compressed air is 
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Experiments with various crystals are reported. 


much less liable to damage the crystal, which is 
especially important for materials that cleave eas- 
ily. 

B. Previous devices experience some diffi- 
culty in reducing the diameter to 0.1-0.2 mm;many 
of them have gaps (e.g., between the fixed and mov- 
ing parts, or between strips of emery paper) in 
which the small spheres become wedged. Also, 
emery paper becomes worn quite rapidly, and re- 
placement represents an unnecessary difficulty. 

We have designed a new instrument free from 
these deficiencies, which also has a special port 
through which the spheres are easily extracted 
without damage. The crystal is rounded by a car- 
borundum grinding wheel of grain size 400 having 
an internal hole of V shape (Figs. 1 and 2). The 
compressed air enters this space via a small tan- 
gential channel. The large central holes on the two 
sides are sealed with close-fitting plates, one of 
which is fixed and has in its center a Lucite window. 
The other can be moved and has two ports, both 
closed by fine wire grids, one of which is flat and 
allows the air to escape. The other is convex and 
is used to insert and remove the crystals. 

The adjustable plate is turnéd to expose the 
large grid; the crystal is placed in it, and the de- 
vice is turned to allow the crystal to fall into the 
interior space. This hole is closed by a fine wire 
grid during the grinding. The crystal is subsequent- 
ly extracted by opening the large grid again and 
tilting the device to cause the crystal to fall into it, 
perhaps with the aid of a slight puff of air. The 
crystal is then readily removed; forceps should not 
be used, especially for soft materials. 
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Fig. 1. Schematic of the device, 
1) Grinding wheel; 2,3) closure 
plates; 4) air inlet; 5) holder for 
wheel. 


Fig. 2. Some components. 1) Carborundum wheel; 2) metal plate, 1.5 


mm thick with (a) flat and (b) convex grids; 3) metal plate, 1.5 mm 
thick; 4) inlet; 5) wheel holder. 


ees tSsTon= Wate milals 


This device has been in use in our institute for 
about a year and has been used on various materi- 
als, as follows. 

1. Careful operation has enabled us to pro- 
duce spheres a few tenths of a millimeter in diam- 


Sphalerite 


SE Cae een 


Rutile — 
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Fig. 3. Size as a function of time for various materials at a 
fixed air flow rate. The slopes of the lines cannot be com- 
pared directly, because the fixed rate was not the same for 

all materials; for instance, a high rate was used with the NaCl. 


TABLE 1 


Material 


Sphalerite (ZnS) 


KyS,Os 1.00 
NaCl, KBr 0.04 
NayS,0, * 5H,0 0.30 
Rutile (TiO,) 0.01 


O15 


010 


005 


Diameter, mm 


20 25 IO 55 
Flow rate 
Fig. 4. Loss in diameter in 5 min for a sphalerite crystal as 
a function of flow rate (cm*/sec). 


eter from nearly all materials, the differences in 
diameter between directions not exceeding 1%. Dif- 
ficulties were encountered only in the following 
cases. 

A. Crystals that cleave readily are liable to 
split off small plates, so these are ground only to 
roughly spherical form (edges and corners re- 
moved). Examples are NaCl and KBr. Even these 
give satisfactory results if the device is operated 
carefully at a low speed with a fine wheel. 

B. If the hardness varies a great deal with di- 
rection. The product may then be an ellipsoid 
rather than a sphere, the major axis corresponding 
to the direction of greatest hardness. This is the 
case for rutile, for example. 

C. Crystals of very low melting point or show- 
ing transitions must be cooled suitably during 
grinding, on account of the heat produced. An ex- 
ample is K,S,O3 * 2H,O, which loses one molecule 


Fig. 5. Glass holder for stor- 
ing crystal spheres, 
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of the water of crystallization between 17 and 30°C, 
and the other at 35°C. 

2. The diameter has a linear relation to grind- 
ing time for a fixed air flow rate, the slope of the 
line being governed by the hardness (Fig. 3). Table 1 
gives the reduction Ad in the diameter produced in 
5 min with an air flow rate of 29 cm*/sec. 

3. Higher flow rates accelerate the reduction 
very greatly (roughly as the square of the rate, 
Fig. 4). Very high rates should not be used, on ac- 
count of excessive danger of damage or even disin- 
tegration. 

4. The small spheres are best stored in glass 
plates with small depressions ground in them, 
which should be fitted with suitable covers (Fig. 5). 
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The Komsomol'skaya Pravda Plant in Moscow 
has produced a plastic kit called "Mosaic" for con- 
structing a variety of mosaics; it is intended main- 
ly for stimulating conceptions of symmetry in 
children, but it is also of great value in instruction 
in crystallography, demonstration of atomic struc- 
tures, and mineralogy. The hexagonal parts come 
in 5 colors, each fitting via a rod into a panel with 
a hexagonal net of holes. This gives a simple and 
cheap way (only 17 roubles, or less than the sim- 
plest crystal structure model) of constructing vari- 
ous patterns based on anion close packing with ca- 
tions in the holes (many compounds and minerals 
consist essentially of patterns of this kind). Each 
hexagon may be considered as the projection along 
a triad axis of a Pauling octahedron, while a pat- 
tern built up from octahedra of different colors 
(i.e., cations) is (see [1] and [2]) present in many 
structures; sometimes it forms the entire structure. 

"Mosaic" enables one to construct (or to have 
a student construct, in an examination) the follow- 
ing patterns given in [1,2]: corundum—mica, il- 
menite, carbonate, spinel; rutile—cassiterite, dia- 
spore —goethite, anatase —brookite —columbite; and 
numerous patterns for minerals with orthorhombic— 
monoclinic sequences: norbergite, chondrodite, 
humite, clinohumite, and olivine. The monticellite 


pattern may also be demonstrated in the latter. The 
pyroxene —ramsayite and amphibole patterns of [1] 
are also well shown, and careful inspection will re- 
veal the "secret" pattern of chrysotile —asbestos 
(see [2], section V, p. 18), although there is noth- 
ing particularly secret about it. 

The models made with "Mosaic" would be even 
better if the parts were not very flat hexagonal 
pyramids, but consisted of a triangle in a hexagon. 
This would give a complete imitation of an octa- 
hedron and would enable one to pass readily to ar- 
rays consisting of tetrahedra. It would probably be 
difficult to get the molds altered for this purpose, 
although it would not affect the beauty of the hexa- 
gons at all. 

This kit will enable one to economize on paper 
printed with a hexagonal millimeter grid (in very 
short supply), which so far has been the basic me- 
dium for demonstrating and interpreting crystal 
structures on the basis of close packing. 
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Rochelle-salt crystals irradiated by the rays 
from Co*’ show altered electrical polarization; the 
usual hysteresis loop (Fig.1la) at a certain dose be- 
comes a double one (Fig.1b) typical of an antiferro- 
electric. Higher doses eliminate the loop; the po- 
larization becomes linear (Fig.1c). 

The following aspects are of interest here. 

1. Anomalous polarization as a function of 
temperature for weak and strong fields. 


-20 -10 0 10 20 50°C 


Fig. 2, Dielectric constant of Rochelle salt as a func- 
tion of temperature: a) unirradiated; b) irradiated for 
2h; c) irradiated for 11 h10 min. Dose rate, 0.1 ° 
1019 eV/h-g; measurements made at 8 V/cm with a 
capacity bridge, 


2. Changes in the domain structure and piezo- 
electric properties caused by irradiation. 

The measurements were made always on the 
same specimens at increasing total dose levels, the 
measurements being performed in the intervals be- 
a, tween successive doses. 

Curve a of Fig. 2 shows the dielectric constant 
as a function of temperature for one of the speci- 
mens, which had the usual two peaks at the Curie 

Fig. 1. Hysteresis loops of points before irradiation. These peaks become 
Rochelle -satts a) unirradiated; lower and move together (curve b of Fig. 2) as the 
oes es an ae dose increases, the final result being a single 
perature +18°C; frequency 50 broad hump (curve c). Curve b of Fig. 2 corre- 
cps; Emax = 1100 V/cm; dose sponds to the double hysteresis loop of Fig. 1b; 
rate 10° ev/h-g. curve c to Fig. lc. 


Cc 
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Fig. 3. Part of the double hysteresis loop 
given by Rochelle salt, 


The hysteresis loop was recorded at the same 
time as the dielectric constant was measured with 
a bridge in weak fields. The dielectric constants 
deduced from the angles ©; and ©, of Figs. 1b,c 
agree closely with the values for curves b and c of 
Ihe re 

Angle 8 on the saturation branch is appreciably 
less than angle & on the linear part (Fig.3), so the 
mechanisms of polarization in the two ranges 
(saturation and linear) are not the same. 

The domains in a specimen giving a double loop 
began to change in orientation only near the critical 
field strength (the field at which the nonlinear rela- 
tion of polarization to field appears). 

It seems clear that the linear part represents 
an extension of the initial part of the polarization 
curve for the ferroelectric (the region of initial di- 
electric constant). 

This elongation of the initial part may be due 
to radiolysis products, which are trapped in the 
lattice and produce stresses (tend to lock the do- 
mains in position). This means that much higher 
field strengths are now needed in order to reorient 
the domains. Prolonged irradiation suppresses 
orientation change altogether. The crystal is thus 


in a different state, in which the antiparallel dis- 
position of the domains gives rise to linear polari- 
zation for fields below the critical value; the simi- 
lar linear region for an antiferroelectric arises 
from the antiparallel array of the dipoles in adja- 
cent unit cells [3]. 

In some cases the domain structure became of 
smaller scale, whereas in others it was almost un- 
altered. 

The piezoelectric constant gi, was measured on 
a 45° X-cut specimen in the intervals between suc- 
cessive doses; there was no appreciable change in 
this quantity. 

This anomalous polarization resembles the po- 
larization of an antiferroelectric and points to a 
close connection between the properties of ferro- 
electrics and antiferroelectrics, although there are 
marked differences also, of course. Rochelle salt 
is of special value here, because the type of polari- 
zation is readily altered by adjustment of the radia- 
tion dose, the two states being ferroelectric and 
similar to antiferroelectric. 

We are indebted to I.S. Zheludev for direction 
of this work. 
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The range 1-30,000 kg/cm? was covered in 
these measurements on indium of 99.995% purity. 
The hydrostatic pressure was produced in a pres- 
sure amplifier designed at this institute [1]. The 
fluid was a mixture of isopentane with n-pentane, 
the pressure being measured to +100 kg/ cm? with 
a manganin manometer. The temperature was 
measured to 1.5°C with an iron—nichrome thermo- 
couple, the melting point being determined by dif- 
ferential temperature analysis with photographic re- 
cording. The methods have previously been de- 
scribed in detail [2]. 


250 


200 


90 0 20 30 
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Figure 1 shows the results; the melting point 
rises smoothly from 156°C at atmospheric pressure 
to 280°C at 30,000 kg/cm’, the mean gradient be- 
ing 4.13 - 107° deg-cm?/ kg in the range 1-30,000 
kg/em?. The curve is somewhat concave to the 
pressure axis, for the slope near the origin is 
4.33 - 107° deg-cm?/ kg. Oelsen et al. [3] give the 
latent heat of fusion at atmospheric pressure as 
781 + 4 cal/g-atom, which with Clapeyron's equa- 
tion gives the change in specific volume on melting 
(from the initial gradient) as 2.93 - 107? em?/g or 
0.336 cm®/g-atom. 

No polymorphic transitions are observed 
throughout this range of temperatures and pres- 
sures, 
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Ultrasound affects crystallization because it 
gives rise to numerous fresh nuclei in the super- 
cooled material if a seed is already present. The 
precise number of nuclei may vary, which is re- 
flected in the final polycrystalline structure. The 
effects of ultrasound on the latter have been ex- 
amined many times, but the effects on single crys- 
tals have not been recorded. A growing crystal 
draws the molecules or groups it needs from the 
adjacent medium, so the major factors here are 
diffusion and concentration currents (for solutions) 
or self-diffusion and thermal conduction (for melts). 
It is debatable whether ultrasound influences ex- 
change of molecules and diffusion, although these 
processes are basic to many of the effects of ultra- 
sound. 

We have examined the effects on the linear 
growth rates of the octahedron faces on single crys- 
tals of potash alum; the rates in the field were com- 
pared with those outside the field, the conditions 
otherwise being the same. The ultrasound was pro- 
vided by a quartz plate driven at 2 Mc with a maxi- 
mum intensity of 0.2 W/cm’. 

The beam passed through the bottom of the 
crystallizer, which was placed at the surface of a 
bath of transformer oil. The crystallizer was a 
double-walled glass vessel, the space between the 
walls containing circulating water at a constant 
temperature. The two crystallizers were con- 
nected in series in the water circuit. The growth 
rates were measured as the rates of displacement 
of the growth fronts as indicated by traveling 
microscopes reading to 0.01 mm. The temperature 
of the solution remained constant to 0.1°C; the 
saturation temperature was measured from concen- 
tration currents as viewed by schlieren methods. 


The solution was run into the preheated crys- 
tallizers at a temperature of 7-8°C above its satu- 
ration point. The vessels were then quickly covered 
with lids having holes to admit the crystal holder 
and thermometer. The solutions were stirred tho- 
roughly before the weighed seeds were inserted. 
Then the ultrasonic beam was directed into one 
vessel, while corresponding faces on both crystals 
were observed via the microscopes. The linear po- 
sitions were recorded, initially every 15-30 min 
and then every hour (total duration of run 13-15 h). 
At the end the crystals were removed, dried, and 
weighed; the change in weight gave a check on the 
growth rate. Dissolution can be examined in the 
same way. 

The ultrasound was found to accelerate disso- 
lution in unsaturated solutions, the mean slope of 
the curve relating linear displacement to time be- 


Face displacement 
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ing about doubled. Again, the growth rate in a 
saturated solution was raised by the ultrasound; 
here the mean slope was increased by a factor of 
2-3. The gain in weight in a period of 10-15 h was 
50-60 mg; Fig. 1 illustrates the general behavior 
of the linear displacement. The ultrasound accel- 
erated the growth of the face turned toward the 


source and also (to some extent) that of the other 
faces, because the relative increase in mass was 
always greater than for the control crystal. 

The ultrasound had this pronounced effect on 
the growth only in weakly supersaturated solutions 
and at low intensities. High intensities merely pro- 
duced the normal disintegration of the crystal. 
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Introduction. §1. X-ray detectors. §2. Main parts of a diffractometer. 83. Inten- 
sity measurement, quantitative phase analysis. $4. Instrumental factors in shapes 

of diffraction peaks, crystallite-size determination. 85. Precision cell-parameter 
measurement. §6. Intensity measurements on single crystals. 87. Texture studies. 
$8. Other applications. 


Introduction Geiger counters are most commonly used in 
diffractometers [7-12] because the above advan- 
tages are here accompanied by good stability even 
in simple equipment. Gas proportional counters 
[12-15] and scintillation counters [15-21] have re- 
cently come into use for other reasons. 


Diffractometers are being increasingly used in 
x-ray structure analysis. These are instruments 
for recording the pattern and measuring the intensi- 
ties with counters. There are types for use with 
powders, types for single-crystal studies, others 
for texture studies, and ones for special purposes. 
These instruments reduce the time needed while in- 
creasing the accuracy and sensitivity, which has 75 
meant more general use of x-ray structure analysis 100 (shy it 
in research and (especially) industry. 

The URS-50I x-ray equipment is in regular pro- 
duction in this country; it has ionization recording 
and is meant mainly for use with powders and in 
orienting piezoelectric crystals [1,2]. The general 
introduction of diffractometers in research and 
routine laboratories requires early exploitation of 
foreign work on devices of this type, together with 
development of methods adapted to the URS-50I and 
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1, X-Ray Detectors 


Fig. 1, Efficiency curves for counters; 1) scintilla- 


The detectors other than films are ionization tion counter with Nal(T1) crystal of thickness 0.7 
chambers, Geiger counters, gas proportional coun- cm, with 127 Be window; 2) the same, 0.1 cm 
ters, and scintillation counters. All of these reduce Nal(TD, 127 p Be; 3) 0.05 cm Nal(TD), 127 p Be; 


4) proportional counter filled with Xe at 300 mm 


the time required and give the intensities directly. Hg, active length 2,7 er re 
’ .7 cm, pt Be window; 5) the 


Counters are more sensitive and accurate than 

fil Puesie at h b 3-61: th f same, 300 mm Hg of Xe, 2.7 cm, window 13 u of 
1ims an ee ion chambers [3-6]; the error o mica +127 p Be; 6) the same, 500 mm Hg of Kr, 
measurement with a counter can be made as low as 2.7 cm, 13 y mica +127 p Be; 7) Geiger counter 


1%. 550 mm Hg of Ar, 10.0 cm, 13 p mica. 
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TABLE 1. Comparison of Types of Counter 


581 


Ratio of 102 line ' 


Efficiency, Dead time, Discr, 
Type yo of graphite to x- 
sec . loss, %o 
ray background 
Geiger 200 2.3 
Proportional i 4.3 
Scintillation 1 4.5 


*Filled with argon, 
**Filled with krypton. 


The most important characteristics of a coun- 
ter are as follows. 

Efficiency. This is defined as the ratio of 
the number of counts recorded to the number of 
quanta passing through the counter slit; it is direct- 
ly related to the sensitivity and inversely related to 
the time needed to attain a given level of error. The 
efficiency of a Geiger counter [10,11,22] or pro- 
portional counter [22] is dependent on the absorp- 
tion in the window material, the gas filling used, 
the active length, and the wavelength of the radia- 
tion. The efficiency is maximal at a certain wave- 
length; it falls toward short wavelengths because 
the radiation becomes more penetrating, and to- 
ward long wavelengths on account of increasing ab- 
sorption in the window and in the dead space [12]. 


Thin windows of materials of low absorption 
are used to increase the efficiency at long wave- 
lengths: 10 » of mica (absorption coefficient for 
CuKa about 40 em7), 20 u of aluminum (40 em”), 
and vacuum-tight beryllium 0.1 mm thick (1.6 cm). 
The dead space is also minimized. The efficiency 
in the short-wave region is raised by filling the 
counter with heavy gases under increased pressure 
and by lengthening the counter. The strongly absorb- 
ing NalI(T1) is used in optimal thickness in scintilla- 
tion counters [22]. 

Figure 1 gives efficiency curves. An argon- 
filled Geiger counter and a xenon-filled proportion - 
al counter give efficiencies that fall rapidly toward 
short wavelengths; the efficiency for 1.54 A (Cuk@) 
is 60% but is only13% for 0.71 A (MoK@). This be- 
havior produces an effect of partial filtration rela- 
tive to films, especially if there is no absorption 
edge in the short-wave region of the continuum to 
delay the fall in efficiency. A krypton-filled pro- 
portional counter here gives much worse perform- 
ance on CuK@, on account of its absorption edge at 
0.87 A, but it can be used with MoK@ (0.71 A) atan 


efficiency of some 30%. The scintillation counter 
gives nearly 100% efficiency in the range 5 to 75 
keV, with some slight fall at long wavelengths due 
to absorption in the window. 

Dead Time 7.! This determines the line- 
ar range of a Geiger counter, which has T= 200 
usec and so loses about 10% of the counts at 500 
sec-!, Proportional and scintillation counters have 
dead times less than 1 usec, so the practical dead 
time is governed by the circuits and may be around 
a microsecond. 

Energy Resolution. A Geiger counter 
produces a pulse whose height is independent of the 
energy deposited in the counter, whereas propor- 
tional and scintillation counters produce pulses pro- 
portional to the energy deposited. Discriminators 
can then be used to eliminate x-ray background [12- 
15]. The energy resolution of a proportional counter 
is 2-2.5 times better than that of a scintillation 
counter, but the latter is still adequate for many 
purposes (Table 1). 

Electronic Systems. These employ 
standard circuits derived from nuclear physics; Fig. 
2 shows the block diagram of the system commonly 
employed with a Geiger counter in recording x rays. 

The other counters need two additional units: a 
linear amplifier and a single-channel analyzer. The 
EHT supply must also be more stable, and the cir- 
cuits must have a dead time of about 1 usec [23-25]. 
The system counts over a period set by the timer 
and also records the mean rate over the averaging 
time of the integrator. 


J VD ineartowOt ae Dili tPaccome ver: 


A counter records the rays only over a narrow 
angular range at any one time; counter and speci- 


1The time following the recording of a quantum for which the coun- 
ter will not respond to a further quantum, 
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Recorder 


HT supply oe 
| register 
Fig. 2. Block diagram of a Geiger-counter system. 


men must move, the angles being read with a goni- 
ometer. 

Bragg-Brentano [7, 8,26] or Zeeman-Bolin [27] 
focusing methods are used with powders to obtain 
adequate intensity in conjunction with good angular 
resolution. The first type is simpler in design, 
while the second enables one to use a broad diverg- 
ent primary beam and is free from the difficulty of 
the Bragg system that all the reflections occur in 
planes parallel to the surface of the specimen. The 
Zeeman-Bolin method is thus preferable for phase 
analysis of specimens showing a preferred orienta- 
tion. 

The time distribution of the measurements im- 
plies that the primary beam must be of high stabil- 
ity. The two methods of stabilization are: (1) cur- 
rent and voltage stabilization at the x-ray tube [8, 
26,28]; (2) a monitor counter working from the 
primary beam, which controls the tube [29,30]. The 
monitor counter can be used without tube control if 
it records an intensity proportional to that of the 
primary beam (e.g., from an aluminum foil placed 
in that beam), the other readings being taken as 
ratios to those from the monitor [31,32]. 


3. Intensity Measurement and 
Quantitative Phase Analysis 


The diffractometer has several differences from 
photographic methods in intensity measurement.The 
first arises from the random time distribution of 
the x-ray quanta, which produces a relative error 
inversely proportional to the square root of the num- 
ber of quanta recorded. The second is that focusing 
methods and flat specimens must be used because 
the intensity is proportional to the number of crys- 
tals in the reflecting position, while the relative er- 
ror in the intensity is inversely proportional to the 
square root of that number. The third is that the 
finite resolving time of the counter gives rise to ad- 
ditional errors in the intensity measurement. 


REVIEW 


Random Counting Error. The standard 
deviation o of the measured intensity is VN if the 
pulses are uniformly distributed, with N =nT, nbe- 
ing the count rate and T the counting time; the rela- 
tive standard deviation n is 1/VN. A background 
Np causes the measured intensity to be the sum of 
the diffraction intensity and the background: 


N=N,+p- 
Then 7 for the diffraction intensity Ny alone is 


,  _VNTNb. KK 4) 
VFN Np KK — 1) Vat” 


(3.1) 


where K = n/ np. 

A line equal in intensity to the background (K = 
2) requires three times the period needed to attain 
a given level of error in the absence of background 
(Fig.3). The 7 applicable to a ratemeter employed 
for the same purpose with the counter at rest is [33] 


1 
"Vn hee, 


‘ 


(3.2) 


and so the effective counting time is double the time- 
constant of the RC circuit. 

There are two cases to be considered in auto- 
matic recording: 

A. RC >> 7, in which case the loss of charge 
during the time spent measuring the line is negli- 
gible, so the relative error in the measurement of 
the peak intensity is 


pi. 4 
Vnre 


B. tT > RC, which is intermediate between the 
first case and the fixed counter case. Here, 


4 (3.3) 


(3.4) 


The error differs very little from that for a fixed 
counter under ordinary conditions: by 9% for RC = 
0.57, and by 4% for RC = 0.25T. The probability 
of a positive deviation of 40 p (in which og is the 
probable deviation) is 0.0035, so it is best to check 
deviations from the mean background intensity in 
excess of 40 p. Only deviations of width equal to 
that of the line, namely,0.67Vn/VT +2RC, should 
be checked, so the deviations to be checked are 
those equal to the line width and of height 

4+ 0.67Vn/VT +2RC [34]. High accuracy and sensi- 
tivity may be attained by increasing the counting 
time and taking readings on the pattern point by 
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Fig. 3, Relation of random counting error to back- 
ground rate; K is the ratio of total rate to back- 
ground rate. 


point [81], which can be done automatically [35, 36]. 
Another way of reducing the error is to minimize 
the background. 

The background has the following components: 
(L) inherent counter background, (2) white radiation 
scattered by the specimen, (3) Compton scattering 
in specimen and air, (4) fluorescence in the speci- 
men, (5) coherent scattering in the air, and (6) co- 
herent scattering in the amorphous component of 
the specimen. 

The first arises from cosmic rays and radio- 
active materials in the counter; it is minimized by 
shielding and by reducing the working or sensitive 
volume (all types of counter). The second is mini- 
mized by the use of a smoothed high-voltage supply 
to the x-ray tube and can be eliminated completely 
if a monochromator is used [37]. Fluorescence may 
be minimized by proper choice of working wave- 
length; in addition, the Compton, fluorescent, and 
white components may all be eliminated by placing 
a monochromator between specimen and counter 
[38]. 

The first four components lie at wavelengths 
other than the working wavelength, so they can be 
suppressed completely with discriminators. Coher- 
ent scattering in the air or amorphous binder can- 
not be eliminated by any energy-discrimination me- 
thod, so particular attention has to be given to other 
techniques. Photographic methods employ a long 
collimator (which comes close to the specimen) and 
a trap for the beam, as well as evacuated cameras; 
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a collimator may also be placed in front of the coun- 
ter, and the volume traversed by the primary and 
diffracted rays may be evacuated [31]. 

The specimen should not contain a large amount 
of amorphous binder. The best materials are xyl- 
ene, 5% Canada balsam in xylene, glycerol, or 
(best of all) no binder at all [39-41]. 

Most diffractometers employ a system that al- 
lows the use of a bent-crystal monochromator of 
symmetrical [81,32,41] or unsymmetrical [43] type. 
The discriminator and crystal monochromator may 
be further assisted by the use of fluorescence to 
produce monochromatic radiation [44]. Differential 
filters have been used to isolate the monochromatic 
part of the diffracted beam [9, 45]. 

Intensity Fluctuations from Varia- 
tion in Number of Reflecting Particles. 
The crystals in a polycrystalline specimen may be 
randomly oriented (no preferred orientation), where- 
upon the relative standard deviation of the number 
in the reflecting position (and,hence,of the intensity) 
is given [46] by 


8ruv, R® sin 0 


os V sew a R80), a)? 


where pu is the absorption coefficient and Vg is the ef- 
fective volume of a crystallite: 


(3.5) 


NM 
Ve= DV hi, 


jf; is the proportion of particles having volume vj, 
Rg is the radius of the goniometer (distance from 
the principal axis to the counter slit), A is the 
cross section of the primary beam in a plane 
through the principal axis, m is the repeat of the 
plane, Wy is the transverse dimension of the pro- 
jection of the focal spot, Hf is the lengthwise di- 
mension of that projection, 6@ is thehalf-width of 
the diffraction line, and Hz is the height of the coun- 
ter slit. 

Figure 4 gives the relative standard error as 
a function of the effective crystallite size and ab- 
sorption coefficient. 

Errors of the order of one percent require 
specimens crushed to a grain size of less than 5 u 
[47]. Another method of reducing the error by a 
factor of 2-3 is to rotate the specimen in its own 
plane [48,49]. 

Lost Counts. These are caused by the fin- 
ite resolving time and represent one of the major 
obstacles to obtaining accurate intensities. 

Two pulses separated by less than the resolv- 
ing time will be counted as one. The pulses have a 
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Fig. 4. Relative standard error in the intensity as af- 
fected by crystallite size and absorption coefficient p. 


uniform time distribution, so 


fone 

Tees (3.6) 
Here N' is the true number of counts, Nj is the 
number observed, and Ty is the resolving time 
(closely equal to the dead time). 

There are several methods of measuring re- 
solving times: the foils method [50-52], multiple- 
source methods [53], with an oscilloscope [54, 55], 
and by adjusting the circuit resolving time [9]. The 
oscilloscope method is the most accurate and reli- 
able. Additional difficulties arise if the strength of 
the x-ray source varies; here, the loss is depend- 
ent on the exact time course of the intensity [4,6,9]. 

Major difficulties arise in quantitative phase 
analysis if the material has particles that readily 
take up a preferred orientation; great care must 
then be taken in specimen preparation, which must 
have no tendency to give rise to a preferred orien- 
tation [41]. It is also possible in this case to em- 
ploy dilution with an amorphous phase, the relative 
intensity being extrapolated to zero content of that 
phase [56]. 

Use of an internal standard in phase analysis 
has reduced the relative error in the assay of a@- 
quartz in industrial dust to about 2% when the quartz 
content is about 50% [57]. Residual austenite in 


steels has been determined in much the same way[4]. 


Black's work [58] on bauxites (evaluation of new 
deposits) gives a good illustration of the use of quan- 
titative phase analysis. Here, use was made of one 
principal line from each of the minerals usually 
present in specimens: gibbsite, bohmite, kaolinite, 
hematite, goethite, and quartz. 
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The peak intensity alone was measured for each 
line, the readings being taken automatically with a 
modified Norelco diffractometer. The operator 
merely sets the counter in the required position, 
readings being taken over periods of 30 sec on line 
and background, the system then returning to the 
initial position. The percent content of each mineral 
(here hematite and goethite are taken as one) is de- 
duced from a set of curves relating the line and 
background intensities to the mineral content. The 
results may be passed to a computer to give the 
chemical composition. The relative standard devia- 
tion for assay of H,O and Al,O; is about 1.5%, but 
rather more for other components. The mean-square 
deviation from chemical analyses is 1-2% for all 
components. 

Phase analyses of polymers and macromolecu- 
lar substances generally are difficult by photograph- 
ic methods, so the diffractometer has become wide- 
ly used here [59-62]. The holder for the polymer 
must allow recording of the radial and azimuthal 
distributions of the scattered radiation as well as 
deformation of the specimen. A thermostat is used 
to cover a range of temperatures [59]. The primary 
beam is rendered monochromatic by reflection from 
a flat crystal of pentaerythritol. A Soller collimator 
produces a wide parallel beam. The stability is pro- 
vided by a ratio system based on a second counter. 
This apparatus has been used to measure the pro- 
portions of crystalline and amorphous phases at 
various temperatures in polyisobutene, natural rub- 
ber, polystyrene, polyethylene, and other such sub- 
stances [60,61]. The proportion of amorphous phase 
is indicated by the general scattering at 9 between 
0 and 60° relative to the integral scattering for an 
entirely amorphous specimen (one at a high tem- 
perature). However, it has been shown [62] that the 
proportion of amorphous phase during the crystal- 
lization of elastomers can be deduced quite accu- 
rately from the main amorphous ring. 

Focused reflections have some disadvantages 
for organic materials. The irradiated area in- 
creases very greatly at small 6, which causes a 
marked shift in the peak intensity and unsymmetri- 
cal line broadening. 

Transmission methods are free from these dis- 
advantages; adequate intensity and resolution are 
obtained from broad divergent beams and a compo- 
site divergent slit in front of the counter [63]. This 
method has been used in research on the @-keratins 
[64]. 

Good resolution with broad beams in transmis- 
sion may also be obtained with a focusing mono- 
chromator that collects the scattered rays [42]. 
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The proportion of crystalline phase in an or- 
ganic material that can be obtained in completely 
crystalline form may be deduced from the ratio of 
some strong crystal reflection to the same for a 
speeimen that is 100% crystalline. This method has 
been applied to specimens of sodium penicillin [65]. 


4. Effects of Instrumental Factors on 
Line Shapes and Crystallite-Size 
Determination 


Instrumental factors tend to alter the shape of 
a peak that otherwise is determined by the size dis- 
tribution of the crystallites: the peak as a whole 
shifts from the Bragg position, the shape becomes 
unsymmetrical, and the peak is broadened. 

Only the shift is important in qualitative phase 
analysis, for it hinders derivation of interplanar 
distances, whereas the broadening is basic to deri- 
vation of the size distribution; but all three effects 
are important in accurate measurements of cell pa- 
rameters, especially the displacement from the 
Bragg position, which impairs identification of the 
center. 

The final line shape is obtained by superimpos- 
ing the following functions [66, 67]. 

1. The diffraction function f (QA@), which is de- 
pendent on the size distribution of the crystals; A@ 
is the deviation from the Bragg angle. 

2. A function g,; dependent on the spectrum of 
the source. 


3. A function g, dependent on the source shape. 


4, A function g3 representing defocusing re- 
sulting from displacement of the plane of the speci- 
men from the focal circle. 

5. A function gy governed by the vertical di- 
vergence of the primary and scattered beams. 

6. A function g; dependent on the penetration 
of the primary beam into the specimen. 

7. A function gg dependent on the counter slit 
width. 


Figure 5 shows the instrumental functions and 
the final shape. 

The calculated shape agrees well with the ob- 
served one. The diffraction and instrumental widths 
are additive only for large 9; a correction curve 
has normally to be used. The high resolving power 
of the diffractometer extends the limit of crystallite 
size determination to 3000 A. The following are 
the factors that broaden and displace the diffraction 
line. 

Source Spectrum [68]. The form of g; 
is dependent on the line shape of the characteristic 
radiation and is given by 


1 
(A6)2-+ 02 

in which K; and ©, are coefficients dependent on 0 
as well as on the half-widthand wavelength. Any 
continuous component is symmetrical with respect 
to the Bragg position and is appreciable only at 
large @. 

Source Shape [2].A Gauss function gives 
a good approximation to the actual distribution along 
the short axis of the projection of the focal spot, and 
SO g) is given by 


ds ars Z: 


(4.1) 


—Ke (A0)? 
’ 


(4.2) 


coe——ae. 


in which Ky is a constant related to the half-width of 
the projection. 

This half-widthis about 0.02° in some diffrac- 
tometers employing projection of the focal spot at 
an angle of 3°. 

Defocusing [33]. The deviation from the 
Bragg position arising from deviation of the surface 
from the focal plane is 


x? sin 20 


» 3 4.3 
2R8, (4.3) 


Here, x3 is the distance from the center of the 
specimen and xsmax = bp/2 is the half-width of the 
Specimen. 

The maximal deviation is given by 


Af = — 


b%, sin 20 


24 ie 4.4 
Sky (4.4) 


The beam may not cover the specimen, in which 
case 


AG, max = — 


4R? 4 4.5 
in? (4.5) 


in which bg is the width of the section of the primary 
beam by a plane through the principal axis of the 
zoniometer, The shape of the 6 function broadened 
by defocusing is given by 


i ARTES BEEN 
83 (A6) dA0 V AO s (4.6) 
Vertical Divergence of Beams [69, 
70]. Here, A@ is governed by the height Xf of the 


focal point, the height Xp of the specimen, and the 
height x, of the counter slit, all above the principal 
plane of the goniometer: 


8 
ai poe 
A ee cane CLES (4.7) 
(tp—4%,)/Rg=%, (4.8) 
(2, —2#,)/Rg= 8, (4.9) 


586 REVIEW 


-00 0 +16 -40 7) 140 200 -4@ 040 -40 O 
a b c d e 


0 “40 


-40 0+40 -40 
f g 


Fig. 5. Functions governing diffraction peak shape. a) Spectrum of source; b) source 
shape; c) defocusing; d) vertical divergence; e) absorption; f) counter slit width; g) 


final shape. 


(4.10) 


If the first and last heights are small relative 
to the second (xf << Xp» Xz << Xp)> 


= —B = tp/R 
2 4, che 0 (4.11) 
Aj = — E 
9 R2 
aR 
Ti, ¢ 
AG ae ee 8k (4.12) 
8R*g 


If the first two are small relative to the third 
(xf S594 Xp << Xy); 


a=0, P=2,/ Kg. 


2 ctg 20 
KG ot As Se 
aS (4.13) 
H2, ctg 20 
AO{ max = Sees 3 
16 Ry (4.14) 


Thus, in the particular cases we have 


4 
§s(A9) ~ 7% - (4.15) 
A@ may be positive or negative in the general case; 
the maximal negative values are given by (4.7)- (4.9) 
with 

HL, is EG 
Ly = a F J P => = 7: > 


and are put as 


Ao; max 
2(1p~p+H,) (A, +H 
(Hp +H)? + (Hz +H py + UE Le P) 
= cos 20 ete 26 
16 R? hake 
8 (4.16) 


The maximum positive values are given by 


0 (AO4 max) 


=, 
Ox» 


and are 


Ags max = 
32R5 


ctg 26. 
(4.17) 


Calculations agree closely with experiment [33]. 
Penetration of Beam [71]. Here, A@ 
is related to the depth of the reflecting layer xs by 


(eee cos § ‘ 
Re (4.18) 


The maximal deviation is 


___ ad, cos 8 
Maman ake a ; (4.19) 
in which dy) = x;max is the specimen thickness. This 
penetration broadens the 6 function in accordance 


with 


25 (A0) ~ S28 (e2h-%s cosec 0) ~ e440. (4.20) 
Finite Counter Slit Width. This trans- 
forms a 6 function to a rectangular function. 
5 naekGecision) Measuremenitots Cend 


Parameters 


It will be clear from the above that the follow- 
ing factors that cause asymmetric broadening will 
displace the peak: (1) defocusing (from flat form of 
specimen), (2) vertical divergence, (3) penetration. 

In addition, the goniometer must be correctly 
adjusted and must be properly set relative to the 
focus of the x-ray tube. 

Defocusing shifts the center of a 6 function by 


A0Omax 
§ g A@d(A6) 
Ao. = 0 __ AO; max 
0 A@émax a 3 c 
gd (A®) 
J (5.1) 


Figure 6 shows the two cases arising for a tri- 
angular line. 


REVIEW 


46, 
‘ LOnqy >4 4 


Fig. 6. Displacement of the peak of a triangular line of half- 
width A due to unsymmetrical broadening. 


The peak is shifted by A@,;max/4 in the first 
case and by \/3 in the second. Alternatively, the 
shift might be deduced by reference to the line join- 
ing the middle of the horizontal intercepts, in which 
case it is A93;max/ 8, and so coincides with the 
shift in the center of gravity. 

It is difficult to deduce the shift due to vertical 
divergence for the general case; but for the cases 
of spot and slit heights small relative to specimen 
height, or spot and specimen heights small relative 
to slit height, we have, respectively, 


2 
8R? 
8 
and 
2 
A6zZ max = — i Neal ; 
4R5 


which are analogous to (4.4) and imply a shift in the 
peak of 


AO, max 
3 


A Soller collimator gives a shift of 


o? ctg 2 6 

12 ? 
in which o is the distance between adjacent plates 
divided by the length of these. The shift due to 
penetration into a specimen of thickness t is given 
by (4.20) and (5.1) as 
sin 20 d, cos. 9 
Auk, aad Roe’ cosec 6 = ys 


Displacement of the surface from the axis by a dis- 
tance s causes the peak to move through an angle 


Ales (5.2) 
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The over-all shift is then 
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The errors arising in the interplanar distances 
from these factors may [4, 71,72] be eliminated? in 
a cos?@ extrapolation (Fig.7). The counter zero 
may be determined precisely in several ways, the 
screen method [73] being the most accurate, as this 
gives the setting to 0.001°. 

The parameters of cubic lattices have been de- 
termined for polycrystalline specimens in absolute 
terms to within the error (+ 0.004%) with which the 
wavelength of CuKe@ is known by applying these me- 
thods of extrapolation and zero setting; the relative 
accuracy was substantially higher, especially when 
thermostatic control and point-by-point readings 
were used, 29 being found to 0.001° [72]. 

Also, 9 in the range 78-88° may be measured 
[74] by using a counter with a ring slit. 

A Laue goniometer may also be used to ex- 
amine lattice symmetry and distortion. 

A proportional counter working into a multi- 
channel analyzer can overcome any uncertainty in 
wavelength measurement and so can give exact lat- 
tice parameters [75]. 


6. Intensity Measurement on Single 
Crystals 


Ionization methods for integral intensity here 
are of high accuracy and sensitivity; the time re- 
quired is also reduced. The last feature is very 
important, on account of advances in direct sign de- 
termination and electron-density calculation by com- 
puter. For instance, the intensities of 2485 reflec- ~ 
tions were measured in 55 h with a diffractometer, 
of which 24 h were accounted for by accessory 
operations [76]. 


16° 


Fig. 7, Relation of line shift to @. 


*This is not sc for the shift due to vertical divergence, though it 


may be that more exact correction for this divergence is responsible 
for the higher accuracy of lattice parameters obtained by extrapo- 
lation, 
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Fig. 8. The equatorial method of 
measuring single-crystal intensities, 
a) Plane of primary beam and prin- 
cipal axis; b) equatorial plane. 


Photographic methods for integral intensity re- 
cord for reciprocal-lattice planes separately: one 
coordinate must thus be known, the other two co- 
ordinates being deduced from the position of the 
spot on the pattern. Beams reflected at different 
angles meet on their way to the film in layer photo- 
graphy, whereas the crystal must be coupled to the 
counter in ionization methods, for which purpose 
all three coordinates must be known in the general 
case. The reflections may be recorded in sequence 
with a’counter, instead of layer by layer. 

Types of layer recording have been described 
with the counter set up on a Weisenberg instrument 
to employ the plane angle method [77] or any type of 
layer recording [78]: equal-inclination, perpendicu- 
lar beam, plane angle, and constant cone. 

A goniometer head rotating around the counter 
axis enables one to record only one reciprocal-lat- 
tice plane without remounting the crystal, this pass- 
ing through the origin node and being perpendicular 
to that axis [79,80]. Simultaneous motion of the 
counter around an axis in the equatorial plane en- 


ables one to use the perpendicular-beam method [81]. 


A special goniometer for use with the equal-inclina- 
tion method provides variation of the angle between 
the primary beam and the common axis of crystal 
and counter [82]. Polar recording requires an addi- 
tional axis of rotation for the crystal in the equa- 
torial plane and parallel to the closest-packed rows 
in the reciprocal lattice [71, 83, 84]. 

This axis in the equatorial method enables one 
to displace the line hmkly (Fig. 8) through an angle 
y in the equatorial plane. The nodes in this series 
are then brought to the reflecting position by rotat- 
ing the crystal through an angle w around the prin- 
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cipal axis of the goniometer together with the ap- 
propriate rotation of the counter around the same 
axis through an angle Y, as in the recording of the 
zero-layer plane. This method is readily used with 
a simple attachment to the goniometer of a powder 
diffractometer [84]. This extra axis introduced into 
an equal-inclination system enables one to follow a 
nodal line parallel to this axis by rotation through 
an angle @ to the position at which it touches the re- 
flection sphere (Fig.9). The equal-inclination me- 
thod (u = —v) applied to nodal series so placed 
merely requires w (counter angle) and Y relatedby 
Y = 2w, which impliesa simple kinematic coupling. 
The latter feature greatly reduces the number of 
special setting operations in this method (called 
tangential). Other possibilities appear if the attach- 
ment to the powder goniometer contains a further 
axis @ (Fig.10) that itself may rotate around an- 
other axis X of the attachment, this axis being per- 
pendicular to ® and lying in the equatorial plane 
(these two additional axes form an Euler system in 
conjunction with the principal axis 2 of the goniom- 
eter). This attachment allows one to record re- 
ciprocal-lattice nodes by the cone method [76], in 
which the crystal is set with one axis (e.g.,b*) 
along the ® axis; the X axis is set perpendicular to 
the plane that bisects 26, and the ® axis is inclined 
by an angle x by rotation around the X axis. Then 
the crystal is turned through ¢ around the ® axis to 
bring into the reflecting position in the equatorial 
plane all nodes lying on a cone with the b* axis as 
its axis and a vertex angle of x, the radius vector 
corresponding to @. 

It is usual to measure the integral intensity in 
one of two ways: with a wide slit or with a converg- 
ent beam. The counter is kept fixed in the first me- 
thod, the slit being opened to accept the entire re- 
flected beam and the crystal being oscillated over a 
narrow range around the reflecting position. Crys- 
tal and counter are stationary in the second method, 
which takes less time; but it requires a uniformly 
emitting spot of appropriate size [9,85]. 

A recent description [82] concerns a diffract- 
ometer that records the intensity, crystal angle, 
and counter angle for a single crystal, this infor- 
mation being sufficient to give the F* solid. 

This instrument employs equal inclination, the 
counter moving through half a slit width during a 
complete rotation of the crystal, which provides all 
combinations of the two angles and so brings all 
nodes of the layer plane into the reflecting position. 

A system of relays is used to slow the motion 
of crystal and counter on reflections and to give’ 
high speeds on background measurements. This 
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Fig. 9. The tangential method of 
measuring the integral intensity 
from a single crystal. a) Plane of 
primary beam and principal axis; 
b) plane P. 


saves much time, and the device thus provides ra- 
pid construction of reciprocal-lattice planes. 

Another automatic system employs a computer 
to derive the setting angles of counter and crystal 
from the known cell parameters. 

A servo system employing selsyns would ap- 
pear best for setting these angles, as in a neutron- 
diffraction goniometer [86]. 


73 Textures 


The accuracy and rapidity of counter methods 
as applied to pole figures makes these essentially 
better than photographic ones here; such a figure 
can be recorded in about an hour for a given plane 
with a powder diffractometer. Texture studies 
merely require some minor modifications to allow 
rotation of the specimen around two mutually per- 


pendicular axes. The intensities provided by the re- 
quired plane are measured with the counter fixed at 


the appropriate 8, which makes the absorption and 
scattering volume corrections much simpler than 
those for the photographic method. There are two 
mutually complementary methods in the construc~ 
tion of pole figures: transmission and reflection. 
The first gives the outer part of the figure, while 
the second gives the central part. 

The principles of the transmission method are 
detailed in [87]. A thin specimen is mounted in a 


holder providing rotation around an axis perpendicu- 
lar to its surface and around an axis coincident with 
the goniometer axis. The counter is fixed at the 
Bragg angle for the appropriate plane (Fig.11), 
while the specimen rotates around the goniometer 
axis. The intensity corresponding to a diameter of 
the pole figure is recorded. The process is then re- 
peated with the specimen turned through an angle 
around the D axis, which gives the intensity distri- 
bution along a diameter at an angle 6 to the previ- 
ous one, etc. 

This method avoids the construction of Wulff 
nets; it needs only polar nets. The change in the po- 
sition of the specimen relative to the beam requires 
correction for the geometrical path of the beam in 
the specimen [88]. The formula given in the original 
paper is correct for rotation angles @ up to 40-50°. 
This method does not give the central part of the 
figure, because the beam is cut off by the holder in 
this range. 

The above scheme has been used with slight 
changes in the slit system to give a marked simpli- 
fication in the calculations [89]. Broad divergent 
beams and a narrow counter slit are used, so only 
the central part of the diffraction peak is recorded. 
The diffracted intensity is then independent of the 
angle of rotation over a fairly wide range (deviation 
from the Bragg angle by up to 40°) for specimens of 
a certain thickness; the absorption and scattering- 
volume corrections are then eliminated. 

The complete pole figure may be recorded, and 
absorption corrections eliminated, by replacing the 
rotation around an axis perpendicular to the plane 


Fig. 10. The cone method for single crystals. a) Equatorial 
plane; b) O(b*)* section. 
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Counter 


Fig. 11. Method of turning specimen in tex- 
ture examination by transmission. 


by examination of a series of cylindrical specimens 
[90], the axis of each lying in the rolling plane,each 
having a different angle to the rolling direction. The 
scattered intensity is measured in a plane perpen- 
dicular to the axis of the specimen. The rotation of 
the specimen is coupled to that of the recorder 
chart. 

An automatic device [91] has been described 
for use with thin specimens in transmission. The 
specimen is rotated continuously around axis D 
(Fig.11), which is synchronously coupled to the ro- 
tation of the polar net in the recorder. The speci- 
men is rotated in steps of 5° around the axis of the 
goniometer, this being coupled to the radial posi- 
tion of the pen in the recorder. The correction for 
the geometrical path of the rays is performed auto- 
matically. The intensity scale has seven steps, the 
pen having accordingly seven colors. Contours can 


then be drawn in directly via lines of the same color. 


The central part of the figure for a massive 
specimen may be recorded in reflection by using 
the same motion of the specimen as in the method 
of Decker et al. (Fig.12a). Here, the path correc- 
tion must be used. The method has the serious dis- 
advantage that the rotation angle @ cannot be larger 
than @ for the given plane, because only a narrow 
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Fig. 12. Methods of rotating the specimen in texture studies in 


reflection, 
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central part of the pole figure is recorded [92]. A 
more generally applicable method was described by 
Schulz [93]. The vertical slit system in the spec- 
trometer is replaced by a set of horizontal slits that 
produce a beam parallel in the vertical plane but di- 
vergent in the horizontal plane. The specimen ro- 
tates around axis D and also around the line of in- 
tersection of its plane with the principal plane of 
the goniometer (C in Fig. 12b). 

The specimen and counter are set to the posi- 
tion corresponding to the interplanar distance. Ro- 
tation on the C axis corresponds to recording along 
a diameter of the figure. The absorption is such 
that the intensity should not be dependent on ®, but 
the intensity is actually found to fall at large , on 
account of inevitable defocusing during rotation on 
C, which causes broadening of the peak [94]. This 
effect may be reduced by widening the counter slit, 
at the expense of loss of resolution. Satisfactory 
results are obtained up to = 40° under these con- 
ditions. 

An automatic specimen holder has been made 
for use in Schulz's method, which overcomes the 
following problems: exact and rapid setting of the 
specimen, increasing the surface area involved in 
the reflection for specimens with rather large crys- 
tals, and automatic intensity recording for speci- 
fied positions of the specimen [95, 96]. 

Schulz's method forms the basis of a goniom- 
eter for examining preferred orientation [97]. The 
drive rotates the specimen continuously around the 
D axis (Fig.12b), one turn around this correspond- 
ing to 10° around the C axis. The normals to the 
chosen plane then describe spiral paths on the 
stereographic projection. The intensity is recorded 
on a synchronized chart, from which it is trans- 
ferred to the pole figure. 

The accuracy may be improved by placing one 
spiral in another to give the turns separations of 5°, 
PA} (OWOe 

Schulz's method in reflection gives good results 
for the central part (up to 40°), while the method of 
Decker et al. [98] gives good results for the outer 
region. 

The complete pole figure may also be con- 
structed by the method of Decker et al. with only 
one holder, which is used in transmission and in re- 
flection. Good overlap between the two regions is 
provided by recording a higher-order reflection 
from the same plane [99]. 

Precision intensity recording has allowed the 
diffractometer to be applied to orientation studies of 
crystallites in cellulose, textile fibers, etc. [100]. 
The holder for use with these [101] enables one to 
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record the azimuthal and radial intensity distribu- 
tions for specimens varying in degree of elongation. 


8. Other Applications 


Quenching (the usual method for examining the 
kinetics of high-temperature phase transitions) is 
not suitable for many materials. The diffractometer 
does not need long exposures and so allows one to 
follow phase transformations at the transition tem- 
perature [102-104]. The main difficulty with high- 
temperature x-ray cameras (protection of the film) 
is absent here, and diffractometers can be used with 
furnaces going up to 1600°C or more [105]. The 
heater may be tungsten wire (1200°C) [91], platinum 
wire (1500°C) [106], or tantalum foil (1600°C) [105]. 
The furnace must be evacuated or filled with inert 
gas, to avoid oxidation or decarburization of the 
specimen. A stream of hot inert gas [107] may be 
used at 600-900°C, the outer jacket of the furnace, 
together with any external heaters and screens, hav- 
ing slots to pass the beams. The top and bottom 
plates are water-cooled, to protect the goniometer 
mechanisms. Specimen temperatures around 600°C 
may be produced simply by heating the holder [110]. 
Examples of this type of use are the isothermal 
decomposition of austenite in steel containing 2.5% 
Ni [103], structure determination for ®-CaO-2Si0O, 
[106], and ordering in CuZn and Cu3Au crystals at 
high temperatures [108]. 

The holder may be of insulating plastic for 
work at low temperatures [109,111]; the metal parts 
leading to the goniometer are fitted with heater 
spirals. The low temperature is produced by direct 
contact with liquid nitrogen [111,112] or by a jet of 
cold gas [109]. 

Single crystals may be oriented rapidly [113], 
the counter being set at 29 and the crystal being 
turned about two mutually perpendicular axes to lo- 
cate the appropriate plane. Two such planes are 
used with a stereographic projection to find the co- 
ordinates of a third plane, which is set up in order 
to check the diffractometer settings. 

Highly radioactive specimens represent a spe- 
cial problem; the detector should respond only to 
the characteristic x rays used. There are two lines 
of attack here: use of a double-crystal system [114, 
115] and differential discrimination with scintilla- 
tion or proportional counters [15,21]. The second 
crystal acts as a focusing monochromator (the first 
crystal is the specimen); it is set up between spe- 
cimen and counter to direct only the characteristic 
radiation to the latter. A piece of lead is placed be- 
tween specimen and counter. A modification that 
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avoids moving the counter with its heavy shield is 
to rotate the x-ray tube and specimen around the 
goniometer axis, but this method generally has the 
disadvantages of great loss of intensity at the se- 
cond crystal and complexity in the goniometer. The 
discriminator methods are free from these; stand- 
ard goniometers may be used. Scintillation coun- 
ters are of lower energy resolution than gas coun- 
ters, but they both provide recording free from 
background caused by the specimen,with loss of only 
15-20% of the characteristic radiation. 


* OK OK 


It will be clear rom this survey that the dif- 
fractometer should find steadily increasing use in 
many forms of structure analysis, as in quantita- 
tive phase analysis, line profile examination, cell 
parameter determination, texture studies, intensity 
measurement for single crystals at high and low 
temperatures, and with radioactive specimens. 

These applications will be facilitated by: 

1. Provision of attachments to the URS-50I to 
adapt it for these applications; 

2. Design of specialized diffractometers, es- 
pecially for use with single crystals and textures. 
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The second edition of the International Tables 
for the Identification of Crystal Structures appeared 
in 1952. This consisted almost entirely of informa- 
tion on the space groups, but diagrams for all the 
cubic groups were omitted by the editors because 
the specialists they had consulted were of the 
opinion that those in the first edition were incom- 
plete but at the same time overcomplicated and ob- 
scure, which made them of little real value. 

This omission gave rise to objections from 
many cryStallographers, and Pauling spoke for 
these in Acta Crystallographica when he stated that, 


if we rule out the incorrect reference of the reader 
(for the cubic groups) for part of the data to one or 
two of the simpler groups, the diagrams ofthe first 
edition were satisfactory (in which I agree with 
him); and further, the diagrams for the space 
groups are of value in instruction as well as for re- 
ference purposes, for without them we cannot have 
a sufficiently detailed study of the space groups as 
the basis of structure analysis. 

The editors agreed with these objections and 
promised to publish the diagrams for the cubic 
groups as a Special supplement to the first volume; 


Fig. 1. Lonsdale and Henry's draft of a diagram of the symmetry elements in 
group Of = Fd8c; the top left quadrant of the projection on a cube face is shown. 
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but this promise has remained unfulfilled although 
four years have elapsed, because the editors have 
not yet been able to decide on the best method of 
representing the cubic diagrams. 

‘Only in August of this year did I receive from 
Professor Lonsdale (chairman of the editorial com- 
mittee) a letter enclosing a diagram for one of the 
most complicated holohedral cubic groups (Fig. 1), 
with a request for comments to be sent to the 
editors (Kathleen Lonsdale and N.F.M. Henry) on 
the method of representation together (if possible) 
with my own suggestions. The letter concluded with 
a request that I bring these draft diagrams to the 
notice of our crystallographers by publishing them 
in this journal together with my reply to Professor 
Lonsdale and my proposed diagram (Fig.2) for this 
group. This diagram has been drawn up largely in 
accordance with the method of representation chosen 
by the British editors. 

It will be seen that the most novel feature is 
that stereographic projections are used at the 
points where the largest numbers of symmetry ele- 
ments coincide, and these the crystallographer 
may most readily employ to envisage the orienta- 
tion of the intersecting symmetry elements. Ex- 
cerpts from these projections even show the indi- 
vidual elements such as inclined triad and diad axes 
(rotation or screw). 

The editors also consider that their diagrams 
have the great advantage of a marked reduction in 
the number of symbols needed; in particular, they 
consider it possible to avoid altogether the 15 sym- 
bols for (mainly) inclined elements that were used 
solely for the cubic groups in the first edition (but 
it is not clear why they considered this a disadvan- 
tage). 

The following is my reply to Professor Lons- 
dale, with the introductory and final sentences de- 
leted. 

I have made a careful study of your proposed 
diagram (of the symmetry elements) for group Of = 
Fd3c and have some comments on the proposed new 
method of representation. 


1. The replacement of an insignificant square 
with a single diagonal by a stereographic projec- 
tion appears very effective. There is no doubt that 
such projections are especially valuable in space 
groups for crystallographers of orthodox back- 
ground in mineralogical crystallography, but Icon- 
sider that the method has not been carried to its 
logical conclusion in the proposed draft, and also 
(and in part because) this has disadvantageous ef- 
fects for other details, as I show below. 


2. How do the signs for inversion axes atlevel 
vf appear at the poles x and x (at y and y on the 
other edge) on the projection representing the inter- 
secting symmetry elements at y/, of the edge at 
level 0? I do not see why at these poles we do not 
have the signs for the 4, and 43; screw axes emerg- 
ing at level 0; and then the signs for the inversion 
axes (at level 4) would have to be displaced to the 
margin. Conversely, the stereographic projection 
around the origin of the cell must be supplemented 
with signs for the inversion or diad axes mentioned 
in the next point. The marginal signs are best given 
as in the International Tables of 1935. 

3. Do we denote the emergence of a horizon- 
tal axis by the sign for an inversion axis or by that 
for a diad axis? The proposed diagram nowhere 
distinguishes the special points of (horizontal) in- 
version axes, and we are forced to envisage them 
as occurring exactly in the middle between two 
points of emergence (from the poles of the stereo- 
graphic projection) of the inversion axis. This is 
incorrect as regards method, on account of the 
special role of these special points, which is so 
important that the origin of the cell is taken atsuch 
a point in many groups. 

I consider that it is absolutely necessary touse 
the mirrors (used in the 19385 edition) that cut 
across the diad axes and so indicate these points, 
and that these should be used in that form. In par- 
ticular, Iam unable to see how you can otherwise 
denote these special points in group Te = 148d (p. 
334 of the 1935 edition) or Of? = Ia3d. 

Mirrors with the heights stated must be shown 
for vertical inversion axes as well as horizontal 
ones. If these points are not overlapped from above 
by signs for other axes (e.g., Screw ones) in the 
projection, the method of symbolizing them can 
give rise to no doubt; if this occurs, the mirrors 
are shown by the side. This displacement of the 
mirrors, if these lie within the stereographic pro- 
jection, cannot give rise to doubt as to the position 
of the points. If mirrors are used, we should de- 
note the emergence of an inversion axis from a 
pole of a stereographic projection by a simple ar- 


row. This is shown to be correct by the next point. 


4, It is usual in group Of, = Fd38c to locate the 


origin not at a center of symmetry (where the sym- 
metry is 6) but at point T= 28, where the sym- 
metry is higher (namely, 12); but this is in no way 
reflected in the stereographic projection around the 
origin, from whose poles emerge neither diad axes 
nor inversion ones, whose existence must be de- 
duced by resort to another stereographic projec- 
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Fig. 2. Draft diagram for space group OF proposed here, 


tion, namely, that for point 4, on an edge. Diad or 
inversion axes essentially must be shown in the 
stereographic projection around the origin in order 
not to confuse nonspecialists, especially those (e.g. 
students) largely unfamiliar with the space groups. 

5. In this context the polyhedra of points in 
general positions (given in the 1935 tables) are not 
unnecessary, as you state in your explanatory text 
for the diagram. They identify directly the com- 
plex of symmetry T= 23 at the 16 corresponding 
points of the elementary cube. Of course, it is 
hardly an advantage of the diagram for oF that one 
can see only half of these polyhedra, and that all of 
the points of these have commas. I consider that 
all that is necessary (leaving out reference to pre- 
vious pages) is to show the polyhedra at levels 0 
and v/ (half with commas and half without) and to 
mention that they overlap at levels 44 and % with 
enantiomorphic polyhedra, which would be absolute- 
ly clear. 

6. Your symbols for triad and diad axes (rota- 
tion or screw) are evidently excerpts from the 
stereographic projections where they do not appear 
in the symmetry complex shown on the complete 
stereographic projection; they are larger for the 
triad axes but very small for the diad ones. This 
is good for the triad axes but not good for the diads. 


The method used in the 1935 tables is better and 
clearer. 

My diagram for group Oh = Fd3c has been 
drawn up in accordance with the above remarks,but 
it conforms generally to your basic principle; in it 
there is an illogical feature. I have not denoted a 
horizontal inversion axis by a distinct sign at all, 
but have everywhere inserted the symbol for a tri- 
ad inversion axis (triangle with blank circle). I 
recognize that this is incorrect, since even the be- 
ginner would find it difficult to imagine to which 
closest triad axis we assign a center of symmetry 
of height 4 or *%. It would be more logical to de- 
lete the blank circles from the signs for the triad 
axes altogether. 

7. Ialso raise the question of the representa- 
tion of d planes, for which I have proposed! a rep- 
resentation without arrows, though this may not be 
very convenient for graphic reproduction. Why not 
represent diamond planes by two bars and two 
points, to distinguish them from ordinary inclined 
planes ? 

8. The symmetry complex T = 23 is used in 
the diagram to distinguish (a) positions without pa- 
rameters; centers of symmetry are used to show 


*Dokl. Akad. Nauk SSSR, 59, 4, 701-702 (1948). 
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(c) positions of symmetry C3j = 3: and mirrors to 
show (d) positions (also without parameters) with 
S,= 4. We are left with no symbols for the no less 
important (b) positions with symmetry D3 = 32; it 
would be very appropriate to indicate these (above 


999 


and below the centers of symmetry) by a sign rep- 
resenting a quartz crystal, though I accept that this 
is only a pious hope. I regret that I am unable to 
discuss all the cubic groups in detail, because this 
would undoubtedly raise other points. 
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